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From Extreme Values of I.I.D. Random Fields 

to Extreme Eigenvalues of Finite-volume Anderson 

Hamiltonian 


A. Astrauska^il 


The aim of this paper is to study asymptotic geometric properties almost surely or/and 
in probability of extreme order statistics of an i.i.d. random field (potential) indexed 
by sites of multidimensional lattice cube, the volume of which unboundedly increases. 
We discuss the following topics: (I) high level exceedances, in particular, clustering of 
exceedances; (II) decay rate of spacings in comparison with increasing rate of extreme 
order statistics; (III) minimum of spacings of successive order statistics; (IV) asymptotic 
behavior of values neighboring to extremes and so on. The conditions of the results are 
formulated in terms of regular variation (RV) of the cumulative hazard function and its 
inverse. A relationship between RV classes of the present paper as well as their links to 
the well-known RV classes (including domains of attraction of max-stable distributions) 
are discussed. 

The asymptotic behavior of functionals (I)-(IV) determines the asymptotic structure 
of the top eigenvalues and the corresponding eigenfunctions of the large-volume discrete 
Schrodinger operators with an i.i.d. potential (Anderson Hamiltonian). Thus, another 
aim of the present paper is to review and comment a recent progress on the extreme 
value theory for eigenvalues of random Schrodinger operators as well as to provide a clear 
and rigorous understanding of the relationship between the top eigenvalues and extreme 
values of i.i.d. random potentials. We also discuss their links to the long-time intermit¬ 
tent behavior of the parabolic problems associated with the Anderson Hamiltonian via 
spectral representation of solutions. 
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1. IISITRODUCTIOISI 


1.1. Extremes of i.i.d. random fields 

In this paper, we assume that ^{x), x & If, are independent identically distributed 
(i.i.d.) random variables on a probability space (H, J", P), indexed by sites of the v- 
dimensional integer lattice If ^ with a distribution function P(^(0) ^ t) =: 1 — 
t G R; here Q denotes the cumulative hazard function of distribution. Define V = 
[—n;n]'^ fl Z'', the cubes in Z''. Let \V\ denote the number of sites in V. We write 
\x\ = |x*| for the lattice /^-distance between x = {x^,..., x") G Z’' and 0 G Z''. 

We consider the variational series (order statistics) 

Cl,v • C('^l,v) ^ ^2,V ■ C('^2,v) ^ ^ ^|V1,V ■ ( 1 - 1 ) 

based on the sample ■= {^( 2 ;): x G V}; here V = {zk^v- 1 ^ fc ^ 1^1}- The first 
\V\^ (0 < e < 1) terms of the variational series (1.1) are referred to as ^y-extremes 
or ^y-peaks. The coordinate Zk^v G V stands for a location of the fcth extreme value 
of ^y; 1 < fc ^ |D|. 

In this paper, letting \V\ —>■ 00 , we study the asymptotic geometric properties of 
^y-extremes almost surely and/or in probability. We are interested in the following 
functionals of order statistics (1.1): 

(EX) Exceedances of the sample ^y over high levels Ly, in particular, clustering 
of exceedances (Theorem 3.1). 

(SP) The decay rate of the spacings y - Ck+i,v and C[|vi=],v “ ^[\v\e],v 
parison with increasing rate of for fixed natural iG G N and 0 ^ e < 0 < 1 
(Theorems 4.3-4.7). 

(MIN) Minimum of the spacings 4,v ~ Ci+i,v) 1 ^ ^ ^ each 0 < £ < 1 

(Theorems 4.8 and 4.9). 

(NEI) Cv -values neighboring to fv -extremes, in particular, ^(z, + y) for 1 ^ ^ 

\V\^ and for fixed y ^ 0 (Lemma 5.1 and Theorems 5.3, 5.4). 

The conditions of the asymptotic results for (EX), (SP), (MIN) and (NEI) are given 
in terms of regular variation (RV) of the inverse function of Q- In Appendix A and 
Section 6, we discuss a relationship between RV classes of the present paper as well as 
their links to the well-known RV classes including domains of attraction of max-stable 
distributions. 

The asymptotic results for (EX), (SP), (MIN) and (NEI) and related RV classes 
were announced without the proof in (Astrauskas 2007; 2008; 2012; 2013). In this 
survey, these results are given in the most general setting with the detailed proof; 
therefore, they present self-contained topics of probability theory and may be consid¬ 
ered of independent interest. 


1.2. Extreme value theory for eigenvalues of large-volume Anderson Hamilto¬ 
nians 

Let us consider the finite-volume Schrodinger operators Hv = kAv -b fv on l‘^{V) 
with periodic boundary conditions (Anderson Hamiltonian); here k > 0 is a diffusion 


3 



constant; Atp{x) := X]|y-a:|=i ^(j/) is the lattice Laplacian, and the i.i.d. random field 
'■= {^( 2 ^)' ^ G V} is the multiplication operator (potential). Denote by y 
the Kth largest eigenvalue of the operators Hy, and let X^^y) {x £ V) be the 
corresponding eigenfunction normalized to have unit P-norm, vY = 1- 

Another aim of this paper is to show in what manner the asymptotic behavior of func¬ 
tionals (EX), (SP), (MIN) and (NEI) determines the asymptotic structure of the top 
eigenvalues ^ and the corresponding eigenfunctions, as V ^ If and AT > 1 fixed. 
In Section 2, we give an overview of rigorous statements on this relationship which are 
proved in the papers by Gartner and Molchanov (1998) and Astrauskas (2007; 2008; 
2012). In Section 6, we review and comment results on the asymptotic expansion 
formulas and Poisson limit theorems for the largest eigenvalues y as well as local¬ 
ization properties of the corresponding eigenfunctions. These results are proved by 
Astrauskas and Molchanov (1992), Gartner and Molchanov (1998), Astrauskas (2007; 
2008; 2012; 2013), Germinet and Klopp (2013), Biskup and Kdnig (2016) and other 
mathematicians. These papers are complemented by the present survey on the asymp¬ 
totic geometric properties of ^y-extremes and related RV classes of distributions. We 
here give proof sketches of the results on the extreme value theory for eigenvalues 
Xk,v (Sections 2 and 6) demonstrating their connections to asymptotic properties of 
^y-extremes. 

Thus, in this survey, we discuss in detail the following important branches of prob¬ 
ability theory: (i) extreme value theory for eigenvalues of the Anderson Hamiltonian 
Tiy = kAy + ^y which is a particular model of random matrices (Sections 1.4, 2 
and 6); (ii) asymptotic geometric properties of random i.i.d. fields (Sections 3-5) and 
(iii) regular variation of distribution functions (Appendix A). On the other hand, we 
briefly comment the links of the extreme value theory for eigenvalues to the follow¬ 
ing important topics of statistical physics: (iv) Anderson localization for the ran¬ 
dom Schrodinger operators "H = kA + ^l-) in the whole lattice Z'^ (Section 1.3), and 
(v) long-time intermittent behavior of solutions u of the parabolic problems associated 
with the Anderson Hamiltonian (PAM) via spectral representation of u (Section 7). 

Of course, asymptotic results for the Anderson models (time-dependent or -indepen¬ 
dent) are heavily necessitated by the asymptotic structure of high ^(•)-values, which 
in turn is determined by conditions on the regularity and tail decay of the distribution 
P(^(0) > t) = e”*^*^*) at its right endpoint tq := esssup ^(0). In the present survey, 
we focus on the case of unbounded from above i.i.d. potentials, i.e. tq := 00 , with dis¬ 
tributional tails heavier than double exponential, i.e., P(^(0) > t) = exp{—e°^*^} as 
t —>■ 00 ; cf. Sections 1.2, 2.2-2.3 and 6.1. For such distributions satisfying additional 
RV and continuity conditions, we will show that with probability one the ^y-peaks 
are spatially separated and differ in height as V t therefore, the top eigenvalue 
A^ y of the operator Hy = nAy -|- ^y is approximated by an isolated ^y-peak, say 
y), plus some corrections of order o(l) involving neighboring ^y-values. More¬ 
over, the ATth eigenfunction ■;/)(•; XK,y) is asymptotically delta like function at the site 
^t(,k) V S the localization center. (In this case, we will say that the ATth eigenvalue 
is associated with the site y, viz. A^^y ££ y). Therefore, we are able to 

apply the standard extreme value theory to prove Poisson limit theorems for the nor¬ 
malized extreme eigenvalues and their localization centers. From these Poisson limit 
theorems one obtains the limiting joint (max-stable) distribution for the normalized 
largest eigenvalues and their spacings, limiting uniform distribution for the normalized 
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localization centers and other important limiting distributions for eigenvalue statis¬ 
tics (Section 6.1). Eigenvalue statistics in turn play a crucial role in studying the 
intermittent behavior of the parabolic Anderson model, PAM (Section 7). 

For the lighter upper tails including the double exponential P(C(0) > t) = exp{— e*} 
and bounded tails i.e. tg < oo, we will prove the rough asymptotic expansion formulas 
for the largest eigenvalues. For such distributional tails, it will turn out that all 
^v-extremes are of comparable amplitude; therefore, the ibth largest eigenvalue is 
associated with a large island of higher ^y-values of a particular preferred shape, 
rather than an isolated ^y-peak. 

To illustrate the relationship between ^y-extremes and the largest eigenvalues Xk,v 
(V t more precisely, we now formulate Propositions 1.1-1.3 which are ’’typical” 
examples of the statements given in Sections 2-6. The first proposition tells us that, 
if the peaks of deterministic (nonrandom) functions fv =■ fvi ‘) are extremely sharp 
and widely spaced, then the ATth largest eigenvalue is approximated by the Kth 
largest value of ^y with sufficiently small error. 

Proposition 1.1 (see Theorem 2.2(ii) in Section 2.2). Fix constants K G N and 
0 < 9 < 1/2, and assume that the deterministic functions fy satisfy the following 
conditions as V f IF: 

min ^;+i_y(^, ^,) —>■ oo (distinct height of peaks), (1-2) 


--;—- min \z.,,—z ,, I —>■ oo (sparseness of peaks) (1.3) 

and, finally, 

< const (61) (1.4) 

for some 0 < const (9) < 1 (negligibility of the lower peaks). Then 
\,V = ^I,v +Oil/fi^v) for all IsflifK. 


We now give an example of i.i.d. random field £,{■) with sufficiently “heavy tails” 
possessing extremes like those in Proposition 1.1. 

Proposition 1.2 (see Theorem 4.3(i) with p = 1, Theorem 3.1 with R = 0 and 
Theorem 4.5). If has the Weihull distribution 

P(^(0) >t)= = e-‘“ [t ^ 0) (1.5) 

with a < 2, then the i.i.d. sample fv (V f satisfies (1.2)-(1.4) with probability 
1 -I- o(l). 

Since, by Propositions 1.1 and 1.2, the eigenvalues are very close to ^ as F f 
, it turns out that Poisson limit theorems (and the corresponding renormalization 
constants) for the largest eigenvalues are the same as that for ^y-extremes according 
to the following proposition. 
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Proposition 1.3 (see Theorem 6.9 and (Astrauskas 2012)). Assume that Qit) = 
t°‘ with a < 2, and write by '■= (log|y|)^/“. Define the point process My on 
[-1/2; 1/2]“^ x^hy 


\v\ 

My ^ (fc) with (k) 


^k,V ^k,V \ 

a-^b\r" J’ 


where 6x denotes the Dirac measure at X G [—l/2;l/2]'^ x R. Then My converges 
weakly to the Poisson process on [—1/2; 1/2]“^ x R with the intensity measure dx x 
e“‘dt. 

Moreover, for fixed K 1, the eigenfunction Xj^ ^) is exponentially localised at 
the site z^ y: 


limsup 

y 


max 
xA^k,v ja; — Zj, 


log |V’(x; Ak,v) 

I log by 


^ -1 


( 1 . 6 ) 


in probability. 

According to Proposition 1.3 the ATth largest eigenvalue Xk,v is associated with 
the ATth largest value of y, viz., Xk,v ^ zxy- For the lighter tails, say, Weibull 
distributions (1.5) with a 2, the landscape of y gets “smoother”, in particular, 
(1.2) fails. Therefore, Xk,v is associated with a lower and “slightly supported” y- 
peak, viz., Xk,v ^ Zt{k),v, where for a > 3, the index t{K) = Ty{K) tends to infinity 
as |P| —>■ oo. This in turn implies that further terms in expansion for ^ become 
essential; see (2.22) and Examples 6.12-6.13. Let us distinguish three classes (J)- 
(JJJ) of light tailed distributions (i.e., universality classes), which ensure a different 
asymptotic behavior of the eigenvalues Xk,v- 

(J) Distribution tails heavier than the double exponential function. Assume that 


logQ(t) = o(t) 


(1.7) 


and Q satisfies additional regularity and continuity conditions as t ^ oo. This class is 
presented by Weibull distributions (1.5) for arbitrary a > 0 and those with fractional- 
double-exponential tails 

P(^(0) > t) = = exp{—e*^} {t ^ to) (1.8) 


for 7 < 1. For such distributions, ^y-extremes possess a strongly pronounced geometric 
structure which can be described as follows: 

For arbitrary sufficiently small constants 0 < £ < d, there exist constants ci > C 2 > 
0 and (large) C > 0 such that almost surely 

min y ^ ^ (distinct height of peaks), (1.9) 


min \z, y — z ,J > (sparseness of peaks) 

l^l<n^\V\» I CV n,V I I I \ 


and, finally. 


^[ivyi.v - ?[|v|e].v ^ ^ (negligibility of the lower peaks) 


( 1 . 10 ) 


( 1 . 11 ) 
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for each large V; see Theorem 4.8 with k = 0, Theorem 3.1 with i? = 0 and 
Theorem 4.6 with p = oo. By the standard finite-rank perturbation arguments 
in (Astrauskas and Molchanov 1992) and (Astrauskas 2008), these properties of 
yield that there is no resonance between -peaks in the Anderson model for 
large V; therefore, the eigenvalues associated with a block of peaks can be deter¬ 
mined by the local eigenvalues associated with separate peaks (i.e., “relevant single 
peak”approximation). More precisely, for fixed natural K and V 'I almost surely 
the eigenvalue Xk,v of Hv = kAv+^v is approximated by the principal (i.e., the first 
largest) eigenvalue of the “single peak”Hamiltonian kA^ -I- ^(-) -I- ^{zr{K),v)^z^^K) v 
where logr(Ar) = o(log |y|). Here ^(-) is the “noise” potential; the site 2T(if),y G ^ 
is a localization center of the A'th eigenfunction ip{-; ^,) of TLv- Thus, Poisson limit 

theorems for the eigenvalues Xk,v of Hy are reduced to those for the principal eigen¬ 
values of the “single peak” Hamiltonians, which in turn are expanded into certain 
(nonlinear) series in ^{x) {x G V); cf. formulas (2.20)-(2.22), Theorems 2.3, 6.2 and 
discussions in Section 6.4. We finally notice that the A'th eigenfunction iIj)-; Xj^ y) 
is exponentially well localized, i.e., there exist non-random constants (decay rates) 
C > 0 and 0 < My —>■ oo such that with probability one 

\i^{x-,Xj,y)\i^Cexp{-Mv\x-{x G V) (1.12) 

for all V large enough. Consequently, ip{-',Xi^y) is asymptotically delta-like func¬ 
tion at Zr(K).v (Astrauskas 2008; 2013). This refers to the correspondence Xk,v 
Zt{k).v- Under assumption (1.7), asymptotic expansion formulas, Poisson limit theo¬ 
rems and localization theorems for the largest eigenvalues are derived by Astrauskas 
and Molchanov (1992), Astrauskas (2007; 2008; 2012; 2013). See also Grenkova et 
al. (1983) and Grenkova et al. (1990) for the case of Weibull distribution (1.5) with 
a <2. 

(JJ) Distribution tails lighter than the double exponential function. Assume that 

t~^ \ogQ{t) ^ oo (1-13) 

and Q satisfies additional regularity conditions as t tends to tg (= the right endpoint of 
Q). This class of potentials contains the important case of f{-) which is bounded from 
above {tg < oo) and those with fractional-double-exponential tails (1.8) for 7 > 1. For 
such f{-), it turns out that ^y-peaks possess a weakly pronounced geometric structure. 
In particular, almost surely v ~ ^[|vi«] v ® 1^1 all 0 < £ < 0 < 1 , 

so that the height of all ^y-extremes is of the same order ^i^y -|-o(l) (see Theorem 4.6 
with p = 0 and Theorem 3.1(i) with arbitrary R ^ 1 and 9{-) = 0 = const). In 
this case, the eigenvalue (AT > 1 fixed) does not longer correspond to an isolated 
potential peak, but to a flat extremely large “relevant island” of high ^y-values. More 
precisely, almost surely the top eigenvalue Xj^ y of the Hamiltonian "Hy approaches 
the local principal eigenvalue of the Hamiltonian restricted to a random connected 
region C V with the following features: The diameter of Ay.^p^. unboundedly 

increases, and f{-) possesses in Ay.^pj values of the order ^i^y -|- o(l) as U f i-e-, 
relevant island of potential values. Moreover, the ATth eigenfunction is expected to 
be highly concentrated in the neighborhood of the region Ay.^p^. In this case, we will 
say that the Alth eigenvalue is associated with A^.^^p^, viz. Xk,v ^ ■'^y-opt- ^or more 
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explanations, see Theorem 6.14 and the proof of Theorem 2.6, where the second order 
expansion formula for ^ is obtained. 

For the Bernoulli i.i.d. random variables ^(x) with tq = 1 (which is the particular 
case of (1.13)), Bishop and Wehr (2012) derived a more accurate expansion formula 
for the principal eigenvalue Aj of the one-dimensional Hamiltonian 'Hy {v = 1). 
They have showed that A^ is associated with the longest consecutive sequence of 
sites X G V with ^(x) = 1, i.e., the “relevant island ”of ^y-extremes, the length 
of which unboundedly increases as V '[ Z. See also (Sznitman 1998) for similar 
asymptotic results in the case of spatially continuous Schrodinger operators with a 
bounded Poisson potential of obstacles. Cf. Section 6.2 below. 

To the best of our knowledge, Poisson limit theorems for the (unfolded) largest 
eigenvalues were proved only in the case y = 1 and bounded ^(0), provided the 
distribution 1— e~^ satisfies additional continuity and tail decay conditions (Germinet 
and Klopp 2013); see also Section 6.2 below. For y ^ 2 or general RV conditions on 
Q satisfying (1.13), the Poissonian convergence of the top eigenvalues still remains an 
open problem. 

(JJJ) Double exponential type tails. Finally, assume that t ^ log (5(1) tends to a 
positive finite constant p~^ as 1 —>■ oo, i.e., the double exponential tails 

P(e(0) >t)= e-'3(‘) = exp{- (1.14) 

satisfying additional RV and continuity conditions. This class of distributions presents 
the intermediate case between (J) and (JJ). For such e”*^, it turns out that all ^y- 
extremes are of comparable amplitude; i.e., almost surely C[|v|e] v ~^[|vi«] v ~ 

\V\ —>■ oo, for any 0 ^ e < d < 1. Therefore, with probability one the top eigenvalue 
A^ y {K ^ 1 fixed) of the Hamiltonian "Hy is approximated by the local principal 
eigenvalue of the Hamiltonian restricted to a random connected region C V 

of bounded diameter, where ^y possesses high values of the optimal shape (so that 
^K,v ^ Ay.^pj). The optimal shape of ^y-values in Ay.^p^ C V is specified by deter¬ 
ministic variational principles. These considerations are referred to as the “relevant 
island” approximation; see Theorems 2.7 and 6.19 for the second order expansion 
formulas for the first largest eigenvalue A^ y, which have been originally derived by 
Gartner and Molchanov (1998). Moreover, the Kih. eigenfunction 7/>(-; Xk,v) is highly 
concentrated in the neighborhood of the region Ay.^pj, as proved by Astrauskas (2008; 
2013) for p large enough, and by Biskup and Konig (2016) for arbitrary p; see also 
Section 6.3 of the present survey. 

Rigorous results on Poisson limit theorems and further localization properties for 
the largest eigenvalues (in the case of double exponential tails) have been proved by 
Astrauskas (2007; 2008; 2013) for p large enough, and by Biskup and Konig (2016) 
for arbitrary p; see also the review paper by Konig (2016) and Sections 6.3-6.4 of the 
present survey for the discussions on their results and the proofs. 

Let us finally summarize the above observations: As the upper tails of potential 
distribution get lighter, the ^y-extremes {V f Zf) get less expressed; therefore, the 
number of higher ^y-values contributing to the asymptotic amount of the top eigen¬ 
values gets larger and concentration properties of the corresponding eigenfunctions 
become weaker. On the other hand, the general theory of Anderson localization 
(cf. Section 1.3 below) suggests that almost surely the eigenfunctions associated with 



the upper spectral edge of the Hamiltonian "Hy decay exponentially like in (1.17), for 
arbitrary potential distribution satisfying certain continuity conditions. 


1.3. Relations to infinite-volume Anderson Hamiltonians 

(I) Anderson localization. The Anderson model on the whole lattice 7A is given by 
the Hamiltonian 

■H = KA-k^(-) 

acting on Here, as above, A is the lattice Laplacian, k > 0 is a diffusion 

constant, and ^{x) {x G 7A) are i.i.d. random variables with a common distribution 
function 1 — e”*^. This is a basic model of disordered quantum systems introduced to 
describe the regions of energy levels (spectrum) of the electron in the random potential 
modelling electrical conductance regimes of alloys, crystals with impurities and so on 
(Anderson 1958). The energy spectrum Spect {%) of the Hamiltonian % = kA + ) 

is almost surely nonrandom: 

Spect (H) = Spect (kA) + Spect (^(-)) = [—2yK; 2^^] -|- supp (1 — e”*^), 

where supp (F) is the support of a probability measure generated by the distribution 
function F, and denotes the algebraic sum of subsets of real line. Therefore, 
with probability one, the spectrum consists of spectral bands situated in the interval 
[Amin; Amax]j where Amin and Amax are respectively the infimum (i.e. bottom) and the 
supremum (i.e. upper edge) of the spectrum. The most important property of the 
Hamiltonian "H on is the presence of pure point spectrum in the neighborhood 

of edges of spectral bands for any k > 0 and any y ^ 1. In particular, there exist 
(nonrandom) real constants A^ = Ai(K, y, Q), Li < L 2 , such that with probability one 

the spectrum in (Amin; Ai) U (A 2 ; Amax) is dense purely point, 

say {Afc}, and the corresponding eigenfunctions decay exponentially: 

\tp{x-,X^)\^Ckexp{-M\x-zJ} {x&U') 

for some random Cfc > 0, M > 0 and z^. £ IF (the localization center), provided e“'5 
is Holder continuous and ^(0) has some finite statistical moments. Moreover, for small 
K or the one-dimensional case y = 1, the whole spectrum Spect ("H) is dense purely 
point with a complete set of eigenfunctions in ?^(Z'^) that decay exponentially with 
probability 1. This phenomenon is known as Anderson localization for disordered 
systems; see, e.g., (Frohlich and Spencer 1983; Frbhlich et al. 1985; Simon and Wolff 
1986; Carmona et al. 1987; Aizenman and Molchanov 1993; Aizenman et al. 2001) for 
the proof of the above assertions under various conditions on potential distributions. 
Recall that, for the periodic ^(•), all the spectrum Spect ("H) is absolutely continuous 
in arbitrary dimension y ^ 1; and this is quite a contrast to the Anderson local¬ 
ization in the case of random potential. See the monographs (Pastur and Figotin 
1992; Stolz 2011; Kirsch 2008) and references therein for more discussions on the 
subject. In the present survey as well as in (Astrauskas 2007; 2008; 2012; 2013) and 
(Biskup and Konig 2016), the phenomenon of Anderson localization is illustrated for 
the top eigenvalues and eigenfunctions of finite-volume models. We here emphasize 
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the relationship between asymptotic geometric properties of ^\/-extremes and local¬ 
ization properties of the leading eigenfunctions of the operator T-Lv regarding their 
localization strength, localization centers, etc., as t 

Let us discuss briefly the basic ideas and methods explored in the study of the 
Anderson localization phenomenon in the multidimensional case v ^ 1. Fix an open 
bounded interval / C M which is covered almost surely by the spectrum Spect ("H). 
The proof of Anderson localization in the spectral intervals / relies heavily on the study 
of the resolvents t/(A -|- is) := {X + ie — 'H)~^ and the corresponding Green functions 
Q{X + i£] X, y) := G{X + i£)6y{x) {x, y € Z’') in the complex domain e > 0, A G /, where 
i = Alternatively, their finite-volume versions f/v,(A -I- is) := (X + is — 'Hy)~^ 

and Gv y) •= Qv (•^ + (^) ( 2 ^ 1 2/ G s-re also explored. The main task 

here is to prove that, for A G /, the Green functions t/(A -I- is] x, y) or Gy (•^ + x, y) 
decay exponentially in \x — y\ uniformly in e > 0, provided I is chosen close to the 
spectral edge or the diffusion constant k. is small. For simplicity, assume throughout 
that the potential distribution has a bounded density p(-) with bounded support. 

Frohlich and Spencer (1983) developed the multiscale method. They constructed 
inductively a sequence of relevant cubes V, V f with the following properties: 
For any fixed A G / and V' f Z'', with high probability the Green function G^, (A -I- 
is\ X, y) decays exponentially in \x — y\ for all y G dV', the “boundary”of V, and all 
X € V' far away from dV', and this estimate holds uniformly in e ^ 0. By finite- 
rank perturbation formulas, this estimate implies an absence of absolutely continuous 
spectrum in I with probability one. Even the stronger form of the construction of 
the relevant cubes V' (“uniformity” in A G /) is applied to prove the exponential 
decay of the eigenfunctions ip{-',X) associated with (generalized) spectral values in I 
and, consequently, the presence of the pure point spectrum in I with probability one 
(Frohlich et al. 1985). See also the survey by Kirsch (2008) for a detailed discussion 
on the multiscale analysis. 

In the above considerations, one should apply the IFe^ner estimate. Loosely speak¬ 
ing, this estimate states that the mean number of eigenvalues of 'Hy in the interval 
/ does not exceed I/IIEIG, provided the density p{-) of potential distribution satisfies 
p{‘) ^ C. In particular, the latter guarantees the bound of probability to find at 
least one eigenvalue X^ y in a small spectral interval. The Wegner estimate and its 
modifications are the basic probabilistic tools in the proof of Anderson localization; 
see, e.g., (Kirsch 2008). 

Simon and Wolff (1986) applied the so-called spectral averaging methods in mul¬ 
tiscale analysis to obtain the effective condition for the Anderson localization in the 
interval I, for the distributional density as above. Recall the Simon- Wolff criterion: 
If for each x G and for Lebesgue-almost every A G / with probability one 

lim E \G{X + is-,x,y)f < 00 , (1-15) 

yez- 

then the operator 'H has only pure point spectrum in / with probability one. If in 
(1.15) one claims a square-summability of the Green function with weights for 

some nonrandom m > 0, then the corresponding eigenfunctions decay exponentially 
with probability one. 

Aizenman and Molchanov (1993) and Aizenman et al. (2001) developed the fractio¬ 
nal moment method to prove the Anderson localization in the spectral intervals I. The 
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task here is to obtain (and explore) the exponential decay of the averaged Green func¬ 
tions, where the average is taken over the random potential. This method enables 
to avoid a complicated dependency of the previous constructions on individual po¬ 
tential configurations in the almost sure setting. Under certain continuity conditions 
on potential distribution, the key statement in (Aizenman and Molchanov 1993) is 
the following fractional-moment criterion: If for fixed 0 < cr < 1, there are constants 
Cl > 0, Ml > 0 such that 

V.{\Q{Xie;x,y)\'^) ^ Cl for all x, y G , (1-16) 

for all A G / and uniformly in £ > 0, then one has, with probability one, the Anderson 
localization in the interval I for the operator "H. 

This implication can be proved by using the Simon-Wolff criterion like in (1.15). 
On the other hand, cr < 1 is chosen to depend only on continuity assumptions for 
potential distribution. In particular, the Holder continuity implies that the left-hand 
side of (1.16) is Hnite. 

Eq. (1.16) can be proved by using suitable Hnite rank perturbation arguments, i.e., 
Krein formulas. These formulas imply that the kernel -|- ie;x,y) is equal to a 
simple rational function in the variable ^(a;) with coefficients depending on potential 
values outside x. Now one can estimate the left-hand side of (1.16) as an integral of 
rational function. These estimates are shown to form a certain iteration procedure, 
from which one deduces (1.16). When applying the iteration scheme, a crucial fact 
is the assumption that the diffusion constant k is small or the interval / is near the 
upper edge of the spectrum. 

Applying similar arguments, Aizenman et al. (2001) deduced a fractional-moment 
finite-volume criteria for the Anderson localization in the interval I. Roughly speak¬ 
ing, these criteria state that, if for some 0 < cr < 1 and some V, the expectation 

E (|I/v(A-I-ie; 0, y)|'^) is sufficiently small 
for all y G dV, uniformly in (A;e) G / x ]R_|_, 

then the exponential decay (1.16) holds true in I. The finite-volume criteria and their 
implications are shown to hold under continuity assumptions on potential distribution 
mentioned above. 

The criteria of (Aizenman et al. 2001) also imply the exponential decay of eigen¬ 
functions of the operators Hv, associated with the eigenvalues in (compact) spectral 
intervals / of Anderson localization, in particular, for / at the upper edge of the spec¬ 
trum Spect {'H). Under the continuity conditions on potential distribution as above, 
we are able to formulate this result in a precise form: There are nonrandom constants 
c > 0 and M > 0 such that with probability one 

IV’(a^;A,,,v^)K |U|^exp{-M|a;-(x G U) (1.17) 

for some y G V (localization center), for all X^. y G I and all V large enough; 
cf. (Klopp 20il). See also the surveys (Hundertmark 2008; Stolz 2011) for detailed 
discussions on the fractional-moment methods and their applications. 

(II) Local fluctuations of eigenvalues in the spectral regions of Anderson localiza¬ 
tion. We denote by /pp := (a; h) an open interval of real axis such that a certain 
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fractional-moment finite-volume criterion is fulfilled in /pp. Therefore, the spectrum 
in the whole of /pp is purely point, so Spect (Ti.) fl /pp C Spect pp('H) with probability 
one (i.e., Anderson localization in /pp). The spectral intervals /pp are distinguished by 
the Poissonian asymptotic behavior of eigenvalues A, , close to a fixed A° G /pp, of the 
finite-volume model T-L^ as V f Z'^. Here and to the end of this section the potential 
distribution is again assumed to have a bounded smooth density with bounded sup¬ 
port. These limit theorems are formulated in terms of the integrated density of states, 
viz. N{X), and the density of states, viz. n(A) := iV'(A) (A G Recall that N{-) is 
the nonrandom distribution function of eigenvalues defined as the almost sure limit 
of the empirical distribution function A^y(A) := #{fc: A^. < A}/|H| as \V\ —>■ oo. 

Moreover, the Wegner estimate implies that N{-) is an absolutely continuous function 
with bounded density n(-), provided the potential distribution has a bounded density. 
The support of n(-) coincides almost surely with the spectrum Spect (H); therefore, 
N{\) 0 (resp., N{\) —>• 1) as A approaches the bottom (resp., the upper edge) of 

the spectrum. See, e.g., (Kirsch 2008) for the definition of functions N{-), n(-) and 
their properties. 

Now pick a number A° from the interval /pp such that n(A°) > 0. We consider the 
normalized eigenvalues 

A0_^:=|H|n(A0)(Afc,v-A0) (1 ^ fc ^ |H|), (1.18) 

and define the corresponding point process Ady on R by 

Al 

A<»v:=E«»!,v 

k=l 

With these assumptions and abbreviations, one needs to show that the point process 
A4y converges weakly, as H to the Poisson point process on R with intensity 

measure d A, i.e. the Lebesgue measure. The first result in this direction was proved 
by Molchanov (1981), who considered the one-dimensional spatially continuous ran¬ 
dom Schrodinger operators. In the case of the Anderson Hamiltonians in Z‘' with 
arbitrary 1, this Poisson limit theorem was established by Minami (1996). Killip 
and Nakano (2007) proved the Poisson convergence of both the normalized spectral 
values (1.18) and localization centers of the corresponding eigenfunctions, extending 
Minami’s result. The proof in the multidimensional case ly ^ 1 relies on applications 
of the fractional-moment criteria for the Anderson localization, combined with the 
Wegner and Minami estimates that control the probability of finding, respectively, at 
least one and two eigenvalues of Tdy in a small interval. In particular, these estimates 
imply the upper and lower bounds of order \V\~^ for the gap between successive 
eigenvalues close to a fixed A° as above. From Poisson limit theorems for (1.18), 
one can extract some information on the statistical properties of eigenvalues in the 
spectral regions of Anderson localization. For example, one can obtain the limiting 
distribution for the normalized spacings of eigenvalues, the limiting joint distribution 
for the normalized eigenvalue and its localization center, and other important limiting 
distributions for eigenvalue statistics. 

Recently, Germinet and Klopp (2013; 2014) presented a comprehensive study on 
limit theorems for eigenvalues of large-volume Hamiltonians in the spectral regions 
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Ipp of Anderson localization. For fixed A° G Ipp as above, the authors considered the 
unfolded eigenvalues 

Tly:=\V\{N{\^y)-N{\^)) (1.19) 

i.e., the eigenvalues under nonlinear renormalization. They proved the following limit 
theorems for eigenvalue statistics: 

1) Poisson limit theorem for the point process based on the unfolded eigenvalues 
(1.19), where the limiting Poisson process coincides with that for (1.18); 

2) Poisson limit theorem for the point process based on both the unfolded eigenval¬ 
ues (1.19) and the normalized localization centers of the corresponding eigenfunctions; 

3) Limit theorems for the empirical distribution function of the normalized spacings 
of eigenvalues close to A°; 

4) Limit theorems for the normalized distance between localization centers of the 
corresponding eigenfunctions; 

and other important limit theorems for various statistics related to the eigenvalues 
and eigenfunctions in the intervals /pp. Moreover, Germinet and Klopp (2013) consid¬ 
ered more general random Hamiltonians % with convolution-type long-range kinetic 
operators instead of the Laplacian: For such Hamiltonians, the Poissonian asymptotic 
results were shown to hold also for the eigenvalues at the spectral edges (in partic¬ 
ular, for the Kih largest eigenvalues with fixed K ^ 1). At the spectral edges of 
the Schrodinger operators, this result was shown to hold for the one-dimensional case 
V = 1] cf. Section 6.2 below. 

The proof of the limit theorems in (Germinet and Klopp 2013; 2014) relies heavily 
on the techniques of the general theory of Anderson localization, including applica¬ 
tions of the fractional-moment criteria in the large-volume setting, as well as various 
versions of the Wegner and Minami estimates. 

In the proof of the above results, the following crucial observation is related to 
localization properties of the corresponding eigenfunctions of the operator T-Ly (and 
is also quite close to the context of our paper treating the top eigenvalues and eigen¬ 
functions of 'Hv)'- Namely, the eigenvalues A^ near A° of the operator 'Hy can be 
approximated, with a “good” error, by independent local eigenvalues of the operators 
restricted to much smaller disjoint cubes V (z) C V centered at z. This approximation 
is feasible due to the facts that the localization centers of the corresponding eigen¬ 
functions of T-Lv are located far away from each other (so the probability of having 
at least two centers in cubes V{z) is asymptotically negligible), and the eigenvalues 
of T-Ly “live” on potential values in small spatial neighborhoods of the corresponding 
localization centers because of the exponential decay of eigenfunctions. One needs to 
consider the small spectral interval ly C Ipp, centered at A°, such that the number 
of eigenvalues in Iv unboundedly increases as K f 'E''. The latter is fulhlled 
if, say, \Iy \ X \V\~‘^ for some 0 < c < 1. The study of the structure of eigenvalues 
A^ inside ly is based on the Wegner and Minami estimates. In particular, using 
the eigenvalue approximation described above, one has with high probability that the 
number of eigenvalues in ly is roughly approximated by N{Iy)\V\ (large deviation 
principle); here iV(/) is the probability measure associated with the integrated density 
of states. Also, note that the above continuity conditions on potential distribution 
imply another important bound: 

Nilv) ^ const (1-20) 
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for some i9 ^ 0 and for all V large enough. These bounds are crucial in estimating 
the probabilities of the occurrence of a single or several eigenvalues in small spectral 
intervals for operators over cubes V and V{z) C V introduced above. 

In (Germinet and Klopp 2013), the results on local fluctuations for eigenvalues ^ 
are extended up to spectral edges. Here the proofs rely on the improved versions of 
Wegner and Minami estimates which ensure the more explicit control of eigenvalues 
in small spectral intervals ly, since the amount N{Iv) is now allowed to be 
exponentially small in instead of (1.20). Thus, this case includes situations, 

where the measure N{-) is extremely small, for example, at the edges of spectral bands 
where the phenomenon of “Lifshits tails” occurs. 

We finally notice that the Poisson limit theorems for the lower eigenvalues (1-18) or 
(1.19) (in the spectral regions Ipp of Anderson localization) agree with the results of 
the present paper and our earlier papers on the extreme value theory for the largest 
eigenvalues oiT-L^, with K ^ 1 fixed. However, in limit theorems for the largest 
eigenvalues, the choice of normalizing constants depends strongly on the regularity 
and tail decay conditions of the potential distribution, in contrast to limit theorems 
for eigenvalues near A*^ G Ipp, where the normalizing constants are simply expressed 
in terms of the (integrated) density of states. In particular, the spacings of the top 
eigenvalues have an asymptotic order, which is much larger than = the order 

of the gaps between successive eigenvalues close to A° G Ipp. It is worth mentioning 
that the proof of Poisson limit theorems for eigenvalues in Jpp relies heavily on the 
methods and ideas of the general theory of Anderson localization; and neither the 
extreme value theory nor links to the asymptotic geometric properties of random 
potential are explored. 

We finally mention the following open problems regarding the Anderson transition 
phenomenon for the infinite-volume Hamiltonians % = /tA-|-^(-): The first conjecture 
is that for y > 3, the spectral bands outside some neighborhood of the spectral edges 
consist of purely absolutely continuous spectrum and the corresponding eigenfunctions 
are delocalized. The second conjecture is that the eigenvalues of Hy in the spectral 
intervals of delocalization obey non-Poissonian asymptotic behavior as fo t i rather 
the limit theorems like in the theory of Wigner random matrices with light-tailed 
entries; cf. Section 1.4 below. The first problem is partially solved for the very special 
models on the Bethe lattice as well as for the Schrddinger operators with sparse 
potential (e.g., Kirsch 2008; Molchanov and Vainberg 1998, 2000). 


1.4. Relations to random matrices 


(I) Wigner random matrices. Another important model of disordered quantum sys¬ 
tems (in particular, heavy nuclei atoms) is presented by real symmetric random ma¬ 
trices 


Hn = {KJ 




with i.i.d. centered entries h. ^ (i ^ j) and N —>■ oo; i.e., large Wigner matrices (Mehta 
2004; Anderson et al. 2010). The extreme eigenvalues (i.e., high energy levels) A^ „ 
and the corresponding P-normalized eigenvectors of Hn are here interpreted as the 
basic states of quantum systems. 

Recently, there has been much progress toward the extreme value theory for the 
eigenvalues A^^ „ of Wigner matrices as N —>• oo and AT > 1 fixed. It has been 
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turned out that there are two different regimes of asymptotic behavior of the largest 
eigenvalues, depending on the tail decay rate of the entries in absolute value: 

(11) For polynomially decaying distributions 

I > t) = + o(l)) as t —>■ oo (1-21) 

with /3 < 4 (very heavy tails), Auffinger et al. (2009) proved that with high probability, 
the Kih largest eigenvalue X^.n of Wigner matrices Hm is approximately equal to 
the ATth largest value among \h^ {1 ^ i ^ j ^ N), K ^ 1 fixed. This in turn 

implies Poisson limit theorems for the normalized eigenvalues where the 

normalizing constants > 0 are chosen the same as in the corresponding limit 
theorems for extremes of ^. | (1 ^ ^ j ^ N). Recall that distributional tails 

(1.21) are in the domain of attraction of the max-stable Frechet law G/ 3 , therefore, 
the eigenvalue „ is of the order A'^ = A^^/^(const -|-o(l)); cf. Example 6.11 below. 
Moreover, with high probability, the Kth. eigenvector is asymptotically concentrated 
on two coordinates, i.e., it behaves like a superposition of two delta functions in limit 
as fV —>■ oo. See also (Soshnikov 2004) for the case (3 <2. To prove these assertions, 
one first observes that, under condition (1.21) with /3 < 4, the larger entries (in 
absolute value) are extremely sparse and strongly pronounced in comparison to other 
entries in Hjq. Thus, the standard perturbation theory for symmetric matrices is 
applied to conclude that the top eigenvalues and eigenvectors of the former matrix 
Hi^ are approximated by the corresponding eigenvalues and eigenvectors of a very 
sparse (symmetric N x N) matrix whose entries are the larger values of h^ ^’s in 
absolute value. This fact in turn enables us to apply the extreme value theory for the 
random variables \hij\ and, as a consequence, to prove Poisson limit theorems for the 
eigenvalues of Wigner matrices H^. 

(1 2 ) Assume that the \hij\’s have lighter tails (including (1.21) with /3 > 4 and 

Bernoulli entries), or the hij^s have some finite statistical moments of higher order 
and satisfy additional conditions on a distributional symmetry. Then the largest eigen¬ 
values A^ „ of Wigner matrices Hn are distinguished by non-Poissonian asymptotic 
behavior, rather the Tracy-Widom limit law; see, e.g., (Soshnikov 1999; Lee and Yin 
2014; Bourgade et al. 2014). In particular, the normalized top eigenvalues Xj^ j^/'/N 
tend almost surely to the nonrandom constant 2(]E/if i-e., the right endpoint 

of the support of their limiting spectral distribution density; cf. (Bai and Yin 1988). 
Moreover, the corresponding ^^-normalized eigenvectors are completely delocalized; 
i.e., with high probability their sup-norm does not exceed At“^/^(log . See, 

e.g., (Tao and Vu 2010; Erdos et al. 2013a; Vu and Wang 2015; Gotze et al. 2015), 
where this delocalization property is extended to all the eigenvectors of provided 
the tails of \hij\ are lighter than the exponential, or the ft-ij’s have a large enough 
number of moments. It is worth mentioning that limiting distributions for eigenvalues 
or eigenvectors (in particular, the Tracy-Widom limit law for the largest eigenvalues) 
can be explicitly compnted for Wigner matrices with Gaussian entries; see, e.g., (An¬ 
derson et al. 2010). Thus, the usual comparison methods (four moments theorem. 
Green function comparison method, etc.) can be used to extend the asymptotic re¬ 
sults for the Gaussian case to general Wigner matrices; e.g., (Tao and Vu 2014). 

The transition from Poisson limit theorems to Tracy-Widom asymptotics for the top 
eigenvalues of Wigner random matrices was discussed in detail by Biroli et al. (2007). 
The value /3 = 4 in (1.21) or, roughly speaking, the fourth statistical moment indi- 
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cates here the threshold separating these two different regimes of asymptotic behav¬ 
ior. Note that the Wigner matrix model with light-tailed entries reflects the global 
or mean-field interaction; thus, the asymptotic geometric properties of entries do not 
play any role in limit behavior of eigenvalues and eigenvectors. The latter is in a 
sharp contrast to the one-dimensional Anderson Hamiltonian in H C Z, which is a 
random band (tridiagonal) matrix reflecting local interaction: the diagonal elements 
are i.i.d. random variables and the deterministic off-diagonal elements are given by 
the Laplacian. In view of discussions on the Anderson model (Section 1.2 above), it 
turns out that the diagonal operator necessitates concentration properties of eigen¬ 
vectors; meanwhile, the Laplacian forces these properties to be less expressed (thus, 
the geometric features of the model play here a crucial role). 

(II) Random band matrices. Recently, there has been a considerable attention 

drawn to symmetric random band matrices of size N ^ oo, where the matrix 

entries vanish if \i — j\ exceeds W, and other entries (above the diagonal) are 
i.i.d. centered random variables; here 0 ^ ^ iV is a band width. Random band 

matrices are natural interpolations between Anderson Hamiltonians {n = 1) and 
Wigner matrices. For band models, the asymptotic behavior of the top eigenvalues 
and eigenvectors depends strongly on the growth rate of the band width W = Wm ^ 
oo as well. Benaych-Georges and Peche (2014) considered the random band matrices 
H^\ whose entries have polynomially decaying distributions (1.21) with arbitrary 
/3 > 0 and band width W = with 0 < /x ^ 1. They established that the band 
model exhibits a phase transition depending on /x and j3, with /3 = 2(1 -|- /i“^) as 
the threshold separating two different regimes of asymptotic behavior of the largest 
eigenvalues „ and the corresponding eigenvectors {K ^ 1 fixed): 

(III) Assume (1.21) and W = such that /? < 2(1 -|- ^“^), i.e., either the dis¬ 
tributional tails are sufficiently heavy or the band of matrix is sufficiently narrow. 
This case includes heavy-tailed Wigner matrices, i.e., /x = 1 and /3 < 4, considered 
in (Ii) above. For 0 < /x < 1 and (3 < 2(1 -|- the asymptotic results are similar 
to that in (R). I.e., with probability 1 -I- o(l), the Kth largest eigenvalue of 

is approximately equal to the Kth largest value of the sample (1 < i ^ N, 
0 ^ j — i ^ W). Therefore, Poisson limit theorems for the normalized eigenvalues 
hold true, where An = (const' -|-o(l)) (cf. Example 6.11 below), 

and the Kth eigenvector is asymptotically localized on two coordinates. The proof 
of these assertions is again heavily based on techniques of the extreme value theory, 
in particular, describing asymptotic geometric properties of entries of band matrices. 
The latter is combined with the perturbation theory for matrices to derive simple 
asymptotic formulas for the largest eigenvalues and eigenfunctions of . 

(IR) Assume the h..^s are symmetrically distributed with tails (1-21) and let 
W = with P > 2(1 -I- i.e., either the distributional tails are sufficiently 

light or the band of matrix is sufficiently wide. In this case, the band models 
posses the mean-field features, like in the light-tailed Wigner models considered in 
(R) above. Thus, each eigenvector associated with the upper spectral edge is asymp¬ 
totically delocalized in the sense that its P mass is more or less uniformly spread over 
N coordinates. Moreover, the extreme eigenvalues do not longer obey Poissonian 
asymptotic behavior; in particular, they tend to a nonrandom positive constant when 
divided by 


16 



Earlier, Sodin (2010) considered the band matrices whose entries are sym¬ 

metrically distributed with tails lighter than the Gaussian tails, including Bernoulli 
entries. He studied the transition from localization to delocalization for eigenvectors 
at the upper spectral edge (as well as the corresponding limit theorems for eigen¬ 
values) with varying degrees of strength and generality. In particular, the results of 
this paper suggest that the eigenvectors associated with the largest eigenvalues are 
delocalized, provided W = with /x > ^cr = 5/6. See also (Erdos et al. 2013b) 
for similar delocalization results for eigenvectors associated with the inner part of the 
spectrum Spect 

In view of these observations on the localization properties at the upper spectral 
edge, we distinguish two paradigmatic models in the theory of random matrices. First, 
the general theory of Anderson localization suggests that, with probability 1, the 
one-dimensional Schrodinger operators "H = kA + ^ (i-C., the particular model of the 
band matrices with W = 1) have exponentially localized eigenvectors at the spectral 
edges, provided the potential distribution is arbitrary satisfying very mild continuity 
conditions (Carmona et al. 1987). On the other hand, the large Wigner matrices Hn 
(i.e., the band matrices with W = N) have localized eigenvectors associated with 
the largest eigenvalues if only entries are very heavy-tailed like in (1.21) with /3 < 4; 
meanwhile, for the light tails, eigenvectors ot are typically delocalized. See also 
(Spencer 2011) for a discussion on Anderson-type models (M = 0(1)), band matrices 
(M = o{N)) and Wigner matrices (M = N). 

Finally, it is worth emphasizing that the above asymptotic results for eigenvalues 
remain valid for the corresponding complex Hermitian random matrices with i.i.d. en¬ 
tries, instead of the real symmetric matrices. 


1.5. The earlier literature on extremes of i.i.d. random fields 

As already mentioned, most statements of the present paper on the ^y-extremes and 
the corresponding RV classes were announced in (Astrauskas 2007; 2008; 2012; 2013). 
We now provide a brief overview of the earlier literature on the related asymptotic 
results for extreme order statistics of i.i.d. random sequences and fields. 

High-level exceedances consisting of single rare ^y-peaks were studied in (As¬ 
trauskas 2001). Related asymptotic results (in particular, the so-called longest head 
runs in coin tossing) for Bernoulli distributed i.i.d. random variables ^(x), a; G Z, 
were discussed, e.g., in (Binswanger and Embrechts 1994). 

In the case of exponentially distributed ?7(0), strong limit theorems for the spacings 
Vk.v ~ Vk+i.v fixed) were proved by Astrauskas (2006). Devroye (1982) derived 
strong and weak limit theorems for mini^x,^|y|(C^ ^ y) where C(0) is uniformly 

distributed. 

In the case of ^(0) with arbitrary distribution, strong asymptotic bounds for y 
are given in (Shorack and Wellner 1986) where K is fixed, and in (Deheuvels 1986) 
where K = Ky —>■ oo. Wellner (1978) derived strong asymptotic bounds for the 
uniform fcth order statistics v (thus, for rj^y) uniformly in fc ^ 1. 

For the Gaussian random fields {^(a;): x G If} with correlated values, Astrauskas 
(2003) studied some asymptotic geometric properties of ^y-extremes almost surely, 
in particular, high level exceedances and minimum of spacings. The geometry of 
high level excursion sets of smooth Gaussian random fields in was investigated in 
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the monograph by Adler and Taylor (2007). See also (Gartner et al. 2000) for some 
geometric aspects of high peaks of smooth Gaussian random potentials related to the 
long-time asymptotics for the spatially continuous parabolic Anderson models. 

For the extreme value theory for random variables, in particular, characterization 
of the domains of attraction of max-stable distributions, we refer to the monographs 
by Resnick (1987), de Haan and Ferreira (2006), Leadbetter et al. (1983), Embrechts 
et al. (1997). See also the monograph by Shorack and Wellner (1986) for a detailed 
account of strong and weak limit theorems for order statistics and their functions 
related to mathematical statistics. Finally, the monograph by Bingham et al. (1987) 
provides a detailed account of the theory of regularly varying functions. 

In the proof of a number of our statements on ^y-extremes, we explore the repre¬ 
sentation y where / := Q'^ is the generalized inverse function of Q and, 

as above, rj^ y stands for the kth extreme value among independent exponentially 
distributed random variables r](x) (x € V) with mean 1. Due to the nice properties 
of 77j, y (for instance, r/^. y is a sum of independent exponentially distributed random 
variables), we first obtain the asymptotic results for y , which are then transferred 
to ^f. y under appropriate conditions on /. These conditions are formulated in terms 
of regular variation (RV) of f(s) as s —>■ oo. We further give a characterization of 
RV classes, in particular, their links to continuity and tail decay of the distribution 
1— e~^ at the right endpoint; see Appendix A. They are also compared with the well- 
known RV classes including the domains of attraction of max-stable laws, 0-regular 
variation, asymptotically balanced, etc; see Appendix A. 

An interesting further problem is an extension of the present asymptotic results for 
functionals (EX), (SP), (MIN), (NEI) to other classes of random fields ^(-) including: 
1) independent non-identically distributed random variables; 2) random fields with 
correlated values, in particular, Gaussian fields (Astrauskas 2003) and moving average 
fields defined as a linear combination of i.i.d. random variables with nonrandom real 
coefficients. See, e.g., the review papers by Elgart et al. (2012), Tautenhahn and 
Veselic (2015) for a detailed background of the random alloy type models kA + ^(•) 
with the moving average potential ^(•). 

1.6. Notation. Representation of i.i.d. random fields 

Let us introduce the further notation and remarks we use throughout the paper. We 
denote by IR+ the positive half-axis and by N positive integers. Let log^ stand for 
the jtimes iterated natural logarithm. For real a, b, we write a V 6 := max(a, b) and 
a Ab := min(a, b), and [a] for the integer part of a. Given a subset U C lA , we write 
\U\ for the number of its elements. Let dist (C/, C/') stand for the lattice l^-distance 
between subsets U, U' dlA. The summation over x €V: a<|x|<6is abbreviated 
to X]a<|a;|<f)- ^ 0 : I^ 0 1: etc. we denote various large numbers, values of which 

may change from one appearance to the next. Similarly, const, const ' etc. stand for 
various positive constants. We write 1/0-1- = oo, log(O-l-) = —oo and l/oo = 0. Let 
g o h = g{h{-)) stand for a composition of real functions g and h. Also, for g > 0 and 
ft. > 0, we write g{t) x h{t) as t —>■ oo, if the ratio g{t)/h{t) is bounded away from 
zero and from above for all large t. 

By Qy (A; Qy ',x,y) {x GV, y £ V) we denote the Green function of the Hamiltonian 
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nAy + in l‘^{V), viz. 

g^{X;Cv;x,y) := g^,{X;Cy)6y{x) := {X-kA^ -C^)~^6y{x). 

Here Sy(- ) is the Kronecker delta function, i.e., Sy{x) := 1 ii x = y, and Sy{x) := 0 if 
x^y. 

Throughout the paper we suppose that all random variables are defined on a com¬ 
mon probability space (H, P). Let E stand for the expectation with respect to P. 

Recall Q(t) := — logP(5(0) > t) is the cumulative hazard function of an i.i.d. ran¬ 
dom field ^(x) = (x) (w G H; X G I/), and let tq denote its right endpoint 

tq := sup{t: Q{t) < oo}. Without loss of generality, we shall assume throughout 
that 0 < tg ^ oo. Clearly Q : {—oo;tq) —>• IR+ U {0} is a right-continuous nonde¬ 
creasing function such that Q{—oo) = 0 and Q{tq) = oo. Most of the conditions of 
our results are formulated in terms of the inverse of the cumulative hazard function 
defined by 

/(s) := := inf {t: Q{t) ^ s} (s G IR+) (1-22) 

(thus / : R+ —>■ {—oo]tq) is a left-continuous nondecreasing function such that f{s) 
tends to tq as s ^ oo). The reason for this is the following useful representation of 
order statisties : 

Ci,v := fiVi,v) > ^ 2 ,v ■= f{V 2 .v) > ■ • • ^ ?ivi.v := /(^ivi,v)> (1-23) 

where 

Vl,V ■ ^ V 2 ,v ■ v{^ 2 ,v') ^ ^ V\v\,v * ^('^|V|,v) (^'24) 

is the variational series based on the sample ijv ■= {ri{x): x G H} of exponential 
i.i.d. random variables with mean 1. 


1.7. Outline 

In Section 2, we collect conditions on deterministic functions fv in terms of func¬ 
tionals (EX), (SP), (MIN), and (NEI), which yield expansion formulas for the largest 
eigenvalues ^ of the discrete Schrodinger operator Hv = kAv + fv on l‘^{V) 
as E t Section 2.1 provides rough bounds for We then study Xj^ ^ in 

the cases of with extremely sharp peaks (Section 2.2), dominating single peaks 
(Section 2.3), dominating large islands of high ^y-values the diameter of which un¬ 
boundedly increases (Section 2.4) and, finally, dominating islands of high ^y-values 
the diameter of which is bounded (Section 2.5). The results of Sections 2.2-2.5 follow 
simply from the more general statements of (Astrauskas 2008; 2012) and Section 2.4 
in (Gartner and Molchanov 1998). 

Sections 3-5 contain the main results of the paper dealing with asymptotic behavior 
of extremes of the i.i.d. random field ) with the distribution function satisfying 
certain RV and continuity conditions at the right endpoint. Functionals (EX), (SP), 
(MIN) and (NEI) are studied in Sections 3, 4.1-4.2, 4.3 and 5, respectively. 

Section 6 provides an overview of current results on extreme value theory for the 
spectrum of the Anderson Hamiltonian TLv = nAy + V f ^ with an i.i.d. po¬ 
tential The issues under discussion include the asymptotic expansion formulas 
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and Poisson limit theorems for the largest eigenvalues and their localization centers. 
We consider separately three cases of the distribution tails e~^ of ^(0): the tails 
are heavier than the double exponential function (Section 6.1); the tails are lighter 
than the double exponential function (Section 6.2), and the double exponential tails 
(Section 6.3). As already mentioned, we give proof sketches of most theorems of this 
section demonstrating their connections to the results of Sections 3-5 on -extremes. 
In Section 6.4, we comment and compare the proofs of Poisson limit theorems stated 
in Sections 6.1 and 6.3 and proved in the earlier papers by Astrauskas and Molchanov 
(1992), Astrauskas (2007; 2008; 2012; 2013) and Biskup and Konig (2016). 

In Section 7, we discuss the long-time intermittent behavior of the solutions to the 
parabolic problems associated with the Anderson Hamiltonian TL = kA -|- ^(•). We 
focus on the representation of the solutions in the spectral terms of the operators 
Hv = wAy -I- ^y. In view of this representation, we discuss some techniques of 
the extreme value theory for eigenvalues of Hy, that can be applied to study the 
intermittency properties of time-dependent Anderson models. 

Finally, in Appendix A, we characterize and compare the RV classes of distributions 
introduced in Sections 3-6. 


2. ASYMPTOTIC EXPANSION FORMULAS FOR THE LARGEST EIGENVA¬ 
LUES OF DETERMINISTIC HAMILTONIANS 

Let V = [—n-,nY H (n G N) be a sequence of cubes. By introducing the periodic 
norm |x| := \x\n '■= ininj,g( 2 n-i-i)Z'' ~ y\^ V may be considered as a sequence of tori 
tending to Z‘'. We are interested in the finite-volume Schrodinger operators "Hy = 
s;Ay-|-^y on P(V) with periodic boundary conditions. Recall that k > 0 is a diffusion 
constant. Ay denotes the lattice Laplacian on Y(V) (i.e., a restriction of the operator 
Atpix) := J2\y-x\=i'^(y) to torus V) and ^y := {^y(x): x G R} G [-oo;oo)tol are 
deterministic functions, i.e., potential. The values —oo of ^y (i.e., “hard obstacles” ) 
are allowed to include the cases which are interesting from a physical point of view; 
see, e.g., (Biskup and Konig 2001; Konig 2016). Write 14 := {x G R: ^y(x) > —oo}. 
Then "Hy is interpreted as an operator on /^(R) with zero boundary conditions outside 
R). The spectral problem 

nvY = >^Y (A G M; V'G (2.1) 

has \Vh\ solutions ^ ^ ... ^ i.e., the ordered eigenvalues of the 

operator Hv- 

In this section, we provide asymptotic expansion formulas for the first K largest 
eigenvalues X^ y under conditions on the first terms of the variational series ^ 
^2 y ^ ^ y of the sample ^y and their coordinates Zi^ y G R defined by 

4,v = ^vi^k.v) (i ^ A: < |R|); here R = {zk,v- 1 < fc ^ |R|}. The results of 
this section follow simply from more general results of (Astrauskas 2008; 2012) and 
Section 2.4 of (Gartner and Molchanov 1998), where one finds more discussions on 
the relationship between ^y-extremes and the top eigenvalues of Hv- 

In this section, the proof of the statements relies on deterministic spectral argu¬ 
ments. It is worth mentioning that the (probabilistic) Feynman-Kac representations 
of the Green function and the principal eigenfunction of Schrodinger operators "Hy as 
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well as the related path decomposition techniques present powerful probabilistic tools 
for deriving the explicit upper bounds for the principal eigenvalue and eigenfunction 
of T-Lv (Gartner and Molchanov 1998; Gartner et al. 2007). However, this method is 
not explained in the present section; see Section 7 below for some aspects of these 
techniques related to the parabolic Anderson models. 

2.1. Preliminaries: rough bounds 

We start with the following simple bounds for eigenvalues provided \Vb\ ^ 2. 
Theorem 2.1. (i) For any V and any , 

< Ci,v and |A, K (2 < ^ < |Vh|). (2.2) 

(ii) For any V, any K ^ \Vh\ and such that mini^fe<;^x — 2 , ^ 2, we 

have that 

Cl V ^ \ V ^ Cl V + 2i^«: for all 1 I ^ K. (2.3) 


Proof. We repeatedly use the fact that the Kth eigenvalue Xk,v = Xk,v{^v) of the 
operator nAy + Cv is a nondecreasing function in each variable fv{x) tending to in¬ 
finity (AT ^ 1, a; G V); i.e., the monotonicity property of eigenvalues (Lankaster 1969, 
Theorem 3.6.3). 

(i) To estimate A^ we abbreviate •= Ci,v if a: = and C{x) := —oo, 
otherwise. Note that ^'(-) < Cv(’) ^ Ci,v i’^ Therefore, A^ is bounded from 
below by Ci.v? i-®-; tho principal eigenvalue of the operator kAv + Cy- Moreover, 
Aj is bounded from above by the principal eigenvalue of the operator kAv + Ci,v 
on P{V), which in turn does not exceed 2vk + Ci v’ since the norm of the Laplacian 
«:Ay is less than 2uk. Similarly, since each eigenvalue A, is bounded from above 
(resp., from below) by the Ah eigenvalue of the diagonal operator + 2vk on A(H) 
(resp., ^y — 2vk), we obtain (2.2) for I ^ 2. 

(ii) We need to show the lower bound in (2.3). Without loss of generality, we 
assume that > 0 (this may be achieved by shift transform of and A in the 
spectral problem (2.1)). Write Sy := {z^y,...,Zj^y}. We introduce the following 
functions: ({x) := fy{x) if a; G Sy , and is zero, otherwise; and further on, ((x) := 0 
if a; G £y, and ((x) := —oo, otherwise. Then ^y(-) ^ ((■) +C(’) in therefore, 
each eigenvalue A, is bounded from below by the corresponding eigenvalue A; y of 
the operator /iAy + fy + (y; here 1 ^ I ^ K. To estimate A; y, we rewrite the 
corresponding spectral problem in the form: 

(A - ^Ay - (y)^ = (y^ {X > 0, G f {V)) (2.4) 

and apply the resolvent operator t/y(A; fy) := (A ——^y)“^ to both sides of (2.4). 
Since GviX; Cv)Sz = X~^Sz for z G Sy , equation (2.4) is transferred to 

tp= Cy(^)V’(^)A"^(5z (A > 0); 
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here 6y{-) is the Kronecker delta function. Clearly, for each 1 ^ I ^ K, the pair 
A; y = and 'ip {-; A; y) = Sz^ ^^ (•) solves this equation. Summarizing, we have that 
^i,v ^ A; y = (1 < / ^ K), as claimed. Theorem 2.1 is proved. □ 

In Sections 2.2-2.5 below, we consider three classes of functions ^y: 

(J) Sparse distinct ^y-peaks dominate in the landscape of ^y as C t i-e-, 
^y possess properties like (1.9)-(1.11). Then the Kth largest eigenvalue A^ ^ is 
associated with an isolated peak ^r(K)yj so that XKy ^ ■ZT(if),v for some t(K) = 
ty{K) ^ 1 (Section 2.3). In particular, if the functions ^y possess extremely sharp 
peaks like (I.2)-(1.4), then the eigenvalue y is associated with the Kth largest 
value of ^y, viz., XKy ^ zxy (Section 2.2). In both cases, the lower bounds in (2.3) 
are achieved as C f 

(JJ) The landscape of ^y is dominated by flat islands of large values with an 
unboundedly increasing diameter. Then the largest eigenvalues are associated with 
such relevant islands. In this case, the upper bounds in (2.2) are achieved as V '[ 'E'' 
(Section 2.4). 

(JJJ) Similarly as in (JJ), bounded islands of large values prevail in the landscape 
of ^y. Then the asymptotic expansion terms of the principal eigenvalue A^ ^ fill the 
gap between its lower and upper bounds in (2.2) (Section 2.5). 

In the case of (J) we obtain the explicit expansion formulas for eigenvalues in terms 
of ^y-values. Meanwhile, for (JJ) and (JJJ) we restrict ourselves to a derivation of 
the second order expansion formulas for eigenvalues. 

2.2. Potentials with extremely sharp single peaks 

For ^ 2, let us write 


■ {'^1, V ’ '^2,V ’ * ■ * ’ V } ^ ^ 


(2.5) 


for the subset of coordinates of the first N largest values of ^y, and 


r„ „ = mm \z, „ - = min \x - y\ 


( 2 . 6 ) 


for the minimum distance between sites in Sy . For natural \ ^ K = Ky < N = 
Ny < |y| and p ^ 0, we introduce the following conditions on functions ^y: 


lim mm £p+i)AK,y(?i.v -?i+i,v) = oo where lim^^ ^ = oo. 


V l^l^K 

C := limsup < I, 
V iK.V 


(2.7) 

( 2 . 8 ) 


T N,V 

y logA^ 


= oo, 


M := limsup max max |‘Cy(2;)| < oo 

y lx—z;_v|=l 


(2.9) 

( 2 . 10 ) 
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and, finally, 


lim min 

y K+l^l^N 


^K,V 



2vk^ 

+ 7— 

S/C,V 




|a:-Zi,v| = l 


^K.v-iv{x) 


( 2 . 11 ) 


We write := oo, and sy(Z) := A - C,+i,v) ior I 

Theorem 2.2. (i) Under (2.7) with p = 0, we have that 


limsup max |A, |s^,(?) ^ const i(k, iz). 

Y l^l^K ’ ’ 


and 


limsup max max 


log|V'(a;; A, ^,) 


^ - 1 . 


y x^zi,V lx — Z, y I l0gS^,(/) 

(ii) Under (2.7)-(2.9) with p = 1, we have that 

const 2 ('t, v) 


limsup max A, ^ 


l-C 


and 


limsup max max 


log 'ip{x;X,Y 


y l^l^K x^zi^v |x — Z, I log 
(iii) If Y satisfies (2.7)-(2.11) withp = 2, then 


^ - 1 . 


( 2 . 12 ) 


limsup max 

y l^l^K 


^l,V ^l,V 


2vk^ 




ilv < 


M- const 3 ( k , v) 

(1 - Cf 


+ const i{K, v). 


and (2.12) holds true. 


Proof. We first note that condition (2.10) implies ^^^ y > —oo for any V D Vq. On 
the other hand, if Cat v = r^.^ y > 1, then A, v for all 1 ^ ^ ^ A^. The 

latter is shown by the same arguments as in the proof of the lower bound in (2.3). 

Now, assuming ^^ y > —oo and letting V f , the assertions of Theorem 2.2(i), 
(ii) and (iii) are derived from Theorem A.l(i), (ii) and (iii), respectively, in (As- 
trauskas 2012, Appendix A) with the abbreviations II := Sy , L := Y v s-od r := r^., y. 
□ 


In the case of conditions (2.7)-(2.9) with p = 2 (part (iii) of the theorem), we 
have imposed additional restrictions (2.10) and (2.11) to control the influence of the 
lower ^y-values on the correspondence Xk,v zk.v- If — ^ic+i,v) = 

0(1), then the lower ^y-values may essentially contribute to the expansion of the 
eigenvalues Xk,v and, therefore, the correspondence Xk,v t-t zk,v fails. See Section 
6.1 of the present paper where we consider the case of i.i.d. samples ^(-) in 1^ t 
with “weakly ” pronounced asymptotic peaks. 

Let ^(-) be an i.i.d. random field with the distribution function 1 — e”*^. We will 
show that, if Q satisfies the condition Q{t) = o(t^'''^) for p = 0,1 and 2 (“heavy tails” 
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and additional RV conditions as t —>■ oo, then with high probability satisfies 
the assumptions of Theorem 2.2(i), (ii) and (iii), respectively, where RT G N is fixed 
and N = [|R|^] for some 0 < 0 < 1/2; see Theorems 4.3(i), 4.5, 3.1 {R = 0), 5.3 and 5.4 
with 0 < e < 0. Therefore, Poisson limit theorems for the largest eigenvalues ^ are 
reduced to those for extreme values of i.i.d. random fields ^(-) or ^(- ) + 2iyK ^/(^(-) Vl) 
(Theorem 6.9). 

2.3. Potentials with dominating single peaks: the general case 

To simplify the proceedings, we need some notation and remarks. For N ^ 2 and 
Sy as in (2.5), we introduce the following function: ^v{x) := 0 if a: G Sy , and 
:= ^v{x) if a: G F \ Sy . Then 

^ S,v{z)5z + fy, 

i.e., ^y is a superposition of ^y-peaks and the noise component /'y. To exclude the 
trivialities, we assume that ■Cw v > each V (for the case when Cw v = 

and y > 1, see the proof of Theorem 2.2 above). For each z gV, let Ay(z) be the 
principal eigenvalue of the “single peak ” Hamiltonian nAy + ^v{z)Sz + ~ ^z) 

on l‘^{V). We associate the sites v G ^ with the variational series 

Ai y . Y ■ '^('^'r(|V|),v) 


based on the sample Ay; here V = y : 1 < ? ^ |F|}. 

Theorem 2.3. Assume that there are natural numbers 1 ^ K = Ky < N = Ny < 
\V\ such that the functions fy satisfy condition (2.9) and the following conditions: 



liin {^Ky - fN,v) = oo 

(2.14) 

and 

■ f • log (Ai,y - Ai+I,y) 

limini min - j -^ ^ (J. 

v rvr.y log (5(;+i)AA:,y - Cwt) 

(2.15) 

Then 

r log Ai^y - Ai^y 

limsup max - ^—j -^^ ^ — 2. 

y nyy log [^ly — ^N,v) 

(2.16) 

and 

log V’(a:;A,_^) 

limsup max max -4-—-G—-r ^ — 1- 

y I^KK x^z^(i),y |a; - 1 log(C,, V - V) 


Proof. 

We write 



y :={z G £y : Xv(z)^fj^ y+2i'K+h} where h:= _ 

(2.17) 
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By the first bound in Theorem 2.1(i), Xv{z) ^ ^v{z) for all -2 G £{/ ■ This combined 
with (2.14) gives 


Kv\>K (2.18) 

for any V D Vq- Finally, according to (Astrauskas 2008, Section 2.2 and Appendix B.2), 

^v{z) ^ Xviz) ^ ^viz) + 2 i'k'^/ h for any (2-19) 

In view of (2.18) and (2.19), the assertion of Theorem 2.3 is derived similarly as in 
the proof of Theorem B.3 in (Astrauskas 2008, Appendix B) with the abbreviations 
L := n := (2.5), B := (2.17) and r := (2.6). □ 


Note that conditions (2.14) and (2.15) of Theorem 2.3 are substantially weaker 
than (2.8) and (2.7), respectively, in Theorem 2.2. According to (2.16) and (2.19) 
with h —>■ oo as in (2.17), we obtain that Xiy = + o(l) uniformly in 1 ^ I ^ AT, 

so that the eigenvalues Xiy achieve their lower bounds in (2.3) as B ■ 

On the other hand, from (Astrauskas 2008, Appendices A and B) we know that, 
for each z & the eigenvalue Ay(z) is the maximal solution to the equation 


Qv{X]^v]z,z) 


1 


( 2 . 20 ) 


here Gv{X;^v',‘,') is the Green function of the Hamiltonian nAy + Cv on P{V), so 
that Gv{X] Sy\ z, z) is expanded over paths: 


Qv{X\ ^v\z,z) = ^ Acfol n 
r v^v 


( 2 . 21 ) 


where the sum is taken over all paths F : vq '■= z ^ vi ^ ^ Vm '■= z in V such 

that \vi — Vi-i \ = 1 for each 1 ^ z ^ m and each m G N, ni,(r) denotes the number 
of times the path F visits the site v € V, |r| := y2vev ^i'(r) — 1^0. Substituting 
(2.21) to the left-hand side of (2.20) and iterating this with respect to the eigenvalue 
A = Ay(z), we obtain the explicit expansion formulas for Ay(z) {z G E^.v) presented 
as a power series in the variables ^v{z) and ^y(x) (|x — z\ ^1), in particular, 


Xv{z) = ^v{z) + 


1 


+0 E 


\r-,yi^viz)-^vix) 
1 


-.1=1 (Cv(^)-Cy(a;))(^y(2;)-^y(j/))(^y(z) - ^v{u)). 


\x-z\ 
\y-z\ = i 


( 2 . 22 ) 


as y t 

Remark 2.4. With notation at the beginning of Section 2.2, assume that there 
are natural numbers 1 < N = Ny < \V\ such that the functions ^y satisfy the 
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conditions: limy = oo and limy(^j ~ ^n,v) — o®- Then there are constants 
consti = consti(K, > 0 such that 

1 , ^1 consti consto 

IA 1 .y-a.yK 7 -7^ +-^ = o(l) (2-23) 

a.y ~ sAf.y ^Af.y 

as y t Z""- 

As shown above, limit (2.23) follows from the assumptions of Theorem 2.3 or The¬ 
orem 2.2(ii) with AT = 1. To prove (2.23) under the (weaker) conditions of Re¬ 
mark 2.4, we first observe from Theorem 2.1 (i) that ^ ^ a v ^■ Second, 

we apply Lemma 2.8 below with R := y —>■ 00 to see that the eigenvalue Aj 
is bounded from above (with accuracy 0{R~^)) by the maximum of local principal 
eigenvalues of the single-peak Hamiltonians over all balls in V of radius R. Using for¬ 
mulas (2.20)-(2.22) for such local principal eigenvalues, we finally obtain the bound 
Ai.y a a.y + consti/(a.y — K.y) + const /R, as claimed. □ 

Remark 2.5. Theorems 2.2, 2.3 and Remark 2.4 include the condition on asymp¬ 
totic sparseness of ^y-peaks: for instance, r^^ y —>■ 00 as U f Molchanov and 
Vainberg (1998; 2000) studied the existence and location of spectral components (pure 
point, absolutely continuous, etc) of the Schrbdinger operators kA -f aO ici P{Z'') 
with sparse deterministic potential ^i') = (• )> where amplitudes Ok are 

bounded and {zk} is a rare subset of Z'^, for example, r„ := min,^^ \zi — z„| —>■ 00 as 
n —>■ 00 . □ 

Let ) be an i.i.d. random field with the distribution function 1 — e”*^. We will 
show that, if the tails e~'^ are heavier than the double exponential function (i.e., 
log(5(t) = o{t)) and satisfy additional regularity and continuity conditions at infinity, 
then with probability one ^y satisfies the assumptions of Theorem 2.3, where AT G N 
is fixed and N = [|U|®] for some 0 < 0 < 1/2; see Theorems 3.1 {R = 0), 4.6 {p = 00 ) 
and 4.8. Therefore, Poisson limit theorems for the largest eigenvalues A^ y are reduced 
to those for extremes of nonlinear functions (2.22) on ^y (Theorem 6.2). 

2.4. Potentials with dominating flat increasing islands of high values 

Let B/j(z) := {y GV: |j/ — zK A} denote the ball in V with center z GV and radius 
A ^ 0. The following theorem gives a simple condition on ^y which ensures that the 
largest eigenvalues A^ achieve their upper bounds in (2.2) as U f 

Theorem 2.6. If 

lim limsup min (A ,, — min fy(x)) = 0, (2.24) 

then, for arbitrarily fixed K gN, 

-K,v) = 2yK- 


Proof. Because of Theorem 2.1(i), we only need to show the lower limit bound 


lirmnf (A^ ^ 2i/k. 


(2.25) 
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From (2.24) we see that there exist a sequence 0 < eij —>■ 0 and sites zy & V such 
that Cy(’) ^ Ci,v ~ for any R ^ Rq and any V D Vo{R). Abbreviate 

:= if a; G MrIzv), and := —oo, otherwise. Since Cy(-) > 

in V and ^ ^i,v ~ ^R fo ^Ri^v) for R and V as above, the monotonicity 

property of eigenvalues implies that Xk,v ^ +^iA ~ ^R^ where Kill 

eigenvalue of the operator kA on with zero boundary conditions. Since 

tends to 2vk letting first V 'I' Z" and then R ^ oo (Kirsch 2008, Section 3.1), 
this estimate implies (2.25), as claimed. □ 

Clearly condition (2.24) is fulfilled if and only if there are a sequence Ry —>■ oo and 
sites zy GV such that 


lim (Ci,v - , min Cv(a;)) = 0. 
V |a:— 


From the proof of the theorem we know that the eigenvalues y of the operator 
Ry = n/Xy + in C are approximated by the corresponding local eigenvalues in the 
regions B/j^,;opt := ^Ryizy) C V where ^y(-) is close to i.e., relevant regions. 

Let ^(-) be an i.i.d. random field with the distribution function 1 — e”*^. We 
will show that, if the tails e~^ are lighter than the double exponential function 
(i.e., log Q(t) —>■ oo) and satisfy additional RV conditions at infinity, then with 
probability one ^y satisfies the assumption of Theorem 2.6; see Theorem 4.6 with 
p = 0 and Theorem 3.1(i) for any large R and 0(-) = 6 = const. In this case, we obtain 
the second order expansion formulas for the largest eigenvalues A^^ (Theorem 6.14). 


2.5. Potentials with dominating bounded islands of high values 

In this section, we describe a class of deterministic functions (potential) ^y : V ^ 
[—oo;oo) for which the asymptotic terms for the principal eigenvalue Aj (V f Z"^) 
fill the gap between its lower and upper bounds in (2.2). We use the variational 
arguments developed by Gartner and Molchanov (1998). To formulate the results, 
we need some abbreviations and remarks related to the variational problems. To 
emphasize the dependence of Ai^y on the sample ^y, we denote by A(^y) := Ai^y the 
principal eigenvalue of the operator Ry = nAy + ^y on l‘^{V). As in Section 2.4, let 
Bij(z) C C be the closed ball of radius i? ^ 0 centered a.t z G V, and let B^j := Bii;(0). 

Given a ball B C G, let (x) := ^y{x) if a; G B, and (x) := —oo, otherwise. As 
before, Ry := nAy + is interpreted as an operator with zero boundary conditions 
outside B. We write 

iv{x) = +hy{x) {xGV), 

where the function hy ^ 0 admits the interpretation as the shape of ^y-values close 
to the maximum . Note that 

X{^y) = + X{hy) and A(^®) = + A(/i® ). (2.26) 

For a fixed constant 0 < p < oo, we are interested in the following supremum of A(/i®) 
over : B —>• [—oo; 0]: 

sup{a(/i®): ^e'‘fo)/^<l}. 
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This variational problem is equivalent to the corresponding variational problem in 
terms of the functionals 


5'”(p):=y^ \/p{x)/S.yfp{x) where p{y) = 0 for y G 


and 

/“(p) := _^p(a;) logp(x) 

for p(-) G P(B), the set of probability measures on B. More precisely, for a sequence of 
balls B/j C the following formulas hold true according to Rayleigh-Ritz theorem 
and (Gartner and Molchanov 1998, Lemmas 2.17 and 1.10): 


sup{A(h®«): ^ < l} 

= sup{A(h®«): = 

= sup ( (p) — p/®” ( 


pe'P(BR) 

where the right-hand side of (2.27) converges (as R —>■ oo) to 
sup (^kS{p) — pl{p)) ='■ 2vKq{p/ k). 

peV{Z'') 


(2.27) 


(2.28) 


Here S{p) and I{p) are the corresponding functionals on 'P{Z‘')= the set of proba¬ 
bility measures on lattice Z*^. It is easy to check that q : 1R+ —>■ (0; 1) is convex, 
strictly decreasing and surjective function; q{0) = 1 and < 7 ( 00 ) = limp_>oo q{p) = 0 - 
Moreover, q{p) = (2plogp)“^(l -I- o(l)) as p —>■ 00 (Astrauskas 2008, Proposition 2.1 
and Corollary 4.5). The supremum on both sides of (2.27) and (2.28) is attained 
(Gartner and Molchanov 1998, Sections 1 and 2.4). Denote by h®pt the maximizer for 
the variational problem on the left-hand side of (2.27). Then p®pt is the maximizer 
for the right-hand side of (2.27) if and only if ft,®p^ = plogp®pt. For this and further 
properties of the maximizes in (2.27) and (2.28) in the limit case R = 00 , see (Gartner 
and den Hollander 1999, Sections 0.3 and 0.4) and (Gartner et al. 2007, Sections 1.3 
and 3). 

The following theorem tells us that, under reasonable conditions on , the principal 
eigenvalue A^ ^ of the operator Hy = kAv -f Ci.v -f hy in H is approximated (letting 
first V '[ and then i? —>■ 00 ) by the local principal eigenvalue of the operator 
restricted to the regions B/j-opt := B/{( 2 y) C V where hy is close to h®pj, i.e., relevant 
regions with optimal potential shape. 

Theorem 2.7. Given a constant 0 < p < 00 and a sequence R —>■ 00 , assume that 
functions ^y satisfy the following conditions: 


lim limsup max 

R—^oo Y zGV 


E 

y^MRiz) 


exp 


( £.v{y) - Ci,v 

p 


} 


^ 1 


(2.29) 
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and 


liminf liminf max min {^viv) - ^i,v “ h^ptiv “ z)) > 0. (2.30) 

R^oo V zev yeBR(z) . / 


Then 


lim {K,v - Ci.v) = 2RKq{p/K). 


(2.31) 


Proof. Limit (2.31) follows from the results of (Gartner and Molchanov 1998, 
Section 2.4) under the stronger conditions on including sparseness of clusters of 
^y-extremes. To prove (2.31) under conditions (2.29) and (2.30), we apply the same 
arguments as in (Gartner and Molchanov 1998, the proof of Theorem 2.16) combined 
with the following lemma by Biskup and Kdnig (2001), which is slightly modified for 
the operator TLy with periodic boundary conditions: 

Lemma 2.8. (Biskup and Konig 2001, Lemma 4.6). For each R G N, V Z) Vo{R) 
and each fv, 

H^v) ^ Ci.v + + const 

for some (universal) const > 0. 

We first obtain the upper bound for Condition (2.29) implies that there is 

a sequence 0 < Er ^ 0 such that 

max ^ ex.p [hviy)/p} < ex.p {er/ p} for each VdVo{R). 

yeMniz) 

In view of (2.26), this estimate and Lemma 2.8 yield that 

M^v) - Ci,v 

^ sup {A(/r®«): h{-) < 0, ^ 

^sup{a(/i®«): h{-)Ef0, ^ <i}+er + ^^ 

ye«R 

for V as above. Taking the limit as first V f and then i? —>■ oo, and using 
(2.27)-(2.28), we arrive at 

limsup (Aj < 2vKq{p/K). (2.32) 

By combining condition (2.30), the monotonicity property of eigenvalues and asser¬ 
tions (2.27)-(2.28), similarly as in the proof of (2.25) we obtain the lower bound 

A(?y) -?i,v > Hhopt)+oil) -t 2vKq{p/K) 

letting first V f IT and then R —>■ oo. This and (2.32) conclude the proof of Theo¬ 
rem 2.7. □ 
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If p = py —>■ oo in Theorem 2.7, then the “relevant” regions ]B_R;opt shrink to single 
sites and, therefore, we are in the situation of Theorem 2.3. Meanwhile, if p = py —>■ 0, 
then we stick to the result of Theorem 2.6. 

Let ^(•) be an i.i.d. random field with the distribution function 1 — e~'^. It will 
be shown that, if the tails e”'^ are the double exponential (i.e., \ogQ{t) 1/p) 
and satisfy additional RV conditions at oo, then with probability one ^y satisfies the 
assumptions of Theorem 2.7; see Theorem 4.7 and Theorem 3.1(i) for arbitrarily large 
R and eiiiy) ~ 1 “ exp{/i“«t(p)/p} (y G IB/^). ^For continuous t/, see Corollaries 2.7, 
2.12 and 2.15 in (Gartner and Molchanov 1998).) In this case, the second order 
expansion formula for Aj ^ holds true (Theorem 6.19). 


3. CLUSTERING OF HIGH-LEVEL EXCEEDANCES OF I.I.D. RANDOM 
FIELDS 

Let ^(x), X G be an i.i.d. random field with the cumulative hazard function Q. The 
main task of the present section is to investigate the almost sure asymptotic structure 
of clusters (“islands” ) of bounded size formed by exceedances of the sample ^y as 
V" t With the abbreviations in Section 1, we also need additional notation. For 
9 < 1, put := /((I — 9) log |R|) where / := Q*~. (Without loss of generality, we 
write := sup{t: Q{t—) = 0}, so that almost surely ^{x) > Ly^ for each x.) Let 
]B/j(z) := {x G Z'': jx — 2 ;| ^ R}, and 1 B_r := Bi^(O). For fixed i? G N U {0} and a 
function 0 _r(-) : ]B_r —>■ (—oo; 1], we denote by Yr the set of balls Bi^(z) C V, and 

£^g := {Bi^(2;) C V: ^(y) > Lvfi^(y-z) for all y G B/j(z)}, 

the subset of clusters of ^y-exceedances in Yr over the level function Ly g^^i^.y We 
abbreviate 

r(£:/(e) := min{dist(B,B'): Bg^:/?^, V>' B^B'} if \£^g\-^2, 

and r{£yg) := if {Sygl ^ 1, by convention; here dist(B,B') stands for the 

lattice l^-distance between balls B,B' C R. If i? = 0 and 9 := 0(0), then Sy^g := Syg 
shrinks to the subset of single ^y-exceedances, so that 

r{£v.e) = ^^H\x-y\- xGSyg, yGSyg, x^y}. 

To formulate the main result of this section, we also need the following abbreviations 

y-R ■= ^niy)) > 0 and 0max,fl := niax {^^(p): 0i^(J/) < l}. 


Theorem 3.1. (cf. Theorems 2.2-2.7). 
following almost sure limits hold true. 

(i) If yLR < 1, then 


liminf 

V 


fog \£% 
log|fo| 


For arbitrarily fixed i? G N U {0}, the 


^ 1 — y-R- 
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(ii) If fiR < 1 and, in addition, Q satisfies the condition 


lim 9^ = 1, 

ntQ Q(t) 

then 

(iii) If Hr > 1 and Q satisfies (3.1), then 

lim \ £vf)\ = 0 . 

(iv) If 0niax,R < Hr < 1 Q satisfies (3.1), then 

\ogr{£^e) 2 hr-1 

lim-^—— =- 

V iog|y| u 


(3.1) 


Remark 3.2. (a) Clearly, for arbitrary Q and hr < !> 

(= the mean number of clusters of exceedances) 

= |Vk| n > ^vMy)) 

y&R 

^ const \V\ exp { — ^ Q{d^v,9Riy)—)} ^ const —>■ oo, 

y&R 

according to Lemma A. 11 (iii) in Appendix. 

(b) On the other hand, by Lemma A.ll(iii), condition (3.1) implies the asymptotic 
formula 

Q(/(®)) = s + o('5) as s —>■ oo, (3.2) 

which in turn yields, for hr ^ Oj the upper bound 

logElf^gl if {I- HR + o(l)) log \V\ 

as \V\ —>■ oo. 

Remark 3.3. (see part (iv)). If 9jna,^,R < /r/j < 1, then 1/2 < hr < I- 

Remark 3.4. In the case R = 0 and 0 < 6* < 1/2, i.e., single rare ^\/-peaks. 
Theorem 3.1 was proved by Astrauskas (2001). For the Gaussian random field f{-) 
with correlated values, the case R = 0 was studied by Astrauskas (2003). Here 
the results depend slightly on the correlation function of ^(•). Finally, assertion (i) 
generalizes Corollary 2.15(b) in (Gartner and Molchanov 1998) where the continuity 
of Q is assumed. 

Proof of Theorem 3.1. To simplify the proof, we assume throughout that 0r{-) = 
0{ ) < 1 in Mr. The general case is treated similarly. 
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(i) We denote by Yr C Yr the maximal subset of nonintersecting balls ]B/j( 2 ) 
in V, so that |V/{| x \V\. The claimed bound is proved by estimating ISyg H Vfl| 
similarly as in the proof of Theorem 1 in Astrauskas (2001), where the exceedances 

^ Lv,e} (x € V) are replaced by mutually independent (multiple) exceedances 
{^(•) ^ in ]Bfl(z)} (]Bfl(z) G Yr). In particular, if Q{tQ—) = oo, we obtain 

that, for any —1 < (5 < 0 and V f 

< (i + <5)E|£:«,nVfl|) 

^ exp { — const (6) ■ '^fl| (1 + o(l))} 

for some const ((5) > 0. Since the right-hand side is summable over V according to the 
assertion of Remark 3.2(a), we conclude the proof of (i) by using the Borel-Cantelli 
lemma. 

(ii) We only need to estimate from above. Fix a function 0'(-) : B/j —>• (—oo; 1) 

such that 0'(-) > d(-) in B/j, and pick a constant <5 > 1 — where fx' := ~ 

0'{y)). We then apply Chebyshev’s inequality and the assertion of Remark 3.2(b) to 
find that, for any V ZX Vo, 

^{\£%,\ > |R |0 ^ < |y|-const (3 3 ) 

where const = const ( 6 l'(-), 5) > 0. Choose a subsequence {V{1)-. / G N} C {V} such 
that 

V{1) monotonously increases and |R(Z)| = 2^(1-I-o(l)) as I ^ oo. (3.4) 

Since the right-hand side of (3.3) is summable over the subsequence {V{1)}, the Borel- 
Cantelli lemma implies that almost surely g/1 ^ |C {l)\^ for all I ^ ^o(w). Because 
of the monotonicity of Syg in we obtain that with probability 1, for any V 

such that V{1 — 1) CV CV{1) and any I ^ lo{uj;9{-),9'{-)), the set £yg is contained 
in i^^yp) gi, therefore, 

< \^vii), 0 ' \ < |C(Z)|‘^ < const|y|^ 

Since 9'{-) > 9(-) and S > 1 — /x' are chosen arbitrarily, this estimate yields the upper 
limit bound for log |i?yti|, as claimed. 

As in part (ii), it suffices to prove the assertions of (iii)-(iv) for the subsequence 
{y(/)} (3.4) instead of {V}. 

(iii) We note that £yg 0 if and only if there exists Rr{z) C V such that ^(•) ^ 
Lyfii^._z) in Mr{z). According to the assertion of Remark 3.2(b), the probability of 
the last event does not exceed E|£lyg| < \V\~p for some 0 < p < —1 + /r. Therefore, 

^ {^v,e ^ 1^1”^- Since the latter is summable over {C(0} (3.4), the Borel- 

Cantelli lemma yields that almost surely £y(^i^ g = 0 for all I ^ lo{<^), as claimed. 

(iv) With Yr C Yr defined in part (i), the almost sure upper bound for r(^£yg fl 
Yr) ^ x{£yg') is derived similarly as in the proof of Theorem 2 of Astrauskas (2001) 
where the exceedances {C(a^) ^ Lvfi} (x G V) are replaced by mutually independent 
(multiple) exceedances {^(•) ^ Ly^g(^._z) in Bh(z)} (Bfl(z) G Yr). 
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To obtain the lower bound for r{£yg), we first note that the event {r{^£yg) = 0, 
I'^vel ^ 2 } implies that there exists C V such that ^(•) ^ Lvfi{.-z) in 8 ^( 2 ) 

and ^{y) ^ for some y £ (B3^(z) \B^(z)) fl V. Therefore, as in the proof 

of (iii) we obtain that, for fixed y \ B/j and for any D Vb, 

= 0 ^ 

^ const |y|P (^(O ^ Lvfi(.) in Br, ^(y) > < 1^1"'’ 

for some 0 < p < y — 0max- Second, similarly as in the proof of part (i), we find that 
P < 2^ ^ |p^|-const V Z) Vq and some const > 0. Summarizing these 

bounds and picking 0 < £ < {2p — \')lv arbitrarily, we get that, for any F D Vb, 

iiF{l ^r{S^g) <\V\^,\£^A + (3.5) 

for some const ^ > 0 , where the last probability is estimated similarly as in the proof 
of Theorem 2 of (Astrauskas 2001) with mutually independent (multiple) exceedances 
{^(•) > T\/,e(.-z) in Bii;(z)} instead of {^(a;) > Lv,e}- Since the right-hand side of (3.5) 
is again summable over {1^(0} (3.4), we conclude from the Borel-Cantelli lemma that 
almost surely r(^£y^i^ g) > 1^(01“^ fo'' nny I > lo{ui',e), as claimed. This completes the 
proof of Theorem 3.1. □ 

Remark 3.5. By the same arguments as in the proof above, the assertions of 
Theorem 3.1 are extended to the following class of high-level exceedances: 

For fixed i?GNU{0},we denote by Sr the set of all subsets U C Z'' , the diameter of 
which does not exceed R. Let Yr := {U £ Sr: U C V}. For a fixed set of functions 
©fl := {^c/.fl(’) G (— 00 ; 1)1^1; U £ §fl}, let SyQ C Yr be the subset of elements 
U G Yr such that ^(•) ^ in U. l.e., SyQ consists of clusters of exceedances 

in Yr over level functions Denote by r {Syo) the minimum distance among 

elements U,U' G £y.e^ ^ ^ ■ Finally, let pR := '^yizui^ ~ Su,Riy)) be a positive 

constant independent of [/ G Sr, and write 9max,R '■= sup^/gg^ uiOKy^ij 0u,Riy)- With 
these notation for SyQ and r (Syo), the almost sure assertions (i)-(iv) of Theorem 3.1 
hold true. 


4. SPACINGS OF ORDER STATISTICS OF I.I.D. RANDOM FIELDS 
4.1. Spacings of consecutive order statistics 

We first formulate the results for the exponential order statistics 77 ^ y and their spac¬ 
ings, which are then transferred to = /(r/^ y) under appropriate conditions for 

/• 

Note that the random variables 

Vl,V - 1)(?7|V1-1.V -»7|vi,v)>l^l^ivi,v 

are mutually independent exponentially distributed with mean 1; see, e.g., (Shorack 
and Wellner 1986, pp. 336). This property immediately implies the first assertion of 
the following lemma. 
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Lemma 4.1. (i) For fixed K 


IP(^l, V ^2, V ^ ? ■ * * 5 Vk-1,V Vk,V ^ ^K-1 ? Vk,V 1^1 

= ( n ^K-l)l I - e-^} ds 

for all ti'^ Q (1 ^ I ^ K — 1) and all t G R. 

(ii) For an arbitrary sequence {K^} such that 1 ^ Ky ^ |V^|, 


limsup WKy max 




< oo in probability. 


Proof. Let us show (ii). Write ?7n,|+i v d. By Kolmogorov’s inequality in 
(Shorack and Wellner 1986, pp. 843), we have that 


P max 




l—k 




> 


c 

1C 


1 / 2 ' 


K / 1 \" 2 

for any C > Cq and any V D Vb(C'). Combining this bound with the following simple 
estimate 


1 l^h 1 


K,, 


3 Co), 


we obtain the claimed assertion of (ii). 


□ 


The almost sure asymptotic behavior of the random variables ij^^ y and rjj^ y —rij^y.^ y 
(|K| —^ oo) is more intricate. 

Lemma 4.2. For any fixed constants K gN and m G N\{1}, the following almost 
sure limits hold true. 


(i) 

(ii) 

and 

(iii) 


. . . log(»7K,v -??K+i.v)+Efc2^1ogil^l 


liminf 

v 


limsup 

V 


= - 1 , 


log„^ iv^l 

Ik.v - Vk+i.v - ^ J_ 

K 


limsup max 
y i^/^|y| 


log™ 1^1 

1 1^1 
^i.v - log — 


1 


log log|K| 


= 1 ; 


here \ ■ := 0. 

Proof. Assertions (i) and (ii) follow from more general results for exponential spac- 
ings in (Astrauskas 2006, Corollary 12). Assertion (iii) follows from the corresponding 
strong limits for the uniform order statistics C|vi-fc+i,v (Shorack and Wellner 1986, 
pp. 408 and pp. 420-424) via transformation rj^ y = — logC|vi-fc+i.v (1 ^ ^ 
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□ 

We now turn to the case ^ = f{ri^ y). For p ^ 0, we denote by A\Y^ the class of 
functions / := Q'^ that satisfy 

lim /(s)^(/(s + c) —/(s)) = oo for any c > 0, (4.1) 

and by AHq the class of functions / that satisfy 

lim /(s)^(/(s + c) —/(s)) = 0 for any c > 0, (4.2) 

and, finally, O^IfP stands for the class of / satisfying 

/(s)^(/(s + c) -/(s)) X 1 as s —>■ oo, for any c > 0. (4.3) 

We see that, if / is in or O^IfP, then the right endpoint tq is infinity or, 

equivalently, /(s) —>■ oo as s —>■ oo. Of course, Ali^ includes the trivial case of finite 
tq > 0. The characterization of AIl^, AIV^ and 0^11^ is given in Lemmas A.6, A.7 
and A.8 of Appendix A respectively. In particular, the functions / := G OAII^’ 
are associated with Weibull type distributions 1 — e”*^, where Q{t) x as t —>■ oo. 

Theorem 4.3. (Cf. (2.7) and (Astrauskas 2012; 2013)). For fixed natural K > Ifi 1 
and real pfiQ, we have the following limits in probability. 

(i) If f G AW^, then 

- ^K,v) = oo. (4.4) 

(ii) If f G AAl, then Mmy Ck,v(^i.v - ^k,v) = 0- 

(hi) If f G OAW, then “ ?k,v) - 1 as 1^1 -t oo. 

Proof. Using notation (1.22)-(1.24), rewrite the left-hand side of (4.4) in the form 

KVk.vT ifiVK.v + {Vi,v - VK,v)) - fiVK.v)) ■ 

The claimed assertions follow by applying Lemma 4.1(i). □ 


To obtain these limits with probability 1, we need the stronger conditions for /. 
Let us abbreviate 

m— 1 

dm, 7 ( 5 ) := s( log, s) (log,„ s)^+^ (s ^ So). 


Theorem 4.4. (Cf. (2.7)). For fixed constants AT G N andp 0, the following almost 
sure limits hold true. 

(i) If lims_>oo f{sY{f{s + I/dm.-yY)) ~ /(s)) = 00 for some m G N and 7 > 0, then 

- ^K+i.v) = 00 . 


(ii)//lims_>oo/(s)^(/(s-l-AT ^ logdm. 7 (s)) —/(s)) = 0 for some m G N andj>0, 
then 

l™CK-i-i,y(CK,v ~ ?K+i,v) = 0- 
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(iii) //lims^oo (/(s + logs) - /(s)) = 0, then 

lim(Cx,v -/(log|V"|)) =0. 

Proof. Assertions (i)-(ii) follow by the same arguments as in the proof of Theo¬ 
rem 4.3, where one applies Lemma 4.2 instead of Lemma 4.1(i). Assertion (iii) follows 
from Lemma 4.2 (iii). □ 

4.2. Spacings of intermediate order statistics 

We denote by P /<2 the class of functions / := Q'^ satisfying the condition 

f ( (\ — 

limsup - -T--T - < 1 for some 0 < e < 1/2. (4.5) 

s—>oo /(s) 

Class (4.5) is characterized in Lemma A. 13 (Appendix A). 

Theorem 4.5. (Cf. (2.8) ). Assume that f G P /<2 (4.5). Then for fixed K gN 
and 9 > e, almost surely 

limsup < const < 1. (4.6) 


Proof. By Lemma 4.2(iii), with probability one the random variable = f{rij^y) 
is bounded from below by /(log \V\ — 2 log log |y|) and ^ is bounded from above 
by/((I —d) log |C|-1-2 loglog |C|) for each C D Vb(w). Substituting these bounds into 
the left-hand side of (4.6) and using (4.5), we obtain the claimed assertion. 

□ 

We denote by RVp the class of nondecreasing functions g : R+ —>■ 1R+ such that, for 
any c > 1 , lims_>oo 5 (cs)/( 7 (s) = d’. l.e., g is regularly varying at infinity with index 
0 ^ p ^ oo. The case p = oo (resp., p = 0) indicates a rapid variation (resp., slow 
variation) of the function g. See Lemma A.3 in Appendix A for a summary of the 
well-known properties of the class RVp. 

Theorem 4.6. (Cf. (2.14) and (2.24)). For some 0 ^ p < oo, assume that G RVp. 
Then, for all eonstants 0 ^ e < 6 < 1, almost surely 

^ ^ 

(^[|V|£].V “ ‘^[|V|9].v) ^ 1-9' 


Proof. To prove this assertion, use Lemma 4.2(iii) and Lemma A.3(ii) similarly as 
in the proof of Theorem 4.5. □ 

The following statement is closely related to the result of Theorem 4.6 with 0 < 
p < oo. As in Section 2.4, let Bi^( 2 ;) denote the closed ball in V with the center z £ V 
and the radius RfiO. 
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Theorem 4.7. (Cf. (2.29) and (Gartner and Molchanov 1998)). For some 0 < p < 
oo, assume that G RVp. Then, for any fixed R gN, almost surely 

limsup max ^ exp { {^{y) - Ci.v)/p} ^ 1 - 

^ veVR{z) 


Proof. For continuous Q, this assertion is a straightforward consequence of Corol¬ 
lary 2.12 in (Gartner and Molchanov 1998) and Theorem 4.4(iii) above. (In view 
of Lemma A.3(ii), the condition of Theorem 4.4(iii) follows from the assumption 
of Theorem 4.7). If the continnity condition on Q is dropped, one applies slightly 
modified arguments based on the technique of function inversion, e.g.. Lemma A. 11 
in Appendix A. □ 

Recall that the conditions of Theorems 4.6 and 4.7 are discussed in Lemma A.3. 
In particular, the assumption of Theorem 4.7 implies that log(5(t) = tf p o{t) as 
t ^ oo, i.e., the double exponential tails e”*^. 


4.3. Minimum of spacings 

We first recall some notation from Section 2.3. For fixed 0<d<l/2, we write 
£,{x) := f{x) if f{x) < /((I — 9) log |R|), and ^{x) := 0, otherwise. For any z gV, let 
\{z) denote the principal eigenvalue of the “single peak” Hamiltonian nAy +^{z)6^ + 
^v(l — SJ in l‘^{V). As in (2.13), let \k,v denote the ATth extreme order statistics of 
the random field Xy. For k = 0 and 0 < £ < d, we know from Theorem 3.1(i)(i? = 0) 
that with probability one Xk,y = fk,v for all 1 < fc ^ \V\‘^ and all large V. For 
any k 0, we are interested in the asymptotic behavior of the minimum of the gaps 
Xk,v - Xk+i,y (1 ^ k ^ \V\^) defined by 

:= min {Xk,v - Xk+i,y ■ 1 ^ k ^ \V\^}- 

Given a constant p > 0, we say that the function F : R —>■ R is log-Holder continu¬ 
ous of order p > 0 at infinity, if F satisfies the following condition: 

|F(t-t s)-F(t-s)| llogsl'^ = 0(1) 
as t —>■ oo and s 0 simultaneously. 


Theorem 4.8. (Cf. (2.15)). Let tq = oo, Q and 0 < e < 6 < 1/2, and assume 
that the distribution tails e~'^ are log-Holder continuous of order p > 0 at infinity. 
For K > 0, assume additionally that e-f G RVoo- Then almost surely 


limsup 

V 


log{-log(5'^,^ A 1)} 

log I HI 




1 -I- £ 
F 


Proof. The assertion follows from Lemmas 3.5 and 4.3 in 
Theorem 3.1(ii) above, where R = 0 and 0 < 0 < 1/2. 


(Astrauskas 2008) and 
□ 
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In (Astrauskas 2003), the results of Theorem 4.8 are extended to the Gaussian 
random helds with correlated values. 

We end this section with some generalization of Theorem 4.3(iii) for the functions 
/ G OAXYP (4.3) associated with Weibull type distributions 1—e“^, where Q{t) x 
as t —>■ oo. 


Theorem 4.9. (Cf. Lemma 4.2 in (Astrauskas 2013)). For some p ^ 0, assume 
that f G OAIIP (4.3). Then, for arbitrarily fixed constants A'gN, 0<e<l and 
any sequence {ny} C N such that Uy = 0(|y|^), we have the following limits in 
probability: 

W as 1^1 ^OO, 

and 

(ii) 


0 < liminf min 


4,v) 


1 


log I 




< oo. 


Proof. Assertion (i) follows from a combination of the formula ^ = fiVkv)’ 
Lemma 4.1(ii) and the limit /(s) x gi/lp+i) as s —>■ oo (the latter follows from 
Lemma A. 8 (iii) with a(-) = const). Assertion (ii) is shown by combining the formula 
fk V ~ fi'hk v) Lemmas 4.1(ii), A. 8 (iii) similarly as in the proof of Theorems 4.3 
and 4.5 above. □ 


5. NEIGHBORING EFFECTS FOR EXTREMES OF I.I.D. RANDOM 
FIELDS 

We finally study the asymptotic properties of -values neighboring to -peaks. It 
is then straightforward to extend the results for ry(-) to 5 (-) = /(?y(-)). 

The following lemma tells us that, for fixed y ^ 0 and for small e > 0, asymptotic 
properties of the random variables ry( 2 ;^ v + J/) (1 ^ A ^ 1 ^ 1 “^) s-nd their extremes are 
the same as in the case of exponential i.i.d. random variables. 

Lemma 5.1. For fixed y G Z‘^\{0}, 0 < £ < 1/2 and a sequence of integers 
K := Kv = 0(|y|®), the following assertions hold true. 

(i) limy P(?y(z^ + y) > f) = e~* for all t > 0 . 

(ii) If, in addition, K := Ky —>■ oo, then 

limP( ^ma>^?y( 2 :, + y) — log AT < t) = exp{— e“*} for all t G M. 


(iii) 


v{zi,v+y) . 

lim limsup max- 1 

M^oo y logi 


0 in probability. 
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Proof. Here and in the sequel, we need the following key statement (which is 
frequently used in (Astrauskas 2013) as well). 

Lemma 5.2. Fix a finite subset U C Z'^\{0}, U 0, and a sequence of nonrandom 
realfunetions {Dftjj): t^ G (IgN). Abbreviate ri{z;l):=Di{{r]{z+x): xGU}) 

for z & JF and I G N. Finally, pick a sequence of integers K := = 0(|y|^) for 

some 0 < e < 2 • Then, for any V and any t G R, 


max nlz ,. 

1<1<K 


l)^t] -Y[FipiO;l)^t) 
^ 1=1 


<3|y|- 


where const > 0 does not depend on V and t. 

Now, part (i) of Lemma 5.1 follows from Lemma 5.2 with U := {y}, where Dx(t,j) = 
ty and Di{t,j) = 0 for I ^ AT. Part (ii) follows from Lemma 5.2, where Dft^j) = ty 
{I G N), combined with Lemma 4.1(i). Finally, by Lemma 5.2 with Dft^j) = ty/ log/ 
(/ ^ 2), we derive that, for any small 5 > 0, 

limsup ?(" max > 1 + (jW V 0 

y log/ ) 


and 


limsup f( max | < 1 - jW TT (l - = 0 

y \M^K\v\'^ log/ ^ iZm 

as M —>■ oo, i.e., assertion (hi) of Lemma 5.1 is proved. □ 

Proof of Lemma 5.2. Fix a constant 9 G (e, 5 ), so that K := ^ f 1^1^ for oach 

V 0 Vq. Write := (1 — 6) log \V\. Denote hy dV the subset consisting of sites 
at which p{-) exceeds the level L^, and let r(£v) be the minimum distance among 
sites in Ey', cf. the notation at the beginning of Section 3. We abbreviate by I the 
intervals (— 00 , t] or (/, 00 ), where t G R. Further, pick 5 to satisfy 0 < 5 < (1 — 29)/v. 
Now 

P(^ max r](z, „; /) G I) 

<P( max^ry(z,,;/)G/, 2-i|P|^ ^ K 2|P|^ r{£,) > \V\^) 

+P(|£:^| < 2-^\V\^)+F{\£,\ > 2|y|^) +P(r(£:,) ^ \V\^) 

= : p{I) +p^^^. 

Using the continuity of exponential distribution, similarly as in the proof of Theo¬ 
rem 3.1(ii)-(iv) with R = 0 and 0 < 0 < 1/2, we obtain that ^ |y|-const £qj. gQjjjg 
const > 0. Thus, to show the assertion of Lemma 5.2, we need to check that 

pil) < P( max^77(2?,;/) G /) + (5.1) 

for a fixed (nonrandom) subset U := {x,: 1 ^ / < K} C Z‘' such that r{V) > lUl*^. 
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Let be the sum over all subsets V' C V with the properties: \\V\^ < 

\V'\ < 2\V\^ and r{V') > We denote by ^ the sum over all permutations 
xi,... of the subset V'. Write Then, for V D Vg, 

V" M 

Vixi)^vix 2 )> Ly, max p{x) < Ly] 

X^V\V' / 

<EE ^ -^>''^(^ 1 ) ^ ’ 7 ( 3 ^ 2 ) ^ ■ • ■ ^ p{X^y,^) ^ Ly, 

max{? 7 (x): x G + C/) U V')} < Ly^, (5.2) 

where V' + U denotes the algebraic sum of the subsets V and U. Since all the random 
variables are mutually independent, the double sum on the right-hand side of (5.2) is 
equal to 

^ma>^?7(a:i; 7) G E^^ _ p^)l'^l-(|t^l+i)l^ I 

^ ^ v 

< v(^m^^^r]{xi;l) G /^ -|- const (V D Vq), 

since |V^^|Pv ^ | via the notation. This completes the proof of (5.1). Lemma 5.2 

is proved. □ 

We now turn to the case ^(-) = f{p {-)). 

Theorem 5.3. (Cf. (2.10)). Fix y G Z^\{0} and iG G N. Then 

limsup \^{Zj^y -I-y) < 00 in probability . 

V 


Proof. The assertion follows from Lemma 5.1(i). □ 

To the end of this section, let us fix constants 0 < e < 6* < 1/2. With Ly g as in 
Section 3, we write f{x) := f,{x) if f(x) < Lyg, and f{x) := 0, otherwise. For natural 
K ^ 1 and I > K, we put 


Xk,v (0 ■=^K,V (^^K,v + - fiy 


(5.3) 


1x1=1 


here Fn' := Fo/(u;) denotes the indicator of Ll' C fl. To study the asymptotic 
behavior of variables (5.3), we introduce the class S'Gin^ of functions / := Q*~ such 
that 


lim inf 

s-100 ae{c, 9 s) 


(^fisf{fis + a)-f{s)) 


/(2q) ^ 


for any 0 < c < 1 and some 0 < ^ < 1/2. 


= 00 


(5.4) 
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The class S'AII^ is a strict subset of (4-1). The following theorem provides 

some generalization of limit (4.4) for p = 2. 

Theorem 5.4. (Cf. (2.11)). Fix K G N. If f belongs to the elasses (5.4) and 

P /<2 (4.5) with s and 9 as above, then 

lim min (l)=oo in probability. 


Proof. We begin with estimating x(0 ■= Xk,v(0 for iX + 1 < 1 < M; M > Mq. 
Let G IF denote the subset of configurations satisfying the following three 

inequalities: 

max max f{z, ^ + x) ^ f{2\ogM), (5.5) 

^K+i,v > 0 and ^ ^ const' < 1 (5.6) 

^K,V ^K,V 

for const' > const (e) specified in Theorem 4.5. According to Lemma 5.1(ii) and 
Theorem 4.5, we obtain that limsup —>■ 0 as M —>■ oo. On the other 

hand, expanding the sum X]|£i;|=i (5-3) over powers of y +a:)/Cx,v with AT+l ^ 

I < M, we get that, for any M > Mq and any V D Vo{M), the inequalities (5.5) 
and (5.6) imply the following estimate 

^ - const /(2logM), 

where const > 0 does not depend on V and M. From this implication and Theo¬ 
rem 4.3(i) with p = 2, we obtain that, for any C > 0, 

limsup P(^ nni^^x(l) < C) < limsup —>• 0 (5.7) 

as M —>• oo. 

It only remains to estimate xiO ■= Xk,v{^) for M < 1 < 2|y|^. Using formu¬ 
las (1.23) and (1.24), we represent y and f{z, y + x) in the form: 

C,,v =/(^log-^+Pv(0^ and ^{z,y+x) = f{p{z,,,+x)), 

where Py{1) := — log(|U|/l). Denote by G T the subset of configurations 

= /(’7^‘^H’)) satisfying (5.6) and the following three inequalities: 

\Pv{K)\if\\ogM, max \p^{l)\<l (5.8) 

and 

V{zi.v+x) .3 

max max -;-;- < -. (5.9) 

|a;| = l log / 2 
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We then apply Lemma 4.1(ii) to Py (1) and Lemma 5.1(iii) to the left-hand side of (5.9) 
to obtain that limsup y —>• 0 as M —>• oo. We now write := ^ V 0 and 

note that, for any M > Mq and any V D Vo{M), inequalities (5.6), (5.8) and (5.9) 
imply the following estimate: 


min y(1) 


^ min 


^ min 


^K,vi^K,v - ^i,v)-const ^ ^+{z^y+x) 

\x\=l 


/ log 


y/M 


,/ log 


M 

y/M 


-/( log^ + 1 


- const'/( 


(5.10) 


Using this implication combined with the fact that, by condition (5.4), the right-hand 
side of (5.10) tends to infinity as \V\ —>■ oo, we obtain that, for any C > 0, 


limsup P f min x(l) < c\ ^ limsup P(n\f2y^^) —0 


as M —>■ oo. This limit and (5.7) yield the assertion of Theorem 5.4. 


□ 


6. POISSON LIMIT THEOREMS FOR THE LARGEST EIGENVALUES 

This section is to provide an overview of current results on the extreme value theory 
for the spectrum of the Anderson Hamiltonian "Hy = kAv + Cvi t with an 
i.i.d. potential /(•). The results under consideration are taken from (Astrauskas and 
Molchanov 1992), (Astrauskas 2007; 2008; 2012; 2013), (Bishop and Wehr 2012), (Ger- 
minet and Klopp 2013), (Gartner and Molchanov 1998) and (Biskup and Konig 2016). 

In Section 6.1, we give Poisson limit theorems for the largest eigenvalues and the 
corresponding localization centers, provided the distribution tails e~^ of /(O) are 
heavier than the double exponential function (Theorems 6.2 and 6.9). These limit 
theorems are then complemented and illustrated by the distributions with polynomi- 
ally decaying tails, Weibull distributions and those with fractional double exponential 
tails (resp.. Examples 6.11, 6.12 and 6.13). 

In Section 6.2, we first give the second order expansion formulas for the largest 
eigenvalues, provided the tails e~^ are lighter than the double exponential function 
(Theorem 6.14). For bounded /(O), further extensions of this result are discussed. 

Section 6.3 provides the second order expansion formulas for the principal eigenvalue 
in the case of double exponential tails (Theorem 6.19), which are further extended up 
to Poisson limit theorems for eigenvalues. 

In Section 6.4, we comment and compare the proofs of Poisson limit theorems stated 
in Sections 6.1 and 6.3. We mention, en passant, that Theorems 6.2, 6.9, 6.14 and 
6.19 simply follow from the corresponding results of Sections 2-5. 


42 



6.1. Distribution tails heavier than the double exponential function 

The extreme value theory for i.i.d. random variables ^(x) deals with the asymptotic 
behavior of the K%\v largest values y of the sample \ W. It is well known 

that for suitable normalizing constants ay > 0 and hy, the non-trivial limiting (max- 
stable) distributions G{-) for P((Ci,v — by)ay ^-) are either Weibull law Dp{t) := 
exp { —(—t)^} {t < 0) or Frechet law Gp{t) = exp{— {t > 0) for some /3 > 0, or 
Gumbel law Gexp(0 = exp{—e“‘} (—oo < t < oo); see, e.g., (Resnick 1987). Note 
that the weak convergence of maxima to Gumbel law is equivalent to the limit 

lim |R|P(^(0) > -I-t/oy) = e“‘ for all t G R. (6.1) 

On the other hand, limit (6.1) implies that the point process Afy on [—1/2; 1/2]'^ x R, 
defined by 

■^v'=J2^^vG) ^v{z) ■= -by)ay), (6.2) 

zGV 

converges weakly (as V f Z"^) to the Poisson process on [—1/2; 1/2]“^ x R with the in¬ 
tensity measure dxx e“‘dt, i.e., the product of Lebesgue measure on [—1/2; 1/2]'^ and 
that defined by the increasing function logGexp(-) on R; see (Leadbetter et al. 1983). 

The necessary and sufficient conditions for (6.1) to hold are generally formulated 
in terms of F-variation of the function at the right endpoint tq or, equivalently, 
in terms of Il-variation of its inverse /olog (Resnick 1987). We say that / := Q'^ is 
in the class Till, if there exists a function a : (—oo; tq) —>■ R+ such that 

lim = c for any c G R+. (6.3) 

a(/(s)) 

Here a(-) is called an auxiliary function. The class AH (6.3) is an argument-additive 
version of the original class of H-varying functions /o log considered, e.g., in (Resnick 
1987, Section 0.4.3). In Lemma A.l of Appendix A, we recall the well-known charac¬ 
terization of the class AH in terms of Q. 

Lemma 6.1. (Resnick 1987, Sections 0.4.3 and 1.1). Limit (6.1) holds true if and 
only if f £ AH for some an auxiliary function a(-). In this ease, the normalizing 
constants can be chosen by = /(log |R|) and Oy = l/a{by). 

We now formulate Poisson limit theorems for the largest eigenvalues y of the 
random Schrodinger operator Tly = KAy+^y introduced in Section 2.1. Throughout 
this subsection, we assume that e-^ G RlGo, so that log Q(t) = o(t) as t —>■ oo by 
Lemma A.3 with p = oo. This class of distributions includes Weibull distributions 
(1.5) for arbitrary a > 0 and those with fractional double exponential tails (1.8) with 
7 < 1. Using the notation from Section 4.3, for fixed small 0 < 0 < 1/2, we write 
/(x) := /(x) if /(x) < Ly g := /((I — 0)log|U|), and /(x) := 0, otherwise. For any 
z € V, denote by A(z) the principal eigenvalue of the “single peak” Hamiltonian KAy-|- 
^(z)S^ -|-5y(l —d^). Let Ak,v be the ATth order statistics of the stationary random field 
A(-) in V, and let ZT(K).y G V stand for its location defined by A{zt{k).v) ■= ^k,v- 
(Recall that the sites z^ y G V (1 ^ I ^ |U|) are associated with the variational 
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series (1.1) based on ^y) Note that, for Z := Zt(k),v and if G N fixed, the eigenvalues 
X{Z) are expanded into a certain power series in the variables ^{Z) and ^{x) (x G V); 
cf. (2.20)-(2.22). 

Theorem 6.2. (see Theorem 4 in (Astrauskas 2007) and Theorem 5.2 in (As- 
trauskas 2008)). Let tg = oo and G RVao, and assume that e~^ is log-Holder 
continuous of order /i > (l + 0)i^/(l —20) at infinity for some small 6 > Q as above, 
i.e., (4.7) holds true. Then the following assertions (I)-(II) hold true: 

(I) (Poisson limit theorem) Assume, additionally, that there exist the normalizing 
constants > 0 and By such that 

lim|y|P(A(0) > By Ay^t) = e“‘ for any t G M, (6.4) 

and define the point process Afy on [—1/2; 1/2]^ x'M. by 

Lz \ 

where A^ (fc) (A, ^ j ; (6.5) 


then Afy converges weakly to the Poisson process Af on [—1/2; 1/2]'^ x K. with the 
intensity measure da; x e“‘dt. 

(11) (Exponential localization) Fix a constant e such that 0 < s < 9 and write 
My := log(L^,^ — Lyg) (so that My —>■ oo); then with probability one 


limsup max max —-^^ — 1. 


( 6 . 6 ) 


Sketch of the proof of Theorem 6.2(1). Using Theorem 4.6 with p = oo. Theorems 3.1 
and 4.8 with R = 0 and O<£<0<l/2, we obtain that almost surely satisfies 
the assumptions of Theorem 2.3, where if G N is fixed and N := [|U|^]. Theorem 2.3 
implies that almost surely Xk,v = Xk,v + 0 (exp{—as |U| —>■ oo, for 
fixed if G N. This asymptotic formula in turn yields that the point process Afy 
is approximated by the corresponding point process Afy where A^ y are replaced by 
Xk,v (1 ^ ^ |U|); see the proof of Theorem 4 in (Astrauskas 2007). The weak 

convergence of Afy to Af is shown by checking Leadbetter’s mixing conditions for 
the random field A(-) (Astrauskas 2007, Lemma 6 ). This concludes the proof of the 
theorem. □ 

In Corollaries 6 .3-6.7 below, we give the alternative conditions on Q (where log Q{t) = 
o{t)) for the Poisson convergence of the largest eigenvalues to hold. 

Corollary 6.3. {Specification of the normalizing constants Ay > 0 and By in (6.4) 
for some examples of potential distributions). Let Q{t) = for < ^ 0 where a > 0 
(Weibull distribution), or Q{t) = e*^ for t ^ to; 0 < 7 < 1 (fractional double 
exponential distribution). Consequently, Q satisfies the regularity and continuity 
conditions of Theorem 6.2. Moreover, the equations for the normalizing constants 
Ay > 0 and By in (6.4) are derived by applying a certain iteration scheme for Ai_y 
as in (2.20)-(2.22) combined with Laplace’s method for the corresponding integrals 
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(Astrauskas 2008, Section 6), (Astrauskas 2016); see also (Astrauskas 2013, Section 3). 
From these equations one derives the explicit expansion formulas for and hence 
those for the eigenvalues ^ up to the random max-stable fluctuations of order 
O(Ay^); cf. Examples 6.12 and 6.13 below. 

Corollary 6.4. {Specification of the normalizing constants Ay > 0 and By in (6.4) 
under general RV conditions on potential distributions). Let tg = oo and, for some 
large to, assume that Q : [to;oo) —t IR-i- is (locally) absolutely continuous with the 
positive density Q' : [to; oo) —>■ obeying the following conditions: 

lim ^ ^ = 1 for any C > 0, (6.7) 

t^oo Q'{t) ^ ' 


and 


liminf Q'{t) > 0. (6.8) 

t —¥00 

Consequently, the function Q satisfies the regularity and continuity conditions of 
Theorem 6.2. Then the centralizing constants By in (6.4) are defined by the equation: 

P(A(0) > = |E|-i (VDVo), (6.9) 

and the normalizing constants Ay > 0 in (6.4) and ay > 0,by in (6.1) are specified 
as follows: 


Ay = ay := Q'{by) where by := /(log|E|) {V D Vq). (6.10) 

Therefore, for any t gR and \V\ -G oo, 

|E|P(C(0) >6v + ^) ^ e-* and |E|p(a( 0) > -b ^ e"*, (6.11) 

where 


By = by + o(l). 


( 6 . 12 ) 


Corollary 6.4 ensures that both the distribution functions P(^(0) ^ t) = 1 — e 
and P(A(0) ^ t) {t G M) are in the domain of attraction of the max-stable Gumbel law 
Gexp(’); cf. Lemma 6.1 and the assertions before this lemma. Note also that the addi¬ 
tional condition (6.8) is to exclude the heavy-tailed (“subexponential”) distributions 
1 — e~^ which are considered in Theorem 6.9(C0) and Example 6.11 below. 

The proof of Corollary 6.4. We observe from Lemmas A.4(1) and A.3 with p = 
oo that conditions (6.7) and (6.8) imply the regularity and continuity conditions 
of Theorem 6.2. Further, from Lemma A.4(111) and Lemma A.l we derive that 
/ G An (6.3) with the auxiliary function a(-) = 1/Q'{-) in [to;oo); therefore, by 
Lemma 6.1 we obtain the first limit in (6.11). 
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To prove the second limit in (6.11), we first notice that, for each V D Vq, there is 
a solution By of equation (6.9) because of the continuity of the distribution function 
of A(0). Let us show (6.12). Since A(0) > C(0); get from (6.9) that 

|y|-i = P(A(0) > B^) > P(^(0) >B^)= 

therefore, By ^ by = /(log|fo|) for V D Vq. If ^(0) ^ Ly ^ := /((I — e)log|fo|) for 
some 0 < e < 9, then we get from (2.20)-(2.22) that almost surely A(0) < ^(0) + Py 
for some (nonrandom) 0 < /ly 0 as |fo| —>■ oo. Thus, for V D Vq, 

|fo|-i = P(A(0) >By)= P(A(0) > By,m > LyJ 

^ P(^(0) > By- py) = e-Q^Bv-Hv)^ 

therefore. By by + Py = by + o(l). These estimates imply (6.12), as claimed. 
To prove the second limit in (6.11), we also need the following observations. First, 
since liminf y ay > 0, it follows that limsup y |fo|P(^(0) p by + M) — >• 0 as M — >• oo. 
Second, for any M p Mq and any V D Vo{M), if ^(0) < by—M, then A(0) < by—M/2. 
These two assertions imply that, for any t G K., 

\V\¥{m > By +tap^) 

= |fo|P(A(0) >By+ tap\ |e(0) - by\ <M)+ Oy^{l), (6.13) 

where o^, „(l) —>■ 0 letting first V ^ If and then M ^ oo. Thus, it suffices to check 
that, for any t G K., any M p M^lt) and V '[ I'^ , 

P(A(0) >By+ tap\ |C(0) - by\ < M) 

= e-*P(A(0) > |e(0) - by\ < M)(l + 0(1)). (6.14) 

To prove (6.14), we follow the arguments of the paper (Biskup and Konig 2016, 
Section 7.1) which are now simplified and adapted to our case log(5(t) = o{t). (Re¬ 
call that this paper considers the case of double exponential tails of potential, i.e., 
log(5(t) « t/p.) The main idea here is the observation that the shift of the eigenvalue 
A(0) by is achieved by the corresponding shift of the single ^y-peak ^(0) on the 
left-hand side of (6.14). Indeed, write := ^(0) — and denote by A^*^ the 

principal eigenvalue of the Hamiltonian AcAy -l-Cy(l — 6g) in P{V). Note that 

= ^(0) and A^°) = A(0). Fix t ^ 0. Comparing expansion formulas (2.20)“(2.22) 
for A^*) with those for A(0), we find that, for any (small) e > 0, any M ^ MQ(t,e) 
and any V D Vo{M,t,s), 

if 1^(0) — by \ < M, then < A(0) — tap^ ^ A^‘\ 

therefore, we obtain the following bounds for the left-hand side of (6.14): 
p(Ab+<^) > -by\< M/2) 

^ P(A(0) >By+ tap\ |C(0) - by\ < M) (6.15) 

^ P(A(‘^ > -by\< 2M). 
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Since e > 0 is arbitrarily small, it suffices to prove limit (6.14) for the upper and 
lower bounds in (6.15). We write to emphasize the dependence of 

on the sample = {^{x)}xev- Let us consider the functions A(‘)(sy) of sy = 
{s(a;)}a;Gy G and the corresponding integrals on the right hand-side of (6.15) 
with respect to the probability measure d^"(s(a:)); here F := 1 — e”*^ stands 

for the distribution function of ^(0) with the density p{-) := F'{-) in [to;c!o). By the 
change of variables u( 0 ) := s( 0 ) — tay^ and u{x) := s{x) for all a; G F \ { 0 }, we get 
that A^*)(sy) = X^°\uv)', therefore, for fixed M > 0 and V f 

P(A(‘) >F^,|C^‘)-5y| <M) 


11 


'|V| 


I A(*) (sv)>Bv,|s(0) — —bv\<M^ 


n <^pisix)) 


xev 


f IT p(u(0) + tUy^) TT ip/ / \\ 

JrIVI {A('»(uv)>Bv.|«(0)-bv|<M} p(u(0)) 


= e/^ 11 p(g(0)+tay^) 

{A(“)(«v)>Bv.|«(0)-bv|<M} p(^(0)) 

= P(A(°)(ey) >By,|^(0)-5y| <M)(e-* + o(l)) 

by applying Lemma A.4(1V) to the ratio of the densities in the last expectation (6.16) 
where ay = Q'(bv)- Formula (6.16) combined with (6.15) implies (6.14) for t ^ 0 , as 
claimed. Since the case t < 0 is treated similarly, this concludes the proof of assertions 
of Corollary 6.4. □ 

Corollary 6.5. {Suppression of the \og-Hdlder continuity of e~'^ in Theorem 6.2). 
Let tg = oo. Assume that G RVoo, and let assumption (6.4) be fulfilled with 

Ay = 0{{Ly^ — ) for some constants 0<£<£'<(5<i. (6.17) 

Then the point process J\fy (6.5) converges weakly to the Poisson process J\f as in The¬ 
orem 6 . 2 ( 1 ) and, moreover, the ATth eigenfunction obeys exponential localization ( 6 . 6 ) 
in probability. 


(6.16) 


To prove Corollary 6.5, we again apply Theorem 2.3. So we need to show that 
the samples ^y satisfy limits (2.9), (2.14) and (2.15) in probability with the same 
abbreviation as in the proof of Theorem 6.2. First, since the condition e^ G RVoo 
implies (3.1) (see Lemma A.3 with p = oo), we may apply Theorem 3.1(iv) with R = 0 
to obtain limit (2.9) with N = [|P|®] where e' < 9 < 6. Second, by Theorem 4.6 with 
p = oo, the condition G RVoo yields (2.14) with fixed K £ N and N as above. It 
remains to prove limit (2.15) in probability with those K and N. As mentioned in 
the proof of Theorem 6.2, assumption (6.4) implies that the point process My, based 
on the sample Ay, converges weakly to the corresponding Poisson process. This 
convergence in turn yields that with probability 1 -I- o(l) the normalized spacings 
^vi^k.v ~ v) bounded away from zero as P f ^^7 for any fixed k £ N 
(Astrauskas 2007, Corollary l(jj)). Combining this with the upper bound (6.17) for 
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Av and observing from Lemma 4.2(iii) and Theorem 3.1(iv) (i? = 0) that almost 
surely “ ^n.v ^ “ ^v.e' ^ foj. ajjy V D Vo, we arrive at 

limit (2.15) in probability with N = [|y|^], as claimed. The assertions of Corollary 6.5 
are proved. □ 

Corollary 6.6. {The second order expansion formula for the top eigenvalues). Let 
again tg = oo, and G RVoo- Assume, in addition, that f{s + logs) — /(s) —>■ 0 as 
s —>■ oo. Then, for fixed K ^ 1, with probability one 

1™ (Vv - /(log Il^l)) = 1™ {^K.v - /(log |l^l)) = 0. 


Notice that the additional condition of Corollary 6.6 is to exclude the heavy-tailed 
(“subexponential”) distributions 1 — e~^, or in other words, the class of i.i.d. poten¬ 
tials whose extremes possess sharp random fluctuations; see Lemma 4.4(iii). 

The proof of Corollary 6.6. We first obtain from Theorem 3.1(iv) with R = 0 and 
Theorem 4.6 with p = oo that the samples fv (V t ^‘^) satisfy almost surely the 
conditions of Theorem 2.1(ii) and Remark 2.4 with N := [|C|^] for some 0 < d < 
1/2. Thus, using the lower bound for ^ (Theorem 2.1(ii)) and the almost sure 
limit (2.23) for (Remark 2.4) combined with Theorem 4.4 (iii), we arrive at the 
assertion of Corollary 6.6. □ 

Corollary 6.7. {Localization centers). Assume that Q satisfies the conditions of 
Theorem 6.2 on the regular increase (i.e. e-^ G RVoo) and the log-Holder continu¬ 
ity at infinity. Then almost surely the eigenfunction ■i/(-;A^ ^,) is asymptotically 
delta-function at the site Zt.(k).v G ^ for each K = o(|C|®); see Theorem 6.2(11). 
Consequently, any site in (6.5) can alternatively be defined as a localization 

center of the eigenfunction if )-; A;. ^,), viz. 

V’(Acc),v;\,v) := for some r(fc) = Ty(fc), (6.18) 

for all 1 ^ fc ^ |C|. □ 

The latter definition of the sites ZT.(k),v in (6.5) is more natural in the context of 
the localization theory for the Anderson Hamiltonians. 

The asymptotic behavior of the localization indices t{K) = tv{K) is studied by 
Astrauskas (2013). 

Lemma 6.8. (Astrauskas 2013, Theorem 2.1). Assume that the condition of Theo¬ 
rem 6.2 on the log-Holder continuity of e~^ at infinity holds true. Fix K gN. 

(i) If f G OAlP (4.3), then 

limsup Ty{K) < oo in probability. 

V 

(ii) If f G AHq (4.2) and e^ G RVoo, then 

log T {K\ 

\imTy{K) = oo and inn — = 0 probability. 
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Recall that, in Theorem 6.2 and Lemma 6.8, the condition £ RVoo implies 
log(5(t) = o{t) (see Lemma A.3 with p = oo ); the condition / £ OAU^ yields 
Q(t) X (see Lemma A.8(ii) with p = 2); finally, the condition / £ ALIq implies 
t~^Q(t) —>■ oo as t —>■ oo (see Lemma A.7 with p = 2). 

In the case Q{t) = o(t^) as f —>■ oo (for example, Weibull distribution (1.5) with 
a < 3), the eigenvalue Xk,v approaches the Kth extreme value of as R f for 
fixed K ^ 1. For such distributions, we obtain a simplified version of Poisson limit 
theorems for the largest eigenvalues: 

Theorem 6.9. (see Theorem 5 in (Astrauskas 2007) and Theorem 2.5 in (As- 
trauskas 2012)). Let tq = oo, and f £ AH (6.3) for some auxiliary function a(-), and 
assume that either of the following conditions (C0)-(C2) holds true: 

(CO) lim^^oo a(s) = oo, 

(Cl) lims_>oo sa{s) = oo and f G PI <2 (4.5) 
or 

(C2) / G (5.4) and f £ F/< 2 . (4.5) 

Write now := /(log|C|), = Uy := l/a(6^,) and 


R, := 


bv, under conditions (CO) or (Cl), 

bv + 2vK^by^, under condition (C2). 


Define the point process Afy on [—1/2; 1/2]'^ xW by 


\v\ 

A/f) := where Ay{k) 

k=l 



(6.19) 


Then Afy converges weakly to the Poisson process Af on [—1/2; 1/2]“^ x M with the 
intensity measure da; x e“*dt. 

Sketch of the proof Conditions (6.3) and (CO) imply that / G (4.1). There¬ 
fore, by Theorem 4.3(i) with p = 0, the samples satisfy the condition of The¬ 
orem 2.2(i), consequently, Xk,v = fK,v + o(l) in probability as \V\ -£ oo, for 
Hxed K £ N. Similarly, (6.3) and (Cl) imply that / is in the classes (4.1) 

and P /<2 (4.5). Therefore, by combining Theorem 4.3(i) for p = 1, Theorems 4.5 
and 3.1(iv) with R = 0 and 0 < £ < 0 < 1/2, we obtain that the samples fv satisfy the 
conditions of Theorem 2.2(ii) with N = [|C|^]. Consequently, Xk,v = + 0(^)/V) 

in probability, for fixed AT G N. Using these asymptotic expansion formulas for Xk,v, 
we obtain that in the cases (CO) and (Cl) the point process Afy (6.19) is approxi¬ 
mated by the corresponding point process Afy with ^ instead of A; see the proof of 
Theorem 2.5(ii) in (Astrauskas 2012). Since Afy converges weakly to Af (Leadbetter 
et al. 1983), this concludes the proof of the theorem for (CO) and (Cl). 

In the case of (C2), we combine Theorem 4.3(i) for p = 2 and Theorems 4.5, 
3.1(iv), 5.3 and 5.4 for i? = 0 and 0 < £ < 0 < 1/2, to find that the samples fv 
satisfy the conditions of Theorem 2.2(iii) in probability, with N and K as above. 
Consequently, Xk,v = v ^ probability, where y is the Kth extreme 

value of the i.i.d. field ^°(-) := ^(•) -I- V 1) in V. Using this limit and 
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applying the same arguments as above with ^ replaced by we obtain the assertion 
of the theorem in the case (C2). □ 

From Lemmas A.6 and A.13 we know that condition (CO) (resp., (Cl) or (C2)) 
implies that Q{t) = o{t) (resp., Q{t) = o(t^) or Q{t) = o(t^)) as t tends to infinity. 

The following corollary provides some limiting distributions for the top eigenvalues 
and the corresponding localization centers, which immediately follow from the Poisson 
convergence results of Theorems 6.2 and 6.9; see (Astrauskas 2007). We also refer the 
reader to (Leadbetter et al. 1983, Chapter 5) for a detailed survey on Poisson limit 
theorems and their applications concerning the extremal properties of random fields. 

Corollary 6.10. {Eigenvalue statisties). Assume that the conditions of Theorem 6.2 
or Theorem 6.9 are fulfilled, with the normalizing constants > 0 and specified 
herein. Then, for fixed AT ^ 1 and V f Z‘', we have the following assertions: 

(i) The normalized spectral gaps 

(•^l.v ~ ^2,V )-^V 1 ■ ■ ■ 1 ~ ^K,V^-^V 1 i^K,V ~ ^v)^V 

are asymptotically mutually independent and have limiting joint distributions with 
the density 

exp {-ti- {K - l)tK-i - Ktx - e"*^ } 

for alHfc ^ 0 (1 ^ fc ^ AT — 1) and all G K.. 

(ii) The normalized localization centers 

A(1),V A(2),V ^t{K),V 

|yp>’ |yp>’■ ■ ■ ’ 

are asymptotically mutually independent, and each of them is asymptotically uni¬ 
formly distributed on [—1/2; 1/2]“^. 

(iii) As a consequence of (ii), the distance between the localization centers is of 

order i.e., for all 1 < / < fc ^ AT, 

- z.m.v I ^ in probability. 


□ 


We now give three examples of distributions 1 — e ^, where log(5(t) = o{t). 

Example 6.11. (Astrauskas 2012). Polynomially decaying distributions. For some 
/3 > 0, assume that / o log G RVijp or, equivalently, e'5 g RVp. The latter is the 
sufficient and necessary condition for the distribution 1 — e~‘^ to be in the domain 
of attraction of the max-stable Frechet law Gp{-), or equivalently, the following limit 
holds true: 

lim |C|P(C(0) > </(log |C|)) = for all t G 1R+; 

see, e.g., (de Haan and Ferreira 2006, Chapter 1). Since / G (4-1), from 

Theorem 4.3(i) with p = 0 and Theorem 2.2(i) we see that \k,v = ^k,v + o(l) 
in probability, for fixed AT G N. Using this limit and denoting By = 0 and ay = 
l//(log|U|), we obtain similarly as in the proof of Theorem 6.9(C0) that the point 
process Afy (6.19) converges weakly to the Poisson process on [—1/2; 1/2]“^ x K.+ with 
the intensity measure /3da; x t~^~^dt. 
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Example 6.12. (Grenkova et al. 1990; Astrauskas and Molchanov 1992; Astrauskas 
2008). Weihull distributions. Let Q{t) = 1“ for t ^ 0, where a > 0. For a ^ 1, 
the function Q satishes conditions (6.7) and (6.8) of Corollary 6.4. For a < 3, the 
inverse function /(s) := Q'*“(s) = (s ^ 0) satisfies conditions (6.3) and (C2) of 
Theorem 6.9. Therefore, by Theorems 6.2, 6.9 and Corollary 6.4, 

lim |fo|P(c(0) >bv + ^) = lim |y|p(A(0) > +-^^ = e"* {t G K). 

Consequently, the point processes Afy (6.2) and Afy (6.5) converge weakly to the same 
Poisson process Af as in Theorem 6.2, where the normalizing constants can be chosen 
as follows: by = (log |fo|)^^“, Ay = ay = Q'{by) = and 

(a) By = by li a < 2, 

(b) By =by + 2vK^by^ if 2 ^ a < 3 
and, as |14| —>■ oo, 

(c) By = by + 2vK^by^ + 0{by°‘~^ ) if a > 3. 

For a > 3, asymptotic equations for By are given in (Astrauskas 2008, Section 6). 

In the case a < 1, we obtain the following almost sure asymptotic bounds for the 
eigenvalues and their spacings for any fixed AT G N and m G N\{1}: 


limsup 

V 


I^K.v - K\' 
l0g2 |l"l 


1 

K' 


and 


liminf 

V 


^ log*|fo|j j log„|fo| = -l 


limsup 

V 



m— 1 


K 


^logjfol /log^lfol 


1 


with ay and by as above; here log^ := log,f,^_i(log) for m ^ 2. For any a > 0, 
these strong limits for ^ instead of A are proved in (Astrauskas 2006, Section 3). 
Therefore, the case of A and a < 1 is derived by the same arguments as in the proof 
of Theorem 6.9(C1) above, where one explores Theorem 4.4(i) {p = 1) instead of 
Theorem 4.3(i). 


Example 6.13. (Astrauskas 2013; 2016). Distributions with fractional double 
exponential tails. Let Q{t) = e*^ for t ^ to, where 0 < 7 < 1. Obviously, Q satisfies 
conditions (6.7) and (6.8) of Corollary 6.4, therefore, 

lim |fo|P(^(0) >by + -^^ = lim |y|p(A(0) > By +-^^ = e"* {t G R). 

Consequently, the point processes Afy (6.2) and Afy (6.5) converge weakly to the same 
Poisson process Af as in Theorem 6.2; here 

by = (loglog|y|)^/^. Ay =ay = Q'{by) = yfe^'Mog |fo| 


and 


By = by + Cl 


b'^ 

Oy 


-1 


log by 


+ C2 


by Moglogfey 


(log 5y)2 


C3 


4 


-1 


(log by)' 


r(l + o(l)) 
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as \V\ —>■ oo, where ci := — 7 ) C 2 := 01(7 — 1) ^ and C 3 := C 2 log ■ 

The last formula and the asymptotic equations for By are derived in (Astrauskas 
2016); see also (Astrauskas 2013, Section 3). 


6.2. Distribution tails lighter than the double exponential function 

Throughout this subsection, we assume that the upper tails e~^ are lighter than the 
double exponential function. This class of distributions includes fractional double 
exponential tails ( 1 . 8 ) with 7 > 1 and bounded tails [tq < 00 ). 

We start with the second order asymptotic expansion formula for the largest eigen¬ 
values of T-Ly = K/S.y + ^y. 

Theorem 6.14. Let G RVq, and fix K gN. Then with probability 1 

lim(A^ - / (log |y D) = 2vn. 


Proof. We apply part (i) of Theorem 3.1, where i? G N is fixed, dfl(-) = 6* is a 
constant, and m := Thus, with probability one there is a ball ^^{zy) such that 
^(•) ^ Lyfi in M}i{zy) for some 9 G and each V D Vo(u;;i?). Therefore, 

by Theorem 4.4(iii) and Lemma A.3(ii) with p = 0, almost surely 

?(•) - ^ Lv,0 - Ly,o + 0 ( 1 ) = 0 ( 1 ) uniformly in Bfl:(zy), 

as \V\ -G 00 , for any i? G N. The latter means that with probability one the samples 
fy satisfy the condition of Theorem 2.6, therefore, y = Ly^ + 2 uk + o(l), as 
claimed. Theorem 6.14 is proved. □ 

From the proof of Theorems 2.6 and 6.14 we see that the top eigenvalue X^v of 
T-Ly is approximated by the corresponding eigenvalue of the Hamiltonian restricted to 
the (random) relevant regions Bopt := "^Rvizy) C V where ^(-) is close to and 
the diameter of which tends to infinity as \V\ —>■ 00 . 

Bishop and Wehr (2012) have obtained more accurate asymptotic bounds for the 
principal eigenvalue A^ ^ of the one-dimensional Schrodinger operators in H C Z, with 
the Bernoulli i.i.d. potential. In particular, their results imply the following 

Theorem 6.15. (Bishop and Wehr 2012). Let iz = 1, and suppose that the random 
sequence ^(•) has a common Bernoulli distribution: a = IP(C(0) = 1) and 1 — a = 
IP(^(0) =0), so that tq = 1. Then with probability one 

= fg + 2 k-I-( log |H|)“^(-£>-I- 0 ( 1 )) as VfL, 
where D = D{K,a) > 0 is the universal constant depending on k and a. □ 

In the proof of this theorem, the authors have established that almost surely the 
relevant region Bopt C H is the longest consecutive sequence of sites in V with ^(•) 
equal to 1, so that the size of Bopt is of order log|H|. See, e.g., the review paper 
by Binswanger and Embrechts (1994) for the strong and weak limit theorems for the 
length |Bopt| as V f Z. 
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Recently, Germinet and Klopp (2013) have proved the Poisson limit theorem for 
the top eigenvalues under nonlinear renormalization, i.e., for the so-called unfolded 
eigenvalues. Write, as above, N{X) (A G R) for the integrated density of states, i.e., 
the nonrandom distribution function of eigenvalues defined as the almost sure limit 
of the empirical distribution function iVy(A) := #{k: ^ ^ ^}/l^l as \V\ —>■ oo 

(Kirsch 2008). 

Theorem 6.16. (Germinet and Klopp 2013, Theorem 2.3). Assume that TLv is 
the one-dimensional Anderson Hamiltonian (v = \), where the potential is bounded 
(1C(0)I ^ const ) with the distribution density p{-) := Assume, in addi¬ 

tion, that the density p{-) is bounded and does not decay too fast at tq (say, p{t) = 
orp{t) = ^ astftq, forO < A < 1/2). Define the point process 

Aiy on the positive half-axis IR+ by 


\v\ 

■^v := 

fc=i 

Then AAy converges weakly to the Poisson process on 1R+ with the intensity measure 
dt, the Lebesgue measure. □ 

For V 2, this Poisson limit theorem was shown to hold if the Laplacian A is 
replaced by some translation invariant operator T'4>{x) = ~ 2/)) where 

T(-) is a real nonrandom function decaying exponentially at infinity. 

The proof of Theorem 6.16 relies on the improved versions of Wegner and Minami 
estimates that control the structure of eigenvalues A^ in a small neighborhood I of 
the upper spectral edge, so the amount N{I) is allowed to be exponentially small in 
see Section 1.3 for more explanations. See also (Minami 2007) for a detailed 
background of Poisson convergence results for unfolded spectral values. It is important 
for applications that this convergence result is given in terms of the integrated density 
of states, the main quantity in the theory of random Schrodinger operators. However, 
in the proof of Theorem 6.16, neither extreme value theory, nor links to the asymptotic 
geometric properties of random potential are explored. 

Corollary 6.17. Assume that Hv is the one-dimensional Anderson Hamiltonian 
{v = 1), where the potential has the distribution density p(-) satisfying the conditions 
of Theorem 6.16. Then, for fixed AT ^ 1 and V f h, we have the following asymptotic 
formulas in probability: 

(i) If p{t) = as 11 then 

Vv =tQ + 2K- (log |H|)"2+°(^). 

(ii) If there is 0 < -d < 1/2 such that p{f) = as t f tg, then 

=tQ+2n- (log |;z|)-2/(i+2^)+o(i)_ 

□ 

Corollary 6.17 can be proved by using the limit theorems for the unfolded eigen¬ 
values and asymptotic expansion formulas for the tails 1 — N{-) at the upper spectral 
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edge derived, e.g., in (Klopp 1998; Biskup and Konig 2001). See also (Klopp 2000) 
and Section 3.5 of (Kirsch and Metzger 2007) for a detailed discussion on the edge 
asymptotics of the integrated density of states. Thus, the bifurcations in the asymp¬ 
totic behavior of the top eigenvalues are caught by those in the tail behavior of the 
integrated density of states, i.e. “Lifshits tails”. Notice that the asymptotics in Corol¬ 
lary 6.17(i) agree with the results in the Bernoulli case (Theorem 6.15). Meanwhile, 
if the tails e~^ decays faster at tg, then the fluctuations of y are much sharper 
(Corollary 6.17(ii)) than those in the Bernoulli case (Theorem 6.15). 

For the tails lighter than the double exponential function including the case tg < oo, 
some heuristics on the asymptotic formulas for Aj ^ 1) and their relations to the 
long-time asymptotic formulas for the parabolic Anderson model have earlier been 
discussed by Biskup and Konig (2001), and van der Hofstad et al. (2006). Their 
assumptions are given in terms of scaling and regularity properties of the cumulant 
generating function logEe‘^*^°^ as t —>■ oo. 

Recently, Biskup et al. (2014) have proved the “homogenized” versions of limit 
theorems for the largest eigenvalues of the (scaled) finite-volume discrete Schrodinger 
operators with a bounded random potential. Their results assert that, as the 
scale parameter £ tends to zero, then 1) the largest eigenvalues of converge in 
probability to the corresponding eigenvalues of the limiting (nonrandom) finite-volume 
continuous Schrodinger operator, and 2) the fluctuations of the largest eigenvalues 
centered by their means are Gaussian in limit. 

We notice that the conditions of the above statements imply that log Q{t) ^ oo 
as t 'f tq > 0 (see Lemma A.3 with p = 0 below). 

We end this subsection with a discussion on the following important model of 
spatially continuous random Schrodinger operators: 

Example 6.18. {Schrodinger operators in with a bounded Poisson potential of 
obstaeles). Let be the ^^-dimensional continuum Laplacian. Define the random 

potential f{-) by 


f^{x) ={x-\-Xi) (cc G R"^); (6.20) 

I 

here {xi} is a Poisson point process in R*^ with the constant intensity p > 0; W{-) 
is a fixed nonnegative compactly supported, bounded measurable function, W{-) is 
non-identically zero Lebesgue-a.e. The potential f{-) is known as a Poisson field of 
“soft” obstacles. Denote by V := [—s; C R"^ the cubes of the volume \V\ such that 
K t Let us consider the principal Dirichlet eigenvalue A“°"* < 0 of the operator 
Acont -g ^(.) in V. The eigenvalue satisfies the following asymptotic formula 
(e.g., Sznitman 1998): As V f R*^, almost surely 

A“r = {\og\V\)-^/''{-Cin,p)+oil)), (6.21) 

where C{v,p) > 0 is the universal constant depending on v and p\ see below. 

Let us sketeh a derivation of the lower bound for A™(1*. By the monotonicity prop¬ 
erty of eigenvalues, we have the bound A“”* ^ A® where A° ^ is the local principal 
Dirichlet eigenvalue of the operator A“°"* + C(’) restricted to the obstacle-free con¬ 
nected open region A G V. We now maximize the eigenvalue A^ ^ over such A. 
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First, since the probability of region A to have no points {xi} is equal to 
and the number of disjoint shifts of A is of order \V\, we obtain from the Borel- 
Cantelli lemma that the volume |A| should be approximately equal to ^“^log|F|; 
cf. the proof of Theorem 3.1 above. On the other hand, we have from the Faber- 
Krahn inequality that the principal Dirichlet eigenvalues ^ of the Laplacian 
in regions A of the constant volume achieve their maximum at the ball. Thus, 
^cont ^ where A% is the principal Dirichlet eigenvalue of the 

operator in the ball Bopt := Miiy{zv) centered at some random zy & V with 

the radius Ry '■= (|Bi log |y and |Bi| is the volume of the unit ball 

Bi C Consequently, as C almost surely 

(1 + 0(1)) = Rv (A?.b, + 0 ( 1 )) = (log n) + 0 ( 1 )), 

where C{v,^) := — A° (|Bi|/i)^/^ > 0. This lower bound can be shown to be equal 
to the upper bound for A“"*, concluding the proof of (6.21). Notice also that the 
rough upper bound for A“"* can be derived by using the spatially continuous version 
of Lemma 2.8 above. 

Summarizing, we conclude that the principal Dirichlet eigenvalue A(™* is approxi¬ 
mated, as C t by the local principal Dirichlet eigenvalue A^.g ^ of the operator 
restricted to the relevant region Bopt := Bi^^,( 2 ;y) C V, so that A(°))* O Bopt. These 
observations and formulas agree with the corresponding formulas for the discrete An¬ 
derson models in Z with the Bernoulli i.i.d. potential (Theorem 6.15). 

As already mentioned, formula (6.21) and the more explicit asymptotic bounds for 
the principal Dirichlet eigenvalue A“”* {V f R^) were proved by Sznitman (1998) 
exploring his original method of enlargement of obstacles. By this method, the ge¬ 
ometry of the spatial regions where ^(-j > 0 and ^(•) = 0 is reduced to the simpler 
geometry of regions associated with the modified potential in the spectral problems, 
without changing the eigenvalues very much. Finally, notice that the asymptotic 
bounds for the principal eigenvalues are crucial for study of the intermittent behavior 
of a Brownian motion in a Poisson field of obstacles; cf. Section 7 below. □ 


6.3. The double exponential tails 

In the double exponential case (1.14), Gartner and Molchanov (1998, Theorem 2.16) 
have obtained the second order expansion formula for the principal eigenvalue Ai^v 
of Ry = kAv -|- ^y by claiming a continuity of Q. We now provide their result with 
the continuity condition removed. 

Theorem 6.19. If G RVp for some 0 < p < oo, then with probability 1 
lim(Ai,v - / (log IV^D) = 2vKq{p/K), 

where the nonrandom function q is defined in Section 2.5. 

Proof. We check the conditions of Theorem 2.7. First, by Theorem 4.7, almost 
surely fy satisfies condition (2.29). To prove (2.30), we fix constants i? G N, d > 0, 
and write 

Oiy) ■■= Oniy) := 1 - exp | (hoptiv) - /p} iv G ®fl)> 
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where the nonrandom function h®pt(-) is defined in Section 2.5. Consequently, 0{-) 
satisfies the assumptions of Theorem 3.1(i). Combining the statements of Theo¬ 
rem 3.1(i), Lemma A.3(ii) and Theorem 4.4(iii), we obtain the following assertion 
with probability one: for any V D Vb(a;; S, R) there is zy & V such that 

?(y) ^ Lvfi(y-z^) > - zv) - 2(5 for all y S Mr{zv). 

Since i? G N and (5 > 0 are arbitrary constants, this estimate concludes the proof of 
the almost sure limit (2.30). Now, Theorems 2.7 and 4.4(iii) imply the assertion of 
Theorem 6.19. 

□ 

From the proof of Theorems 2.7 and 6.19 we see that almost surely the eigenvalue 
Ai,y approaches (as \V\ —>■ oo) the local principal eigenvalue in the random region, 
where C(-) «CiT + /ioA(-) for R arbitrarily large, so that Ai,y is associated with the 
(random) “relevant island” of high ^y-values of optimal shape, the diameter of which is 
asymptotically bounded. From Theorem 3.1(iv), Remark 3.5 and the last assertion of 
Lemma A.3, it follows that the “islands” of ^y-extremes are located asymptotically far 
away from each other. Moreover, if the constant p /k is large enough, these “islands” 
are located in the neighborhood of single extremely high ^y-peaks; see (Astrauskas 
2008, Theorem 4.4 and Corollary 4.5) and (Astrauskas 2013, Theorem 2.1(iii)). 

For arbitrary 0 < p < oo, Poisson limit theorems for the largest eigenvalues and the 
corresponding localization centers are proved by Biskup and Konig (2016); see also the 
survey by Konig (2016) on these limit theorems and related topics. To formulate their 
result, we again define the sites ZT(fe),y G P by (6.18), i.e., the localization centers of 
the fcth eigenfunctions ■(/'(•; ^k,v) (1 ^ ^ ^ |P|) of the Hamiltonian "Hy = «:Ay 4- ^y. 

Theorem 6.20. (Biskup and Konig 2016, Theorem 1.2). Let tq = oo, and assume 

that Q is a continuously differentiable function such that 

V Q'W 1 , 

lim , = — tor some 0 < p < oo. 

t^oo Q{t) p 

Then the following assertions (I) and (II) hold true: 

(I) (Poisson limit theorem) There are constants By = /(log |P|) -|-2y«:(7(p/K) -|-o(l) 
and Ay = p“^log|P| such that the point process Afy (6.5) converges weakly to the 
Poisson process M on [—1/2; 1/2]^ x R with the intensity measure da; x e~*dt. 

(II) (Exponential localization) AsV f IT and K ^ 1 fixed, we have with probability 
1 -f o(l) that there exist non-random constants C > 0, M > 0, C > 0 and M' > 0 
such that 

^ c exp{-M\x - for all x G V, 

and 

Ak,v)I <C''exp{-M'(loglog|P|)|a;-z^(^,_^|} for ja; - | ^logjPj, 

i.e., the Kth eigenfunction is highly concentrated in the neighborhood of its localization 
center. 

For sufficiently large p, i.e. p > po, Poisson limit theorems and localization theorems 
for the top eigenvalues were earlier proved by Astrauskas (2007; 2008; 2013). For 
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p > Po, the corresponding localization properties present an interesting intermediate 
case between the single site concentration property, i.e. Xk,v ^ ^t(k)v’ 
case p = oo (Theorem 6.2(11)) and the non-single site concentration property, i.e. 
Xk,v ^ for 0 ^ p < Po (Section 6.2 and Theorem 6.20). 

From Remark A.5(i) and Lemma A.3, we notice that the conditions of Theo¬ 
rem 6.20 imply e-f S RVp, i.e., the assumption of Theorem 6.19. Moreover, from 
Remark A.5(ii), Lemma A.l and Lemma 6.1 with a(-) = 1/Q'{-) we also see that the 
conditions of Theorem 6.20 yield the limit (6.1) with by = /(log |R|) and ay = Ay = 
p~^ log \ V\. Consequently, the distribution 1 — e~^ is in the domain of attraction of 
the max-stable Gumbel law Gexp(’)- We finally notice that the conditions of Theo¬ 
rem 6.19 or Theorem 6.20 imply the limit /(s) = (p -I- o(l)) log s as s —>■ oo, which in 
turn is equivalent to log(5(t) = p~^t + o{t) as t —>■ oo; see Lemma A.3. 

6.4. Some comments on the proofs 

In this section, we briefly comment and compare the proof of Theorems 2.2, 2.3, 6.2, 
and 6.9 by Astrauskas and Molchanov (1992) and Astrauskas (2007; 2008; 2012; 2013) 
(“relevant single peak”approximation) and the proof of Theorem 6.20 by Biskup and 
Konig (2016) (“relevant island”approximation). 

(RSP) “Relevant single peak”approximation. As already mentioned in Section 1.2, 
the proof of Theorems 2.2, 2.3, 6.2, and 6.9 is based on the finite-rank perturbation 
arguments and the analysis of Green functions involving the cluster expansion over 
paths. To be more precise, fix Z := Zt(^k),v G tho localization center of the 
Kih. eigenfunction, and denote by A(Z) the extreme eigenvalue of Ry = nAy + ^y 
associated with the site Z. Let Gy RA; •, •) be the Green function of the Hamiltonian 
KAy + {l — 6z)^v on R{V). Under the conditions of Theorem 2.2 or 2.3 (i.e. sparseness 
and difference in height of ^y-peaks as U f the eigenvalue X{Z) is a solution to 
the dispersion equation 


g{f\X; Z,Z) 


1 

W) 


( 6 . 22 ) 


(7\ 

and the corresponding eigenfunction is Gy {X{Z); ■, Z). By expanding the Green 
function G\f'^ (A; •, •) over paths, one proves that equation (6.22) is approximated by 
the corresponding equation GviX;Z,Z) = 1/^(Z) for the principal eigenvalue of the 
“single peak”Hamiltonian nAy + ^y + ^{Z)5z', here Gv{X] •, •) stands for the Green 
function of the operator nAy + ^y. Again expanding l/y (A; •, •) over paths, one finds 
that the eigenvalue X{Z) of Ry is approximated by a certain (nonlinear) function 
on ^y and ^(Z); cf. (2.20)-(2.22). Moreover, because of the sparseness of ^y-peaks, 
the extreme eigenvalues X{Z) become asymptotically independent, so that they obey 
asymptotic Poisson behavior a.s V 'f ’E" (see Theorems 6.2 and 6.9). 

We notice that the analysis of the Green functions combined with the finite-rank 
perturbation theory is essential to study the largest eigenvalues of the finite-volume 
operators Ry in the “relevant single peak” approximation. Recall that these tech¬ 
niques also play a crucial role in the proof of the Anderson localization for the infinite- 
volume Hamiltonian R (Kirsch 2008; Stolz 2011); see also Section 1.3 above. 
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(RI) “Relevant island ’’approximation. Recently, Biskup and Konig (2016) have 
developed novel arguments to prove Poisson limit theorems for the largest eigenval¬ 
ues in the case of double exponential tails (see Theorem 6.20 above). As in the 
single-peak approximation, the analysis of the extreme eigenvalues is here based 
on controlling the dependence of an eigenvalue on the geometric properties of 
peaks and the associated regions in V. This enables to identify the “relevant” regions 
:= MRy{zy) C V (where ^(•) is high and of the optimal shape) such that the 
Kih. largest eigenvalue Xk,v of k/Xv-\-S,v is approximated by the local principal eigen¬ 
value Aj of the Hamiltonian restricted to /^(B^j) (cf. also Theorems 2.7 and 6.19 
and their proofs in the present survey). In other words, the eigenvalues associated 
with a block of “relevant islands ”of high ,Jy-values can be determined by the local 
principal eigenvalues associated with separate “relevant islands ”. It is worth noticing 
that the conditions of Theorem 6.20 imply that the islands of high ^y-values are lo¬ 
cated extremely far from each other as H (See also Theorem 3.1 for the related 
limits under the continuity assumption (3.1)). The proof of Theorem 6.20 involves 
the following procedures on a simplihcation of potential configurations: 1) those re¬ 
gions, where the potential possesses the lower values, are deleted from V {domain 
truncation and component trimming)-, 2) for the radius Ry tending to infinity slowly, 
the analysis of the local principal eigenvalues in all balls {z) C H is reduced to 
the consideration of independent identically distributed local principal eigenvalues in 
disjoint balls in V {coupling to i.i.d. variables)-, 3) the local principal eigenvalue in 
the region Bfpt is separated from other local eigenvalues in B,(^^ {reduction to one 
eigenvalue per component)-, 4) extremal type limit theorems for the local principal 
eigenvalues in {z) are comparable to each other for the different increase rate of 
Rv —>■ oo, with the same normalizing constants Ay and By {stability with respect to 
partition side), and so on. 

Summarizing, the main idea of the proof of Theorem 6.20 explores the straight¬ 
forward geometric arguments controlling the dependence of eigenvalues on potential 
configurations, rather than the techniques of resolvents or Green functions. This is in 
contrast to the relevant single peak approximation in (RSP), where the Green func¬ 
tions are the main object of analysis. On the other hand, although most of the proof 
of Theorem 6.20 is based on deterministic arguments, we are not able to reformulate 
this assertion in terms of ^y-extremes (like in Theorems 2.2-2.7 above), except for 
the case of sufficiently large p considered in (Astrauskas 2008, Theorem B.3). 


7. APPLICATIONS TO THE PARABOLIC ANDERSON MODEL 
7.1. The parabolic Anderson model 

The parabolic Anderson model (PAM) is the Cauchy problem for the following heat 
equation with random potential: 

—^—^- = k'S^ {u{s,x+y) — u{s,x))+^{x)u{s,x), s ^ 0, x&lA, 

l^i (7.1) 

u{0,x) = 5q{x), x&lA-, 
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here, as above, ^(•) is an i.i.d. random field (potential) with distribution P(^(0) > t) = 
g-Q(t). jg |.j^g Kronecker delta function at the origin (i.e., the localized initial datum 
of the problem); the variable s ^ 0 is referred to as a time. The equation has almost 
surely a unique nonnegative solution, provided ^(0) V 0 has a finite moment of order 
> ly (Gartner and Molchanov 1990). 

The PAM appears in the context of population dynamics, chemical kinetics, mag¬ 
netism and turbulence, etc. (e.g., Gartner and Molchanov 1990; Molchanov 1994). 
The following interpretation of the solution u is well-known in the mathematical liter¬ 
ature (e.g., Molchanov 1994): Let {X (s): s ^ 0) be a continuous-time random walk in 
with a generator «:Adif, where AdifV’(a^) := Let ^+{-) ^ 0 

and ^“(•) ^ 0 be independent random i.i.d. fields on Z", and write ^(•) := C*"(■)“■?” (’)• 
For a fixed realization (C"''(')i'?”(■))) consider a system of particles which obey the 
following diffusion and branching mechanism: 

1 ) at time s = 0, there is a single particle at the origin; 

2 ) particles move independently of each other according to the random walk A(-); 

3) at the site x a particle disappears with intensity ^~{x) and splits into two new 

particles with intensity ^“'"(a:), which further move according to A(-). 

Then, for a fixed realization the solution u{s, x) to (7.1) is the expected 

number of particles at the site x at time s, where the expectation is taken over 
a branching mechanism and diffusion (but not over random medium $"''(■)) C~(’))' 
Thus, the sum U{s) := u{s,x) is the expected total mass of particles at time 

s. We see from the Feynman-Kac formula (7.2) below that U{s) is equal to U{s,0), 
where U{s,-) is the solution to equation (7.1) with the homogeneous initial datum 
{7(0, •) = 1 instead of the localized one. 

For i.i.d. random potentials, the PAM exhibits an intermittency effect: As s —>■ oo, 
the overwhelming contribution to the total mass U{s) = '^^u{s,x) of the solution 
u to (7.1) comes from a small number of spatially separated and relatively small 
islands of large u{s, •)-values, i.e., intermittent islands. This is in contrast to the case 
of constant potential ^(•) = const, for which the solution u(s, •) is spread over the 
spatial ball of radius O(y^) as s —>■ oo, i.e., diffusion effect. 

Various aspects of long-time intermittent behavior of the PAM (asymptotic ex¬ 
pansion formulas for the total mass U{s) and its statistical moments, concentration 
properties for the solutions u(s, •), etc.) have been intensively studied, during the 
last two decades, by mathematicians Molchanov, Gartner, Sznitman, Konig, Biskup, 
Mdrters, den Hollander, Sidorova, van der Hofstad, and their colleagues. See (Konig 
2016) for a recent survey on the subject and references therein. The main techni¬ 
cal tools of intermittency theory are a spectral representation and the Feynman-Kac 
formula for the solutions u, U. Recall the latter formula: 


s,x) = 


exp /y ^(A(a)) dal t6(0, A(s)) 


= 


expjy C(A(a))da|do(A:(s)) 


, s > 0, X G Z"^; (7.2) 


here the expectation E^, is taken with respect to the random walk A(-) in Z*^ as above, 
conditioned by A(0) = x. Looking at (7.2), we see that the intermittent behavior is 
determined by competition between two factors: (i) extremely large exponential fac- 
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tor in (7.2) associated with a portion of the trajectories X[0; s] := {X{a): 0 < a ^ s) 
spending much time at spatial regions where ^(•) is high, and (ii) very small proba¬ 
bilities of such trajectories in (7.2). In the model, the potential necessitates concen¬ 
tration properties of u; meanwhile, the Laplacian forces these properties to be less 
expressed. It turns out that, as the upper tails of potential distribution get heavier, 
the ^y-extremes get more pronounced as 1^ = t^(s) t (cf. Section 1.2); therefore, 
the long-time intermittent properties (in particular, mass concentration properties) 
of the PAM become stronger. In (Gartner and Molchanov 1998; Gartner and den 
Hollander 1999; Gartner et al. 2007), the emphasize has been made on the double 
exponential tails (1.14). Such tails indicate the critical situation between formation 
of widely-spaced single peaks of u(s, •) in the case of tails heavier than in (1.14) (e.g., 
Gartner et al. 2007; Konig et al. 2009; Sidorova and Twarowski 2014; Fiodorov and 
Muirhead 2014), and formation of widely-spaced extremely large “islands” of higher 
values in the behavior of u{s, ■) for the tails lighter than in (1.14) (Biskup and Konig 
2001; van der Hofstad et al. 2006). The results of these papers suggest that the opti¬ 
mal strategy of particles is to move quickly to the spatial region where the potential 
values are high and of preferred shape, and to stay here for the remaining time. 

In this section, we will focus on the representation of the solutions u and U in the 
spectral terms of the Anderson Hamiltonian Hv = ^Ay -|- ^y where V = V{s)\ If. 
In view of this representation, we will discuss some techniques of the extreme value 
theory for eigenvalues, which can be applied to study intermittent properties of the 
PAM; cf. Theorems 7.1-7.2 below. 


7.2. Asymptotic expansion formulas for the total mass 

The first result in this direction was obtained by Gartner and Molchanov (1998), who 
particularly derived the second-order asymptotic formula for the logarithm of the to¬ 
tal mass U{s) (s —>■ oo) with probability one (and, as a by-product, the corresponding 
result for the principal eigenvalue of the operator "Hy), provided the potential distri¬ 
bution satisfies mild RV conditions and has all positive exponential moments finite; 
cf. the assumptions of Theorem 7.1 below. Let us sketch the proof of the almost sure 
asymptotic formula for the PAM following the arguments of their paper. The first 
key observation is that the solutions u and U are approximated, as s —>■ oo, by their 
finite-volume analogues and Uv{s)^ respectively. I.e., and are solu¬ 
tions to the corresponding equations in V (s) with the Dirichlet boundary condition; 
here F(s) C I'^ denote cubes centered at the origin, whose size length is of order 
s(logs)“ for some constant c > 1. In particular, almost surely 

U(s) = Uy^^^{s) + o{l) and u(s,x) = (s,x)-|-o(l) as s —>■ oo (7.3) 

uniformly in x, by using the standard cut-off procedure for the solutions u and U, 
based on the following facts: 1) the overwhelming asymptotic contribution to the 
Feynman-Kac representation of u and U is given by trajectories X(-) which stay 
inside the box F(s) during the whole time interval [0; s] i.e. A'[0;s] C F(s); and 
2 ) the contribution from trajectories A(-) visiting the complement of F(s) during 
the time interval [0; s] is much smaller. On the other hand, the solution Uy with 
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V := V(s) admits the spectral representation 
|V| 

^ each s ^ 0 and x G F, (7.4) 

fe=i 

and, therefore, the total mass Uy (s) = '^v (S) has the following representation 

^v(s) = X! l)yV'(0; A, for each s ^ 0; (7.5) 

k^l 

here [ip,(f)y denotes the inner product of the functions ip and (pinl'^pV), and 1 stands 
for the function taking everywhere value 1. Recall that and ipl-y) are the kth 
eigenvalue and the corresponding eigenfunction of the Hamiltonian 'Hy = kAv + S,v- 
Moreover, the eigenfunctions are chosen to form an orthonormal basis of P{V), and 
the principal eigenfunction to be strictly positive in V. We may also assume, without 
loss of generality, that ^/’(O; A^ ^,) > 0 for all k. 

Let us prove that the main asymptotic contribution to the logarithm of Uy (s) in 
(7.5) comes from the first term associated with the principal eigenvalue. Looking at 
(7.5) and applying the Cauchy-Schwarz inequality and Parseval’s identity, we easily 
obtain the upper bound 

Uvis) < (7.6) 


for each s ^ 0 and each V. To derive the lower bound for C/(s), one needs more 
sophisticated arguments: Let Y = H(s) denote centered cubes of side length of or¬ 
der s(logs)“'^ for some c > 1, and the site zi G H as a localization center of the 
principal eigenfunction of the operator "Hy = kAv + ^v- Recall the Feynman-Kac 
representation for U{s): 


U{s) = Eo 


exp ^{X{a)) da 


s ^ 0, 


so that U{s) ^ LL{s), where C/(s) is the same expectation Eo[-] when restricted to the 
particle trajectories X[0;s], which initially move from the origin to the site zi until 
time 1, then stay in Y (s) time interval of length s— 1, at the end of which the particles 
return to zi and the remaining time move freely. Assuming that ^(•) is bounded from 
below (so percolation effects of very low values of ^(•) are neglected), we obtain from 
the strong Markov property that C/(s) = v{s — 1, Zi) almost surely, where v{s, x) 
is the expectation over trajectories in Y (s) starting from x G H(s) and ending at zi 
during the time interval [0; s], viz. 


v{s, x) = E; 



Thus, v(s,-) is the solution to equation (7.1) in H(s) with the initial datum 
instead of Sg. Using the spectral representation for v(s — l,zi) (where all terms 
are nonnegative!) combined with the previous estimates, we finally obtain the lower 
bound for U{s): As s —>■ oo, almost surely 

U(s) ^ x(s-l,zi)e°(") > e^i-^®- 2 '^'^®+°("V(- 2 i;Ai.y)^ 

> 0'^l.yS“2l/KS-|-o(s) | —1 _ gAi.yS —2i/KS-|-o(s) 
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since V'(- 2 i; Ai,y)^ > \V\ ^ by the definition of the site zi G V. From this formula and 

(7.6) , we get that almost surely 

Aj — 2vk + o(l) ^ s“^ log U(s) < Aj — 2iyK + o(l) as s —>■ oo. 

Applying the almost sure asymptotic formulas for the principal eigenvalue Aj in 
Corollary 6.6 (p = oo), Theorem 6.19 (0 < p < oo) and Theorem 6.14 (p = 0) of the 
present paper, we can now derive the corresponding asymptotics for the total mass 
Uis): 

Theorem 7.1. (Gartner and Molchanov 1998; Section 2.1). Let essinf ^(0) > —oo. 
Assume that f := Q'^ satisfies the following RV conditions: there is a constant 
0 ^ p ^ oo such that G RVp and, additionally, /(a + log a) — /(a) —>■ 0 as a ^ oo. 
Then with probability 1 

og U{s) _ ^ ^ 2i/k(1 — qlp/n)) + o(l) as s —>■ oo. ( 7 . 7 ) 

s 

Here the nonrandom constants q(p) are specified in Section 2.5; in particular, q{p) 
are strictly decreasing in p; < 7 ( 0 ) = 1 and q{oo) = 0 . 

Recall that the condition G RVp with 0 < p < 00 (resp., p = 00 and p = 0) 
ensures the double exponential upper tails (1.14) (resp., heavier and lighter upper tails 
than the double exponential) of the potential distribution 1 — e”*^; see Lemma A.3 
of Appendix A. The additional RV condition of the theorem is to exclude heavy¬ 
tailed distributions of potential. Therefore, the first term on the right-hand side of 

(7.7) is equal (with the accuracy o(l)) to the largest values of the potential in V(s); 

see Theorem 4.4(iii). The second term describes the shape of the potential in the 
neighborhood of its maxima and is specified by deterministic variational principles; 
see Section 2.5. From Sections 1.2, 2 and 6 , we know that the principal eigenvalue 
Ai v{b) of the operator Hvis) = is approximated (as s —>■ 00 ) by the local 

principal eigenvalue in the connected region Aopt(s) C R(s) where the potential fv(s) 
possesses high values of a particular preferred shape. The logarithmic asymptotics of 
U (s) is therefore fully specified by these high values of the potential. 

Gartner and Molchanov (1998) derived also the second-order expansion formulas 
for statistical moments of the total mass U{s) as s —>■ 00 . For the upper distribu¬ 
tional tails of ^(0) lighter than the double exponential, Biskup and Kdnig (2001), van 
der Hofstad et al. (2006) obtained more accurate expansion formulas for statistical 
moments and almost sure behavior of C/(s). The spatial correlation structure for the 
PAM was investigated by Gartner and den Hollander (1999). See also (Molchanov and 
Zhang 2012) for the Anderson parabolic model with Adif replaced by the fractional 
Laplacian —(—Adif)^, 0 < 0 < 1, where the potential has Weibull type tails. In these 
papers, refined variational arguments were involved to obtain additional information 
on intermittent islands of solutions u{s,-) and the asymptotic structure of related 
^v(s)-extremes (their size, optimal shape, etc., as s —>■ 00 ). See also (Konig 2016) for 
a recent survey on the subject. 

Van der Hofstad et al. (2008) considered the case of i.i.d. potentials with heavy 
upper tails, i.e., polynomially decaying (Pareto) distributions and Weibull distribu¬ 
tions (1.5) with a < 1. Thus, all positive exponential moments of the potential are 
infinite, in contrary to the assumptions of Theorem 7.1. For such classes of potentials. 
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they proved extremal type limit theorems and almost sure asymptotic bounds for the 
logarithm of the total mass U{s). 

7.3. Asymptotic concentration formulas 

However, the rough asymptotic expansion formulas for the PAM (like in Theorem 7.1) 
provide only appropriate information on the geometric structure of intermittent is¬ 
lands of the solutions u{s, •) to (7.1). Recall that the intermittent islands are formed 
by those highly concentrated m(s, •)-values which give the main contribution to the 
total mass U{s) = and the contribution from the complement of these 

islands is negligible as s —>■ oo. Recently, there has been a considerable attention 
to the following mathematical problems regarding a geometric characterization of 
intermittency effect: 

1) description of the shape and location of intermittent islands; 

2) description of the shape of potential values which generate intermittent islands; 

3) specification of the minimal number of these islands, etc. 

Thus, taking into account (7.3), one needs to prove the following concentration 
formula for the total mass U{s): As s —>■ oo, 

n{s) 

U{s) -- (7.8) 

in the sense of asymptotic equivalence almost surely or in probability, where ASpt(s) 
are believed to present random connected regions (i.e. intermittent islands) in V (s) 
at a large distance from each other, such that the diameter of each ASpt(s) is much 
smaller than this distance; n(s) are relatively small numbers, and R(s) t are 
centered cubes as above. 

Let us give a heuristic explanation of formula (7.8) in the spectral terms of the 
operator = nAy -|- in R = H(s), provided the conditions of Theorem 7.1 are 
fulfilled. To this end, we need more careful inspection of spectral representation for¬ 
mulas (7.4)-(7.5), by applying extreme value theory for eigenvalues including Poisson 
limit theorems and localization properties discussed in Sections 1.2, 2, and 6. First, 
notice that the exponents of the top eigenvalues are essentially larger than those asso¬ 
ciated with lower eigenvalues. Therefore, it suffices to consider the sum of a few first 
terms of (7.4)-(7.5) associated with the largest eigenvalues Aj, 1 < fc < n := n(s); 
the other terms in (7.4)-(7.5) associated with the lower eigenvalues are asymptotically 
negligible. Thus, the random field u(s, •) (s —>■ oo) may be interpreted as a superposi¬ 
tion of a few wave functions ^(- y) for 1 ^ k ^ n. The general theory of Anderson 
localization suggests that the fcth eigenfunction \jj{- ;X^ y) is exponentially localized 
at the s-dependent center y € V, and it is highly concentrated in a s-dependent 
neighborhood C R of the site . The regions C R possess a relatively 

small size and are asymptotically far from each other according to Poissonian behav¬ 
ior of the localization centers z^^^^ k ^ 1. Moreover, ^(•) possesses in A^pt 

the deterministic optimal shape specified by the variational principles; cf. Section 2.5. 
From these observations when applied to (7.4) as s —>■ oo, we see that the function 
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is a very good approximation for uv(s,-) in the region for each 1 ^ fc ^ n. 
This in turn suggests that the mass concentration formula (7.8) holds true, where 
V = V(s) f Z'' as above, n = n(s) are relatively small numbers, and A^p^. C V are the 
relevant regions defined above. Thus, the intermittent islands ASpt(s) have relatively 
small size and are far away from each other as s —>■ oo. This concludes the heuristic 
explanation of formula (7.8). 

Gartner et al. (2007) proved the concentration formula (7.8) almost surely, for 
potential distributions satisfying the conditions of Theorem 7.1 with 0 < p ^ oo, 
i.e., the double exponential upper tails and heavier than the double exponential. 
They showed that almost surely n{s) = and the connected regions ASpt(s) c 

t^(s) are asymptotically bounded only when defined properly (see above), and the 
distance between them is of order For p = oo, the regions ASpt(s) shrink to 

singletons. Moreover, the shape of the solutions u(s, ■) and potential values in Aqpj(s) 
are specified (via the variational formulas) by the local principal eigenfunction and 
the principal eigenvalue in ASpt(s). This agrees with the heuristics given above; 
cf. also Sections 1.2, 2, and 6 treating the extreme value theory for eigenvalues of 
the Hamiltonian "Hv. However, the proof of the asymptotic concentration becomes 
complicated by applying straightforwardly the asymptotic results for the top spectrum 
of "Hv, because of the possibly different signs of the kih eigenfunctions in (7.4)- 
(7.5) where the factor ■i/'(0; A^, ^^) should be also taken into account. Instead, the 
authors explore the Feynman-Kac formulas for u{s, •), U (s) as well as for the principal 
eigenfunctions with a slightly modified potential. To prove the exponential decay of 
the principal eigenfunctions, they apply a decomposition technique for the trajectories 
of the random walk in the corresponding Feynman-Kac representations. 

Sznitman (1998) earlier proved similar mass concentration results for a Brown¬ 
ian motion in among Poisson obstacles; here the potential ^(•) is given by for¬ 
mula (6.20). In particular, the spatial regions ASpt(s) C in (7.8) were shown to 
have no obstacles and unboundedly increase almost surely as s —>■ oo. The optimal 
strategy of the Brownian particle during the time period [0; s] is to move quickly to 
one of the obstacle-free regions Aqp^(s) of the optimal shape, i.e. the ball of radius 
const (logs)^A^ and to stay here for the remaining time. Notice that the intermittent 
behavior of this model is rather similar to that of the spatially discrete PAM with 
the potential bounded from above. The related asymptotic results for the principal 
Dirichlet eigenvalues are discussed in Example 6.18 of the present survey. 

However, in (Gartner et al. 2007) and (Sznitman 1998), the problem of the minimal 
number of intermittent islands was not considered. This problem in precise setting 
was solved by several mathematicians for Pareto distributions 

P(e(0) >t)= e-'3(‘) = (t ^ I) (7.9) 

with (3 > u, and Weibull distributions (1.5) with arbitrary a > 0. Recall that the 
choice oi P > V in (7.9) is to guarantee the existence and uniqueness of the solution 
u(s,-) to equation (7.1). Write, as above, U(s) := ^^u{s,x) for the total mass. 

Theorem 7.2. (e.g., Konig et al. 2009; Sidorova and Twarowski 2014; Fiodorov 
and Muirhead 2014). Assume that the potential has either Pareto distribution with 
fj > V, or Weibull distribution (1.5) with arbitrary a > 0 . Then there exists a random 
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process Zopt(s) (s ^ 0) with values in 1'' such that 


lim ^ — I probability. 

s—>00 U (sj 


(7.10) 


This theorem states the complete localization property for u{s, •) as s —>■ 00 , which 
is the strongest case of mass concentration formula (7.8) in probability with n{s) = 1 
and A;(pt(s) = {Zopt(s)}, a singleton. 

For the Weibull distribution with 0 < a < 2, asymptotic formula (7.10) was proved 
by Sidorova and Twarowski (2014). This result was extended by Fiodorov and Muir- 
head (2014) to an arbitrary a > 0. See also (Muirhead and Pymar 2014) for the proof 
of the single-site concentration property for a random walk in a random environment 
(instead of the standard random walk like in the PAM) for the Weibull-distributed 

1.1. d. potential. 

Earlier, Konig et al. (2009) proved (7.10) for the Pareto-distributed potential. They 
also established a two-site concentration property for the PAM: Almost surely U (s) ~ 
u (s,Z(i)(s)) + u (s,Z(^^(s)) as s —>■ 00 , where Z^^^s) are distinct random processes 
with values in TZ. I.e., one obtains (7.8) with probability one, where n(s) = 2, 
AQpt(s) and Aqpj(s) are two singletons. This is the strongest almost sure version of 
the localization property, since for i.i.d. random potentials, the random field u{s, •) is 
asymptotically concentrated on at least two distinct sites in IZ. The reason of this 
fact lays on the observation that the localization sites are changing infinitely often 
as s —>■ 00 ; thus, at very rare time moments when particles move from the previous 
localization site to a new one, the mass of particles should be concentrated on at least 
two different sites. 

We recall from Sections 3-5 of the present survey that, for potential distributions as 
in Theorem 7.2, the ^y-peaks are strongly pronounced as E t therefore, the single¬ 
site concentration formula, no surprise, holds true according to the general picture 
of intermittency based on the Feynman-Kac formulas. On the other hand, for such 
potential distributions, there is a very precise extreme value theory for eigenvalues of 
the operators "Hy = KAy -t^y, which can (and does) provide powerful techniques for 
the investigation of single-site concentration properties; cf. Sections 1.2, 2.2-2.3 and 

6.1. See also (Konig 2016) for a recent survey on the subject. 

Let us sketch the proof of Theorem 1.2 for the Weibull-distributed potentials by 
applying the extreme value theory for eigenvalues. We follow the terminology and 
ideas of (Fiodorov and Muirhead 2014). First, one obtains finite-volume approxi¬ 
mation formulas (7.3) with “good” accuracy, where V = E(s) are centered cubes in 
IZ with the size length of order s(logs)^/“. Then, let us look at the spectral rep¬ 
resentation formulas (7.4) and (7.5) for My(s, •) and 17y(s). Because of the factor 
(^/>(-; Aj, ^,), u(0, •))^ = V’(0;Afc_v) ^ 0 (1 ^ fc < |E|) in (7.4), we need to study the 
penalised spectrum 

T(s; k) := Afc,y -I- s“4og ^’(0; A;, ^,) - 2nK {1 k ^ \V\) 
instead of the usual Spect ("Hy) = {Aj, : 1 < fc < |E|}. Let us rewrite My (7.4) in 
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the terms of the penalised spectrum d'(s) := {^'(s; k): 1 < fc < l^|}: 


Uv{s,x) = 

fc=l 


{x€V). 


(7.11) 


For 1 ^ ^ \V\, denote by ^,(s) the ^th largest value in the sample ^'(s). It will 

turn out that the gap 'I'i_v(®) ~ 'I^2,v(®) between the the first largest ^,(s) and the 
second largest 'I'^ y{s) in tE'(s) is sufficiently large; and moreover, the eigenfunctions 
ip{-;Xk v) decay exponentially (so that ■i/'(0;A^ > 0), for each k ^ \V\^ with e > 0 

small enough; cf. also Theorem 6.2, Corollary 6.10 and Example 6.12 of the present 
paper. We shall prove that the right-hand side of (7.4)-(7.5) is dominated by just one 
term associated with y{s); therefore, the localization center of the corresponding 
eigenfunction should be the concentration site for the random field uv(s,-) as s —>■ 
oo. To be more precise, let Aopt := A^. ^ and V’opt(’) denote the eigenvalue and 
eigenfunction of the operator Hy associated with the first largest value 4/^ ^ (s) among 
the penalised spectrum 4/(5), i.e. 

^i.y(s) = Xkyv + log'i/'opt(0) - 2yK; 


here k* < \V\'^ with e > 0 small enough. Let Zopt(s) £ V stand for the localization 
center of the eigenfunction V’opt(-)- Also, for random processes V(s) and W{s), we 
write Y{s) « IF(s) as s —>■ 00, if the difference F(s) — W{s) tends to zero sufficiently 
fast in probability. Using this abbreviation and applying Cauchy-Schwarz inequality 
for eigenfunctions in (7.11), we obtain that 


max 

X 


Uv{s, x) 

gS'I'i,v(s) 


V^opt {x'j 


^ \V\ exp{-s i'i’i,v{s) - 4 ' 2 ,v(s))} « 0, 


therefore, for the total mass Uv{s) we have that 


Uv{s) 


(V'opt7 l)y 


^ |Up exp {-s (4'i,y (s) - 4'2.y(s))} 


0 , 


as s —>■ 00. Consequently, 

Uy (s, .^opt (^)) ^ V^opt (-^opt('^)) ^ ^ 

Uv{s) ('*/'opt,l)y 

as s —>■ 00, because of the sharp exponential decay of '0opt(’)- Since u(s, •) « t6y(s, •) 
in V and U{s) « Uv{s), the last formula implies (7.10), as claimed. 

However, the penalised spectral values are too complicated to handle. In order 
to study the spacings of the largest values in ^/(s) as well as further properties of 
the concentration site Zopt(s), one needs a good approximation for 4'^ y(s ) by a 
simpler function on potential configurations. To this end, we introduce the following 
auxiliary quantities: Write J := [(a — l)/2] (= the integer part) for a > 1, and J = 0 
otherwise. Given z € V, let (z) denote the principal eigenvalue of the single-peak 
Hamiltonian 

kA^. + Y + on 

v- 
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I.e., A(-')(z) is the local principal eigenvalue on the lattice ball of radius J with a single 
^y-peak at z surrounded by the island of lower ^y-values. Let := X^'^\zr{k),v) 
denote the /cth largest value of the sample {A*^'^^(z): z € V}. We now observe from 
Section 6.1 of the present survey that the eigenvalue Xk,v is approximately equal to 
the local eigenvalue A^'^y, i.e., A^^y ~ x\fy as F = y(s) t and k < \V\^ with e > 0 
small enough. This observation and more careful inspection of the exponential decay 
of the fcth eigenfunction ip{-;X^^y) at the origin (cf. Astrauskas 2008; Section 6) suggest 
that, with high probability, the Kth largest value of the penalised spectrum 

d>(s) is approximately equal to the A'th largest value Ty^y(s) of the penalisation 
functional 


Ti'^)(s,z) :=A(‘^)(z)-^-i^^^-2y«; {zeV). (7.12) 

s a 

Comparing the penalised spectral values to (7.12), we observe, in addition, that the 
quantity log logs = log 6 y(s) + 0 ( 1 ) is the nonrandom rate of the exponential 
decay of the fcth eigenfunction '4’i-;X^y) described in Proposition 1.3; moreover, 
log 7 /’( 0 ; A;. ^,) asymptotically behaves like — | log &y as C f Recall also 

that, for each z £ V, the eigenvalue A*^'^^(z) is a certain (nonlinear) function of the 
sample {^(z + x): 0 < |x| < J}; cf. (2.20)-(2.22). Now the concentration site Zopt{s) 
can be defined as the maximizer of the random field tO) (g^ •) in R; here 0 ) is a 

certain function of both the sample {^(z + x): 0 ^ |x| ^ J} and the site z. Similarly 
as in Theorem 6.2 above, one obtains Poisson limit theorems for the (normalized) 
penalisation functionals T*^'^^(s,-) and their locations. This limit theorem implies the 
limiting distributions for the normalized concentration site .^opt(s) as well as for the 
spacings T^y(s) — y(s) of the random field T('^^(s, •) in V. The latter in turn 

implies the existence of a sufficiently large gap between the largest values in the pe¬ 
nalised spectrum 'I'(s), as claimed. This assertion concludes the heuristic proof of the 
single-site concentration property (7.10) of the solution u{s, ■). □ 

As stated in (Fiodorov and Muirhead 2014), the above considerations are a starting 
point in obtaining more information on the asymptotic behavior of the concentration 
site Zopt(s) and the shape of potential in its neighborhood as s —>■ 00 , provided the 
random potential has Weibull distributions. In particular, Poisson limit theorems for 
the penalisation functional (7.12) imply that the site Zopt(s) G R(s) is of order s up 
to logarithmic corrections and has the limiting distribution as a product of univariate 
Laplace distributions. Recall also that the random field T’^'^^(s, •) has the finite range 
(=J) of dependency; see (Fiodorov and Muirhead 2014) where J is called as the 
radius of influence. For sufficiently heavy tails (Weibull distributions with a < 2), 
the eigenvalue A*^'^^(z) in (7.12) may be replaced by ^(z), so that T('^)(s, z) (z G V) 
are independent non-identically distributed random variables. (Recall that the case 
a < 2 was studied in (Sidorova and Twarowski 2014) by exploring the Feynman-Kac 
representations). These observations are crucial for describing the shape of ^(•) in the 
neighborhood of the concentration site Zopt(s): As s —>■ 00 , with probability 1 -I- o(l) 
the single peak f (Zopt(s)) is extremely high: 

f (^opt(s)) = &v(«) (1 + 0(1))) 
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where by := (log cf. Example 6.12 above. Meanwhile, the neighboring values 

^(x) (1 ^ |x — Zopt(s)| ^ J) are essentially lower: there exists a strictly decreasing 
nonrandom function d{-): [1; J] !->■ [0; 1) such that 

^ - 6d(l--Zopt(s)|) fo, 1 ^ la; _ Zopt(s)| < J. 

This characterization of the concentration site Zopt(s) agrees with the asymptotic 
results for the largest eigenvalues and eigenfunctions 'tp {-; of the Hamilto¬ 
nian Hy = kAv+^v as V f Z'' and K ^ 1 fixed, provided ^(0) has Weibull distribu¬ 
tion; cf. Sections 1.2, 2.2-2.3, and 6.1 of the present paper. In particular, let us look 
at the asymptotic expansion formulas for A^ (Example 6.12) to observe that the 
leading term by comes from an isolated high peak y); meanwhile, the further 

terms of order o(l) come from the neighboring values ^(x) (1 ^ |x — | ^ J) 

with the same influence radius J as in the PAM above; cf. also (Astrauskas 2008; 
Section 6 ). Moreover, the eigenfunction ; A^ ^,) is highly concentrated at the site 
^t{k) V! the localization center (cf. Theorem 6.2(11)); and the neighboring values ^(x) 
(1 ^ |x — I ^ J) have significant influence on the asymptotic behavior of the 

localization index t{K) = tv{K) and on the concentration degree of 'ip{- ;Xk,v) 
the neighborhood of v’ (Astrauskas 2013). These observations suggest the 
following conclusion: the lighter are the tails of potential, the larger is the influence 
radius in both models; thus, the weaker are the localization properties of both models, 
PAM and (time-independent) Anderson Hamiltonian. 

It is worth mentioning, at the heuristic level of rigor, that Section 2 (resp., Section 
6 ) of the present paper exhibits all classes of “typical” configurations of the potential 
(resp., RV classes of potential distributions) which are thought to guarantee certain 
concentration properties for the solutions u to equation (7.1). For instance, the results 
of Sections 2.3 and 6.1 (i.e., “relevant single-peak”approximation) should be prelim¬ 
inaries to establish a single-site concentration property in probability for solutions 
u, provided the upper tails of potential are heavier than the double exponential. In 
particular, the results of Section 2.2 and Theorem 6.9 (“sharp single-peak ’’approxi¬ 
mation) are related to the single-site concentration property for u with zero influence 
radius J = 0, provided the distributional tails are heavier than Weibull’s tails with 
a = 3. In view of the results of Sections 2.5 and 6.3 (“relevant island”approximation 
in the double exponential case), it can be conjectured that, with high probability, the 
solution u{s, •) exhibits entire concentration on a single island, the diameter of which 
is asymptotically bounded as s —>■ oo. See (Konig 2016) for a heuristic explanation of 
this conjecture. 

If the upper tails of potential are lighter than those of double exponential including 
bounded tails (tg < oo), no results on concentration properties for the PAM are 
known, with the exception of the very special (but important) spatially continuous 
model of Brownian motion in a Poisson field of obstacles studied by Sznitman (1998). 


A APPENDIX A: REGULAR VARIATION 

In this section, we study the classes of functions / := Q*~ (the left-continuous inverse 
of the cumulative hazard function) introduced in Sections 4-6. These classes are 
characterized in terms of Q. The tail behavior of Q{t) as t '[ tq is also treated. In 
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Section A.l, we recall the classical results on the domain of attraction of max-stable 
Gumbel law and regular variation RVp. The classes A\Y^ (4-1), Allg (4.2) and OAII^’ 
(4.3) are studied in Section A.2, and P /<2 (4.5) in Section A.3. Finally, examples 
and counterexamples are given in Section A.4. 

A.l. The domain of attraction of max-stable Gumbel distribution and regular 
variation 

We now give the well-known characterization statements for the distribution function 
to be in the domain of attraction of the max-stable Gumbel law Gexp(’)- 

Lemma A.l. (Resnick 1987; de Haan and Ferreira 2006). The following assertions 
(i)-(iii) are equivalent: 

(i) / € An (6.3) with an auxiliary function a{-) > 0; 

(ii) there exists another auxiliary function ai : {—coAq) —>■ M+ such that 

Q{t + cai{t)) — Q{t) ^ c as tftq, for any c G K; 

(iii) there exist functions b : (—oo;tQ) —>■ R and 02 : (—oo;tQ) —>■ M+ such that 

Q{t) = b{t)+ f 1 / 02 ( 5 ) ds {tKtq), 

Jto 

here b{t) —>■ 6 G R (t t ^q), the function 02 is locally absolutely continuous with the 
density 03 (t) —>-0 (tf tq) and, for tq < 00 , 02 (t) —>-0 (tftq). 

In this case, aof(s) = oio/(s)(l-|-o(l)) = a 2 o/(s)(l-|-o(l)) as s 00 . Moreover, 
the limit in (ii) with an auxiliary function ai(-) >0 implies that f G AH (6.3) with 
the same auxiliary function a {-) = ai(-). 

Example A.2. For tq = 00 , p ^ 0 and B > 0, consider the subclass AH^ C AH 
associated with the auxiliary function 02 ( 5 ) := B{p + l)“^s“^ in Lemma A.l(iii). In 
this case, Q{t) = B~^t^^^ -I-const +o(l), i.e., 1 — e~^ are Weibull type distributions. 
In the next section, we extend the subclass AH^ to the boundary cases B = 00 , B = 0 
and 0-type asymptotics. 

Let us discuss the class RVp of (nondecreasing) functions, which are regularly vary¬ 
ing at infinity with index p. Recall that, for 0 < /3 < 00 , the condition /olog G RVxjp 
is sufficient and necessary for the distribution 1 — e“'5 to be in the domain of attrac¬ 
tion of max-stable Frechet law G^(t) := exp{—(t > 0); cf. Example 6.II. We 
now explore the class RVp to characterize the double exponential type distributions. 
Lemma A.3. For tq = co and 0 ^ p ^ 00 , the following assertions are equivalent: 

(i) ef G RVp,- 

(ii) /(s) — f{6s) —>■ —p\ogS as s ^ 00 , for any 0 < 5 < 1; 

(iii) Q{t + C)/Q(t) —t as t ^ co, for any G > 0. 

For any 0 ^ p ^ 00 , either of (i)-(iii) implies that hmt_>oo log (5(t) = p~^. 
Finally, for any 0 < p ^ 00 , either of (i)-(iii) yields that Qit—)/Q{t) —t I as t —t 00 , 
i.e., the continuity condition (3.1). 

Proof. The equivalence of (i)-(iii) follows from Theorems 1.5.12, 2.4.7 and Propo¬ 
sitions 2.4.4(iv) and 1.3.6(i) in (Bingham et al. 1987) combined with the observation 
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in (Resnick 1987, Sect. 0.2) that Q{t — e) ^ ^ Q{t) for all t G K and all 

£ > 0. For p = oo and p = 0, the equivalence of (i)-(iii) is also proved, respectively, 
in Lemma A.6 {p = 0) and Lemma A.7 {p = 0) of the present paper adapted for 
the argument-additive functions; see also (de Haan and Ferreira 2006, Chapter 1 and 
Appendix B.l) for the case 0 < p < oo. □ 

The following lemma is compounded of Lemma A.l and Lemma A.3 with p = oo, 
provided that there exists the density of the distribution 1 — . 

Lemma A.4. (Cf. Corollary 6.4). Let tq = oo. For some large to, assume that 
Q : [to;oo) —>■ K+ is (locally) absolutely continuous with the positive density Q' : 
[to;oo) —>■ K+ obeying the following conditions: 

^ ^ C > 0, (A.l) 

i^oo Q'[t) 

and 


liminf Q'{t) > 0. 


(A.2) 


Then the following limits (I)-(IV) hold true: 

(I) limt_>.oo Q(t + u)/Q(t) = limt_>oo Q'(t + u)/Q'(t) = 1 uniformly in compact sets 
of u G K.; 

( II ) liminf 4 _>oo Qit)/t > 0 ; 

( III ) hmt_>oo {Q{t -I- vai{t)) — Q{t)) = v uniformly in compact sets o/u G K., with 

ai(-) = l/(5'(-) in [to;oo); 

(IV) with p{t) := e~^'d)Q'[pj (t ^ to) as the distribution density and ai(-) as in 
part ( 111 ), 


lim 

t—^OO 


p{t + U + vai(t)) 


= e 


p{t + u) 

uniformly in compact sets ofv,u G M. 

Proof. (1) By L’Hopital’s rule, we obtain the first limit for any u G 


The 


uniform convergence follows from Theorem 1.2.1 in (Bingham et al. 1987) adapted 
for the argument-additive functions. 

(II) The assertion follows from (A.2). 

(III) Writing 

Q{t + vai{t)) — Q{t) — V = V J (^ai{t)Q'{t + 9vai{t)) — d9 for t ^ to 

and applying assertion (1) and condition (A.2), we easily obtain the claimed limit. 

(IV) Let us rewrite the ratio under the limit in the form: 


p{t + U + vai{t)) 
pit + u) 


= exp I — (Q(t + U + vaift)) — Qit + m)) | x 
Q'it -I- u -I- vaiit)) /Q'it + u) 


Q'it) 


\ Q'it) 


(A.3) 
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t ^ to- Since ai(-) is a bounded function, assertion (I) implies that the last two ratios 
on the right-hand side of (A.3) converge to 1 locally uniformly in u,v G M. It remains 
to prove the uniform convergence of the exponent on the right-hand side of (A.3). By 
the theorem of continuous convergence (see, e.g., p. 2 in (Resnick 1987)), it suffices 
to check that, for arbitrary functions u{t) —>■ u and v{t) —>■ v, the following limit holds 
true: 

Q(t + u{t) + v{t)ai(t)) — Q(t + u{t)) — v(t) ^ 0 as 1 1 oo. 

This is shown similarly as in part (III), so we omit the details. Lemma A.4 is proved. □ 

Remark A.5. (Cf. Theorem 6.20). Let tq = oo. For some to, assume that Q : 
[to; oo) —>■ 1R+ is (locally) absolutely continuous with the positive density Q' satisfying 
the following condition: (log(5)'(t) p~^ as t —>■ oo, for some 0 < p < oo. Then the 

following limits hold true: 

(i) hmt_>oo Q{t + C)/Q{t) = limt_>oo Q'it + C)/Q'{t) = uniformly in compact 
sets of C G R; 

(ii) with a(-) = 1 /Q'(-) in [to;oo), 

lim ((5(t -I- Ca{t)) — Q{t)) = C for any C G K. 


The proof of the assertions of Remark A.5. (i) Write logQ in the form: 

t /■* 

logQ(t) = const -I-h / £(s)ds (t > to), 

P Jto 

where £{t) := (log(5)'(t) — p~^ —>■ 0 as t —>■ oo. Using this representation and the 
conditions of Remark A.5, we obtain the claimed limits for any (7 G M. The uniform 
convergence follows from Theorem 1.5.2 in (Bingham et al. 1987) adapted for the 
argument-additive functions. 

(ii) Since a(t) = o(l), the claimed limit is derived similarly as in the proof of 
Lemma A.4(III). □ 

A.2. Classes AU^, ABg and OAW 

Recall that, for p ^ 0, the classes AII^, AIIq and OAII^ consist of functions / := 
satisfying, respectively, /(s)^(/(s + c) — /(s)) —>■ oo, —>■ 0 and x 1 as s —>■ oo, for 
any c > 0; cf. (4.1)-(4.3). We first formulate the results of this section. To avoid 


trivialities, we restrict ourselves to the case tq = oo. 

Lemma A.6. For any p ^ 0, / G AXAf^ if and only if 


lim {Q{t + ct~^) — Q{t)) = 0 for any c > 0. 

t^OO ^ ' 

(A.4) 

In this case, 


Q(t) = o(t^"''^) as t ^ oo. 

(A.5) 
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(A.6) 


Lemma A.7. For any 0, f G AIIq if and only if 

lim (Q{t + ct~^) — Q{ffj = oo for any c > 0. 


In this case, 


Q{t)t ^ ^ oo as t ^ oo. 


(A.7) 


Lemma A.8. For any p ^ 0 and f G OAW, the following assertions hold true: 

(i) there is a constant c > 0 such that Q{t + ct~^) — Q{t) 1 as t ^ oo; 

(ii) Qft) X as t ^ oo; 

(iii) if a function a : K+ —>■ K+ is chosen to satisfy liminf s_>oo o(s) ^ ci > 0 and 
liminf s_>oo(s — a(s)) ^ C 2 > 0 , then 

const (s — < f{s) — f{a{s)) 

< const' (^s^/(P+^'> - a(s)i/(P+i) + a(s)-P/(P+i)) 

for any s ^ sq md for some const' ^ const > 0 . 

Before proving Lemmas A. 6 -A. 8 , we provide an example of Q satisfying assertion (i) 
of Lemma A .8 such that / := does not belong to OAIfP for p ^ 0, and further 
on, two technical lemmas for later use. 

Example A.9. For p ^ 0, write Q{t) := + [t], t ^ 0. Note that, for each c > 0, 

Q {t + ct~P) — Q{t) = c{p + 1) + g{t) + o{l) as t —>■ oo, 

where 0 ^ g{t) := [t+ct~P] — [t] = 0(1). I.e., Q satisfies the assertion of Lemma A. 8 (i) 
for any c > 0. However, for each t := n G N, we get Q{n) — Q{n—) = 1, therefore, 
/ := ^ OAH^ according to Lemma A.10 below. 

Lemma A.10. If liminfs_>oo /(s)^(/(s + c) — /(s)) > 0 for each c > 0 and for 
some p ^ 0, then 

lim (Q{t) - Q{t-)) = 0, 
i.e., Q is continuous at infinity. 

Proof. Assume for a moment that there exists a sequence —>■ oo such that 
Q{tn) — Q{tn—) —>■ c° > 0 . This limit implies that s„ := Q{tn) -G- oo and, in addition, 
that f{sn — c) = f{sn) for any 0 < c < c° and any n ^ no(c), contradicting the 
assumption of the lemma. This completes the proof of the claimed assertion. □ 

Lemma A.11. (Resnick 1987, pp. 4). For all s G K+ and t G (—oo;tQ), the 
following assertions hold true: 

(i) /(s) ^ t if and only if s ^ Qit); 

(ii) /(s) > t if and only if s > Q{t); 

(iii) Q(/(s)-) ^ s < Qifis)). 
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We now are in a position to prove Lemmas A.6-A.8. To simplify the proceedings, 
we need the following abbreviations: 

fp{s;c) := f{sf{f{s + c)-f{s)) and Qp{t; c) := Q{t + ct~P) - Q{t). 

Proof of Lemma A. 6. Assume first that (A.4) holds true. I.e., for each £ > 0 there 
is So = so(£) > 0 such that 

Qp{f{s);s~^)<e for all s ^ sq. (A.8 ) 

Since (A.4) implies a continuity of Q at infinity, from Lemma A.ll(iii) we have that 
Qifi^)) ^ for each s ^ sq. This and (A.8) yield that Q(/(s)+e“^/(s)“^) < s+2£: 
for each s ^ sq. Inverting Q (see Lemma A.ll(ii)), we get that fp{s]2e) > 1/e for 
each e > 0 and each s > so(£). I.e., / G 

To prove the inverse implication, assume for a moment that there are a sequence 
of reals —>■ oo and constants c > 0, £ > 0 such that Q{t„ + ct~^) > Q{tn) + £ for all 
n > no(£,c). Inverting Q (see Lemma A.ll(i)), we get that f {Q{tn) + s) ^ tn + ctf^, 
which combined with < f{Q{tn) + £/2) (see Lemma A.ll(ii)) gives that, for each 
n > no(£,c), 

fp{sn; e/2) < c with s„ := Q(t„) + e/2. 

Since s„ —>■ oo, the latter violates the assumption / G concluding the proof of 

the first part of Lemma A.6. 

To prove (A.5), we note that, for any natural M ^ 2 and any s ^ 2M, 

/(s + I)P+^ - /(s)P+^ > fp{s; I) ^ I{M) := inf /p(s; I), 

s^M 

and, therefore, 

/(sf+i - /(Mf+i > (s - M - 1)/(M). 

Hence liminf s_>oo^ /(M)^Ap+i)_ Since I{M) —>■ oo (as M —>• oo) by 
the assumption, the latter implies that /(s)s“^Ap+i) qo as s —>■ oo, which in turn 
yields (A.5). Lemma A.6 is proved. □ 

Proof of Lemma A.7. Assume first that (A.6) holds true, i.e., for each £ > 0 there 
is to = to{e) such that Qp{t; e) 1/e for each t ^ to- By Lemma A.ll(i), the latter is 
equivalent to f[Q{t) + l/e) < t+et~P. Substituting t := f{s) —>■ oo into this inequality 
and then applying Q(/(s)) ^ s (see Lemma A.ll(iii)), we obtain fp{s;l/e) < £ for 
each £ > 0 and each s ^ so(£)) i-e., / G AHq. 

To prove the inverse implication, assume for a moment that there are a sequence 
t„ —>■ oo and constants d > 0, c > 0 such that Qp{tn]5) < c for each n ^ no(c,5). 
Here, inverting Q (see Lemma A.II(i),(ii)) and denoting s„ := Q{t„) —>■ oo, we obtain 
that fp{sn] c) > 5 for each n ^ no(c, d), contradicting the assumption / G AHg. This 
completes the proof of the first part of Lemma A.7. 

We will prove (A.7) under the weaker condition by assuming (A.6) for some c > 0. 
(The forthcoming arguments are applied to prove assertion (ii) of Lemma A.8 as well). 
Write := ( 5 (c^Ap+i)i) and observe that 

lim I) = lim Qp{t;c) = oo, 

f —Von ^ f —Von 
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i.e., limit (A. 6 ) for c > 0 is reduced to that for c = 1. With the abbreviation 
Rik) := we obtain that, for fixed natural M ^ 3 and any natural t ^ 2M, 


Qk'){t)-Qk)(M) 

t-l R(k)-1 


E E 


k^M Z=0 
t—1 -R(fc) —1 


fc + 


i + i 

R{k) 


g(-)( 


fc + 


R{k) 


t-i 




k^M Z=0 


k^M 


/P +1 

= mt^O!,'>(r;l )—(1 + 0(1)) 


as t —>■ oo. Here, by (A.6), the infimum tends to infinity as M —>■ oo, therefore, (A.7) 
is fulfilled. This completes the proof of Lemma A.7. □ 

Proof of Lemma A. 8. (i) We first prove that if, for each c > 0, the function /p(s; c) 
is asymptotically bounded away from zero as s —>■ oo, then Qp{t; S) = 0(1) as t —>■ oo, 
for some <5 > 0. Assume otherwise that, for each 5 > 0, there exists a sequence 
—>■ oo such that Qp{tn',5) > 2M for any M > 0 and any n > no{M). Here, 
inverting Q similarly as in the proof of Lemma A. 6 , we obtain that fp{sn] M) ^ 6 
with Sn '■= Q{tn) + M —>■ oo as n —>■ oo. Therefore, liminfs_>oo fpis; M) ^ 6. Since 
5 > 0 is arbitrary, we obtain the contradiction proving the desired implication. We 
next observe that, if Qp{t] 6) = 0(1) for some 5 > 0, then 


Qp{t]k6) = 0(1) as t ^ oo, for any fc € N. (A.9) 


This implication is easily proved by induction in k. We omit the details. 
With the abbreviation 


M := 2 limsup /p(s; 1) > 0, 

S—^OQ 

we finally show that the function Qp(t', M) is asymptotically bounded away from zero 
as t —>■ oo. For this, fix an arbitrary sequence oo, and define a sequence {««} 

by f{sn+) ^ tn ^ f{sn) (^^ G N). Write t„ := /(««). Combining Lemmas A.10 
and A.ll(iii), we have that Q{tn) — Qijn) = o(l) and, consequently, 

Qp{tn] M)Qp{Tn,M) p o{l) as n ^ oo. (A.IO) 

On the other hand, from the definition of M, it follows that /(sn) + Mf{sn)~^ > 
/(sn + 1) for any n ^ np. Applying Q to both sides of this inequality and then using 
Lemmas A.IO and A.ll(iii), we obtain that Qp{Tn',M) ^ 1/2 for n ^ uq. The latter 
combined with (A.IO) implies that the sequence Qp{tn',M) is asymptotically bounded 
away from zero, as claimed. This and (A.9) conclude the proof of part (i). 

(ii) The assertion is shown by the same arguments as in the proof of limits (A.5) 
and (A.7). We omit the details. 

(in) If a(s) or s — a{s) are bounded from above for any large s, then the bounds in 
(iii) simply follow from part (ii) and condition (4.3). 
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For simplicity we abbreviate a := a{s), and assume that both a and s — a tend to 
infinity as s —>■ oo. By combining assumption (4.3) and the limit /(s) x we 

obtain that, for s ^ sq, 

/(s)-/(a)«; (/(a + fc + l)-/(a + fc)) 

O^k^s—a 

^ const E (a + fc)-P/(P+i) 

O^k^s—a 

ps — a 

^ const + const / {a + dk 

Jo 

and 

/(s)-/(a)> Y if{a + k + l)- fia + k)) 

O^k^s—a —1 

> const E ia + k)-P/<~P+^'> 

O^k^s—a—1 

> const {s — a — 1 )s~p/^P'^^\ 

as claimed. Lemma A .8 is proved. □ 

Remark A.12. (A relationship with classical regular variation), (i) Consider the 
casep = 0. Obviously, for /3 = oo or /3 = 0, / is in AII^ if and only if g := expo/olog G 
RV 0 . Therefore, for p = 0, Lemmas A .6 and A.7 follow from the well-known results 
for the class AVg with /3 = oo and /3 = 0, respectively (Bingham et al. 1987). 

The class OAII*^ links to the exponential type distributions 1 — e”*^, with Q{t) x f 
as t — >■ oo. Moreover, if / G OAn°, then g := expo/ o log is in ORV, the class of 
0-regularly varying functions studied, e.g., in (Bingham et al. 1987, Section 2). 

(ii) In the case of p ^ 0, if / G OAH^, then / o log is asymptotically balanced or, 
equivalently, the maximum of i.i.d. sample is stochastically compact (Bingham 
et al. 1987, Sections 3.11 and 8.13.12). 


A.3. Class P /<2 

Recall that the class P /<2 consists of functions / := such that 

f (os') 

liminf > 1 for some 1 < c < 2; 

s^oo /(s) 

cf. (4.5). This is the subclass of positive increase class PI considered, e.g., in (Bingham 
et al. 1987, Section 2.1.2). 

Lemma A.13. / G PIc 2 if and only if 

limsup Q{Cot)/Q(t) < 2 for some Cq > 1. (A.11) 

t—¥(X) 

In this case, 


\ogQ{t) 

limsup —^^ const 

t—^oo log t 


log 2 
log Co ■ 


(A.12) 
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Proof. Assume in contrary to (A. 11) that, for any C > 1, there is a sequence 
tn = tn{C) —>■ oo such that Q{Ctn) ^ (2 — e)Q{tn) for any e > 0 and any n > no(e). 
Similarly as in the proof of Lemma A.6, inverting Q (see Lemma A.ll(i),(ii)) and 
writing s„ := Q{tn) + £ —>■ oo, we obtain that 


limsup f{Ssn)/f{sn) ^ C for any 1 < 5 < 2 and any C > 1 

n 


or, equivalently, lim„/(5s„)//(s„) = 1, contradicting the assumption / G PI< 2 - 
According to these arguments, the inclusion / G P /<2 implies (A.11). 

To show the inverse implication, we suppose otherwise that / := Q'^ ^ PI< 2 , 
i.e., for each 1 < c < 2, there exists a sequence s„ = s„(c) —>■ oo such that 
f{csn) ^ (1 + £)/(sn) for any £ > 0 and any n ^ no(s). In this inequality, we 
invert / (see Lemma A.ll(i)) to obtain cs„ < (5((1 + £)/(sn)). On the other hand, 
by Lemma A. 11 (hi), 


CSn > cQ{f{Sn) “ ) > cQ {{1 - £)/(s„)). 


Summarizing these estimates and using the abbreviations := (1 — s)f{sn) and 
5 := (1 + £)/(l — £), we have that Q{5tn) ^ cQ{tn). Since 1 < c < 2 is an arbitrary 
constant but close to 2, the latter implies the limit limsup j Q{5f)/Q{t) ^ 2 for each 
5 > 1, contradicting assumption (A. 11). This concludes the first part of the lemma. 

Let us show (A.12). By (A.11), there exist numbers C > 1 and to = ^o(C') such 
that Q{Cf)/Q{f) < 2 for all t > to- Applying this estimate, we obtain that, for any 
n G N and any t G [(^^"to; C'”+^to), 


Q{t) QjC^to) 

Q(C"to) ■ Q(C"-ito) 


QjCto) 

Q{to) 


■ Qito) 


^ Q(to)2”+i ^ constt(^°s2)/^°s^, 
i.e., (A.12) is done. Lemma A.13 is proved. 


□ 


A. 4 . Comparison of the classes Ali^, ALE and P/< 2 - Examples 

In view of limit theorems for eigenvalues (see Theorems 6.2 and 6.9), we need to 
compare the classes (4-1), SAA^ (5-4), AA (6.3) and P /<2 (4.5) of functions 
f:=Q^. 

Lemma A. 14 . (i) For any p ^ 0 , there exist examples fi G PI^2\AAP^ and 

/2 G AAP^\PI^2- Consequently, there is f2 G S'An^\P/< 2 . 

(ii) There exist examples fs G P/<2\An and fi G An\P/<2 with an auxiliary 
function 04 ^ I. 

(iii) For any p > 0, there exists an example /s G nP/<2)\An and, therefore, 

there is /s G (S'AII^ fl P/< 2 )\An. 

(iv) For p ^ 0 , if f € AA with an auxiliary function a : R+ ^ R+ such that 
s^a^s) —>■ 00 as s ^ 00 , then f G AAf^. 

(v) For 0 ^ p < 00, if is in RVp, then f{s + logs) — /(s) —>■ 0 as s —>■ 00 . 

Proof, (i) Write Q{t) := [t], t ^ 0; i.e., 1 — e~^ is the geometric distribution. Let 
us show that /i := Q'^ G P/<2\An()^. Indeed, since Q{n) — Q(n—) = 1 for all n G N, 
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Lemma A. 10 implies that /i ^ for any p ^ 0. However, since Q satisfies (A. 11), 
we have that /i G PI< 2 , as claimed. 

Consider the function / 2 (s) := b(t) dt, where b(t) := n if 2^” < t ^ and 

b{t) := 2” if 2^”+^ < t < 2^”+^ for all n G NU{0}. Let us show that /2 G An5^\P/<2. 
Obviously /2 G AH)^. With s„ := 2^”+^ and 1/2 < 5 < 1, we see that 

f 2 {Ssn) fsll Ht) dt 

/2(Sn) 

where /J// 6(t) dt = const 4”n and 

rs„ ^2sj \ 

/ &(t) dt = ^ I / Idt + 2* dt I + / n dt 

■tl ;^g \Jsi/2 Jsi ) Js„/2 

n—1 „ 

= ^ (P 4' + 2 • 8') + n • 4" = I • 8”(1 + o(l)) 

1=0 


as n —>■ oo. Summarizing, we find that /2(<5sn)//2(s„) —>■ 1 (as n —>■ oo) for each 1/2 < 
(5 < 1, i.e., /2 ^ PI< 2 - Since An;(^ C S'AH^, we also obtain that /2 G S'An^\P/< 2 . 

(ii) As in part (i) above, let /a be the inverse of the cumulative hazard function of 
the geometric distribution. Since /a is not in AH (Resnick 1987, Corollary 1.6), we 
obtain that /a G P/< 2 \An. 

To prove the existence of /4 G An\P/<2, it suffices (via Lemmas A.l and A. 13) to 
find a continuous function a : [l;oo) —>■ [l;oo), with derivative a'(t) —>■ 0 as t f oo, 
such that the function Q{t) := /j*l/a(s)ds does not satisfy (A.11). For this, we 
abbreviate := (logn)”, m„ := (logn)^. 


1 f tn+l 
\ tn 


and b^ .— 1 -t- 


and consider the functions 


a{t) ■= 


i/in if P < t ^ P+1 

bn if P+1 UIti ^ ^ P+1, 


for n P 2, 


(A.13) 


and Q{t) := J* l/a(s)ds. Obviously, a ^ 1 is continuous in [l;oo) and a'{t) —>■ 0 as 
1 1 oo. Therefore, /i := Q'^ G AH with the auxiliary function a. Let us show that, 
for each c > 1, 


Qip^n) 

Q{tn) 


1 Itl" 

1 “1“ —T- —y 00 

/r l/a(i) dt 


as n —>■ oo. 


so that /4 ^ P/< 2 - Indeed, from (A.13) we see that, for n ^ no(c). 


(A.14) 




■ dt = tn log C. 


(A.15) 
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To estimate the integral in (A. 14), we again use (A. 13) and the bound a ^ 1. 
Thus, 



dt 

a(t) 


E 


k^2 



('tk+i—m-k 


'tk 


dt 

t 


r*fc+i 


1 dt 




n—1 

< X! (tk^ogitk+i/tk) +mk) < t„(logn)“^/^ 

k=2 


(A.16) 


for any n ^ no- Now (A.15) and (A.16) imply (A.14), as claimed. 

(iii) Consider the example Q{t) = t + sint (t ^ 0) given by Von Mises. Obviously, 
Q satisfies (A.4) and (A.11), consequently, := is in AII^ fl P /<2 for any 
p > 0. However, /s ^ AH. We observe that the function fe{s) := s + sins (s > 0) 
is also in (AH^ fl P/< 2 )\An for any p > 0. (This is verihed by straightforward 
calculations). Consequently, since AU^ C S'AH^, the functions /s and /e are in 
(5An^nP/<2)\An. 

(iv) The assertion follows from the definition of AH (6.3) and AH^ (4A). 

(v) The assertion follows from Lemma A.3(ii). Lemma A.14 is proved. □ 

We finally provide two examples of distributions which represent RV classes con¬ 
sidered in Sections 3-6. 


Example A.15. For a > 0, let Qa{t) = {t ^ 0), i.e., 1 — e”*^" is Weibull 
distribution. Clearly fa{s) := Q'^{s) = for s ^ 0. By straightforward calcula¬ 
tions, we obtain that if a < p -I- 1 (resp., a > p -I- 1 or a = p -|- 1), then fa £ AH^^ 
(resp., fa G AHq or /„ G OAH^). Also, for a < 3, /q, G SAH^. Finally, for any 
a > 0, fa £ PI< 2 , expo/o, £ RVoo and fa is in AH with the auxiliary function 
a{t) := The latter means that the distribution 1 — e“'5“ is in the domain 

of attraction of the max-stable Gumbel law; cf. Section 6.1. 


Example A.16. Given 7 > 0 and p > 0, let {t ^ to), i.e., the 

fractional double exponential distribution. Then fj^p{s) = (plogs)^/'*' for s ^ sq. 
Obviously, if 0 < 7 < 1 (resp., 7 > 1 or 7 = 1 ), then expo/^ p is in RVoo (resp., RVq 
or RVp); cf. Section 6 . 

Acknowledgments. The author would like to thank the referee who reads the first 
version of this paper for his insightful remarks, critical comments and suggestions. 
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